UNCLASSIFIED 


r 


Rep^uxiuced 
Ixf  ike 


ARMED  SERVICES  TECHNICAL  INFORMATKW  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTW  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  \diat3oever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  emy 
patented  invention  that  may  in  any  way  be  related 
thereto. 


AFCRL-63-238 


U1  I 


N 

CD 

lO 

O 

O 


IKSTABILITY  OF  STRATIFIED  SHEAR  FLOW 


by 


JdRCEN  HOLMBOE 

Pr*|*ct  DlMctor 


Department  of  Meteorology 
University  of  California 
Los  Angeles 


Contract  No.  AF19(604)-7999 
Project  No.  8628 
Task  No.  86280 


FINAL  REPORT 


March  1963 


Prepared 

for 


A  S  T  I  A 


Mr  r 


‘I  L:.- 


GEOPHYSICS  RESEARCH  DIRECTORATE 
AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  STATES  AIR  FORCE 
BEDFORD, MASSACHUSETTS 


AFCRL-63-238 


INSTABILITY  OF  SmTITIED  SHEAR  FLOW 


by 


Jorgen  Ilolmboe 
Project  Director 


Department  of  Meteorology 
University  of  California 
Los  Angeles 


Contract  No.  An9(6o4)-7999 
Project  No,  8628 
Task  No.  86280 


FL’WL  REPORT 
March  1963 


Prepared 

for 


GEOKIYSICS  RESEARCH  DIRECTORATE 
AIR  FORCE  CAMBRIDGE  RESEARCH  UBORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS 


SequeBts  for  additional  copies  by  Agencies  of  the  Department 
of  Defense,  their  contractors,  and  other  Government  agencies  should 
be  directed  to  the: 

ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 

Department  of  Defense  contractors  must  be  established  for  ASTIA 
services  or  have  their  'need- to -know*  certified  by  the  cognizant 
military  agency  of  their  project  or  contract. 

All  other  persons  and  organization  should  apply  to  the: 

U.S.  DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TECHNICAL  SERVICES 
WASHnraTON  25,  D.C. 


(i) 


INSTABILITY  OF  A  STRATIFIED  SHEAR  FLOW 
Introduction 

T^ls  report  contains  a  discussion  of  Instability  In  a  number  of 
hydrodynamlcal  models  which  have  variations  of  density  and  velocity 
along  the  vertical.  Most  of  the  results  which  are  presented  here  are 
not  new.  However,  an  attempt  has  been  made  to  exhibit  the  similarity 
of  the  physical  mechanism  of  the  Instability  In  the  different  models, 
which  Is  sometimes  hidden  In  the  classical  treatment  of  the  subject.  If 
the  report  has  succeeded  in  this  objective.  It  Is  at  least  In  part  because 
the  so-called  method  of  symmetric  waves  has  been  used  In  the  mathematical 
analysis  of  the  models.  This  method  studies  the  Instability  as  an  Initial 
value  problem,  instead  of  deriving  the  instability  from  the  properties  of 
of  the  normal  modes  of  the  system. 

The  report  may  be  regarded  as  a  systematic  introduction  to  the 
theory  of  long  unstable  cyclone  waves  In  atmospheric  models  which  take 
Into  account  the  compressibility  of  the  air  and  the  rotation  of  the  eai'th. 
Much  of  the  preliminary  studies  of  these  baroclinic  waves  has  been 
carried  out  under  the  sponsorship  of  this  contract,  but  is  not  yet  ready 
for  publication.  This  work  will  be  continued  and  reported  on  imder  a 
current  three  year  contract  which  is  sponsored  by  the  National  Science 
Foundation. 
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Shearing  Instability 


1.  Inatability  of  a  vortex  sheet.  -  Conaider  a  system  consisting 
of  two  unbounded  homogeneous  fluid  layers  with  the  eame  density  which  have 
the  uniform  relative  translation  2U  parallel  to  their  plane  surface  of  separ¬ 
ation.  This  surface  is  called  a  vortex  sheet.  IVe  shall  use  a  frame  of  refereii 
in  ivhich  the  fluid  layers  move  in  opposite  directions  with  the  same  speed  U, 

We  call  this  the  symmetric  frame. 

The  vortex  sheet  way  be  regarded  as  a  continuously  distributed  row 
of  vortex  filamonts  of  equal  strength  and  sense.  The  strength  of  the  fila¬ 
ments  in  terms  of  the  velocity  discontinuity  2U  is  obtained  by  applying  Stolyes' 
theorem  to  a  rcctemgular  curve  with  one  pair  of  sides  parallel  to  the  vortex 
sheet.  If  the  length  of  these  sides  is  L,  the  circulation  of  the  curve  is 
2UL,  which  is  the  combined  strengths  of  the  vortex  filaments  distributed  along 
the  length  L  of  the  vortex  sheet.  Since  the  distribution  is  uniform,  the 
vortex  strength  per  unit  length  of  the  vortex  sheet  is  2U.  This  quantity 
is  called  the  sliding  vorticity  of  the  vortex  sheet. 

It  is  readily  seen  that  the  vortex  filaments  are  stationary  in  the 
symmetric  frame.  Consider  any  one  of  them.  This  filament  is  not  moved  by 
its  own  field.  All  the  other  filaments  may  be  grouped  in  pairs  whose  neutral 
stagnation  points  coincide  with  the  filament  in  question,  so  none  of  the 
fields  move  the  filament.  The  same  is  true  for  all  the  vortex  filaments. 
Accordingly  the  undisturbed  plane  vortex  sheet  is  stationary. 

Let  us  next  assume  that  the  vortex  sheet  is  deformed  into  a  sinusoidal 
surface  of  wave  length  Ls2nA  nnd  u  small  amplitude  A^,  such  that  IcA  <1. 


Consider  now  the  motion  of  the  filaments.  The  filaments  at  the  nodal  points 
have  no  motion  because  all  the  other  filaments,  when  grouped  symmetrically 
in  pairs  with  respect  to  a  node,  have  fields  whose  neutreil  stagnation  points 
coincide  with  the  node.  Consider  next  the  filament  on  one  of  the  crests. 
Again  all  the  other  filaments  may  be  grouped  in  pairs  symmetric  with  respect 
to  the  crest.  The  field  from  each  pair  gives  the  crest  filament  a  partial 
motion  in  the  direction  of  the  flow  in  the  upper  layer.  In  the  same  way  we 
find  that  the  trough  filaments  move  in  the  direction  of  the  flow  in  the 
lower  layer.  This  means  tliat  vorticity  is  being  concentrated  at  the  down¬ 
wind  nodes  and  depleted  from  the  upwind  nodes.  It  is  readily  seen  that 
this  periodic  redistribution  from  the  initial  uniform  distribution  of  the 
vortex  filaments  results  in  the  evolution  of  a  new  field  with  ascending 
motion  at  the  crests  and  descending  motion  in  the  troughs.  So  the  amplitude 
of  the  vortex  sheet  deformation  begins  to  grow. 

The  vertical  motion  associated  with  the  periodic  vorticity  con¬ 
centrations  and  depletions  at  the  vortex  sheet  nodes  has  the  same  sense 
in  both  layers  (see  the  middle  diagram  in  fig.  l).  This  field  adds  to  the 
initial  vertical  motion  field  v/hich  is  associated  with  the  deformed  interfac 
Since  the  first  stream] incs  in  both  layers  coincide  with  the  vortex  sheet, 
the  initial  vertical  motion  is  everywhere  opposite  in  the  two  layers  with 
maximum  amplitude  at  the  nodes.  The  resultant  of  this  initial  field  and 
the  evolving  field  from  the  periodic  vorticity  concentrations  at  the  nodes 
is  a  field  which  moves  downwind  in  both  layers  in  the  symmetric  frame. 

During  this  evolution  the  deformation  of  the  vortex  sheet  has  no  progressive 
motion  in  the  symmetric  frame,  so  at  all  times  the  flow  in  both  layers  is 
parallel  to  the  vortex  sheet  at  the  nodes.  This  means  that  the  amplitude 
of  the  vertical  motion  must  grow  faster  than  the  amplitude  of  the  vortex 
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sheet  as  long  as  the  progressive  motion  of  the  wave  continues. 

The  wave  must  stop  before  the  upper  and  the  lower  field  arrive 

in  phase.  For  the  streamlines  in  both  layers  are  parallel  to  the  vortex 

sheet  at  the  nodes.  If  the  upper  and  lower  wave  are  placed  in  the  same 

phase  (lower  diagram  in  fig.  1),  the  upper  and  lower  streamlineb  are  both 

horizontal  at  the  points  where  the  vortex  sheet  deformation  should  have  its 

nodes.  So  the  vortex  sheet  can  have  no  deformation  in  this  state.  Since 

the  vortex  sheet  amplitude  must  keep  on  growing  after  the  initial  state  with 

the  upper  and  the  lower  field  in  opposite  phase  (top  diagram  in  fig.  1), 

it  is  evident  that  the  wave  cannot  reach  the  state  of  equal  phase  which 

calls  for  no  deformation  of  the  interface.  The  wave  must  become  stationary 

in  some  intermediate  state  (middle  diagram  in  fig.  1).  When  this  state 

of  stationary  phase  is  reached* bo th  the  wave  amplitude  and  the  vortex  sheet 

amplitude  will  keep  on  growing  at  the  same  relative  rate,  as  long  as  the 

simple  sinusoidal  periodicity  is  not  destroyed,  that  is  as  long  as  the 

product  kA  remains  much  less  than  unity.  During  this  early  part  of  the 
s 

evolution  the  growth  of  the  vortex  sheet  eunplitude  is  equal  to  the  vertical 
velocity  at  the  crest  of  the  vortex  sheet  deformation.  This  vertical 
velocity  is  proportional  to  the  vortex  sheet  sunplitude,  so  A ^A^«  const. 

In  other  words  during  the  early  linear  part  of  the  evolution  the  wave  grows 
at  an  exponential  rate  in  the  state  of  stationary  phase. 

Consider  the  evolution  of  the  wave  from  the  state  of  equal  phase 
(bottom  diagram  in  fig.  1).  Ihis  state  has  added  vortlcities  distributed 
sinusoidally  along  the  xmdeformed  vortex  sheet,  with  the  maximum  positive 
and  negative  vorticity  added  to  the  vortex  filaments  shown  in  the  figure. 

The  vertical  component  of  this  added  field  has  the  same  phase  in  both 
layers.  It  deforms  the  vortex  sheet  into  a  sinusoidal  surface  with  nodes 
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at  the  points  of  maximum  added  vortlclty.  As  this  deformation  develops, 
it  follows  from  the  symmetry  that  it  cauinot  have  any  progressive  motion 
in  the  symmetric  frame.  So  the  streamlines  are  parallel  to  the  vortex 
sheet  at  the  nodes  at  ail  times.  But  this  kinematic  constraint  forces  both 
the  upper  and  the  lower  field  to  move  symmetrically  upwind  from  the  initial 
non-deformed  state.  Again  it  is  clear  that  the  wave  in  this  evolution  canjio 
reach  the  state  of  opposite  phase  since  it  moves  downwind  through  that  state 

The  two  non-tilting  states  shown  at  the  bottom  and  the  top  in  fig.  1 
are  accordingly  mutually  exclusive.  In  the  evolution  from  either  one  of 
these  states  the  wave  cannot  reach  the  other.  From  either  non-tilting 
state  the  wave  approaches  the  same  state  of  stationary  phase  which  has  a 
downwind  tilt  from  the  non-tilting  state  with  the  deformed  interface  and 
an  upwind  tilt  from  the  non-tilting  state  with  the  undeformed  interface. 

The  wave  in  the  vortex  sheet  which  we  have  discussed  here  is  called 
after  its  discoverer  a  symmetric  Helmholtz  wave.  It  is  convenient  to  have 
short  names  for  the  two  non-tilting  states  of  the  wave.  The  state  with 
deformed  vortex  siieet  and  the  upper  and  lower  wave  (vertical  motion)  in 
opposite  phase  is  called  the  a-state.  Ihe  state  with  a  non-deformed 
vortex  sheet  and  the  upper  and  lov/er  wave  ^  phase  is  called  the  b-state. 
From  either  one  of  these  the  symmetric  Helmholtz  wave  will  approach  asymp¬ 
totically  the  same  state  of  stationary  phase  with  a  downwind  tilt  from 

the  a-state  and  exponential  grov/th  of  the  amplitudes  as  long  as  A  kCl. 

s 

The  phase  velocities  with  which  the  wave  leaves  the  non-tilting 
states  may  be  anticipated  from  dimensional  considerations.  The  vortex 
sheet  in  the  undisturbed  state  is  defined  by  a  single  physical  parsuneter 
namely  the  velocity  U  which  measures  one  half  of  the  sliding  vorticity. 

Vi/hen  the  wave  is  introduced  the  wave  length  is  a  second  parameter.  The 
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phase  velocities  in  the  non-tilting  states  must  be  related  to  these  basic 
paurameters,  so  they  must  both  be  proportional  to  U.  The  simplest  choice 
suggests  that  the  wave  moves  in  opposite  direction  through  the  non-tilting 
states  with  the  speed  U  and  approaches  a  state  of  stationary  phase  half 
way  between  these  states.  The  exponential  growth  rate  in  this  asymptotic 
state  has  the  dimension  of  inverse  time,  so  it  must  be  proportional  to  kU. 

It  is  in  fact  equal  to  kU. 

The  values  of  the  stationary  phase  angle  and  the  growth  rate  of  the 
wave  amplitude  are  obtained  mathematically  by  applying  the  kinematic  and 
dynamic  boundary  conditions  to  the  vortex  sheet :  Let  us  first  consider  the 
kinematic  conditions,  which  state  that  the  fluid  particles  in  the  vortex 
sheet  has  the  same  vertical  motion  as  the  fluid  particles  immediately  above 
and  below,  thus 

Dz 

(1.1)  "W  ~  "'J  “  *  (Kinematic  conditions) 

Here  z  is  the  ordinate  of  the  vortex  sheet  measured  from  its  undeformed 
s 

level  (z=0),  and  the  subscripts  1  and  2  refer  to  the  upper  and  the  lower 
layer.  Since  the  motion  in  the  layers  is  solenoidal,  the  velocity  field 
^  of  the  wave  disturbance  may  be  represented  by  a  stream function  f  whose 
sign  we  choose  such  that 

(1.2)  V  =  ui  +  wk  =  X  i. 

Here  k  are  the  orthogonal  triple  of  unit  vectors  of  a  rectangular 

coordinate  system  (x,y,z)  fixed  to  the  symmetric  frame  with  k  pointing 
upward  from  the  lower  to  the  upper  layer,  and  ^  pointing  along  the  basic 
flow  in  the  upper  layer.  ^  points  along  the  axis  of  the  vortex  filaments 
in  the  vortex  sheet.  The  ordinate  z  is  measured  from  the  interface  level. 

Since  the  motion  in  the  layers  is  irrotational,  the  streamfimction  in 
(1.2)  satisfies  Laplaces 's  equation.  Since  further  the  strearafunction  has 
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a  ainuBoldal  variation  in  the  x-direction  with  the  wave  number  k,  its 
Laplace  equation  may  be  written 

7®*  =  =  0. 

The  bounded  solutions  of  this  equation  in  the  two  layers  are 


♦  = 


(in  the  upper  layer) 


kz 

t  =  ^^e  .  (in  the  lower  layer) 

Here  and  denote  the  streamf unction  values  at  the  interface  level 

(z=0).  These  have  the  same  amplitude  in  the  symmetric  wave,  and  by 

inspection  of  fig.  1  we  see  that  the  streamfunction  in  the  upper  layer 

has  the  same  sign  as  z  in  the  a-state.  In  an  arbitrary  tilting  state 

s 

of  the  wave  with  the  phase  8  measured  positive  downwind  from  the  a-state 
the  wave  elements  are  therefore 

z  =  A  cos  kx  =  Real  part  (A  e^^) 

(1.4)  =  UAco6(kx-0)  =  Reed  psurt  (Ze^^) 


-  ♦,  =  UAcos(kx  +  9)  ■  Real  part  (Z*e^^)  . 

On  the  right  are  introduced  the  complex  wave  parameter  Z  euid  its  conjugate 
Z*,  defined  by 

(1.5)  Z  =  UAe"^®  »  X  -iY. 

By  including  U  in  the  amplitude  factor  UA  of  the  streamfunction  the  amplitude 
A  has  the  dimension  of  length.  It  is  readily  seen  that  A  is  the  interface 
level  streamline  amplitude  of  the  resultant  field  in  the  symmetric  frame. 

< -  Let  Zj  denote  the  ordinate  of  this  streamline  in  the  upper  layer. 

The  streamline  slope  at  an  arbitrary  tine,  to  the  lineetr  approximation,  is 


which,  using  (1.4),  gives 


z,  =  A  coe(kx  -  0) 
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lliis  is  a  sine  curve  with  the  amplitude  A  and  the  crest  displaced  the  phase 

9  dovmwind  from  the  interface  crest  ais  shown  in.  the  middle  diagram,  fig.  1. 

Since  the  interface  defcrmaticn  is  staticnary,  the  streamline  is  parallel 

to  the  interface  at  the  node,  so  its  amplitude  is  given  by  A  =  A  cos  9. 

6 

This  condition  is  of  course  implicitly  contained  in  the  formal  kinematic 
condition  (1.1). 

We  note  that  the  symmetric  v;ave  in  (1.4),  besides  the  two  unchanging 
parameters  U  and  k,  is  defined  by  the  two  amplitudes  A^,  A  and  the  phase  9- 
These  three  parameters  are  functions  of  time  only.  To  find  their  evolution 
in  time,  we  need  three  independent  conditions.  Tv;o  of  these  are  the 
kinematic  conditions  in  (1.1).  The  third  is  the  dynamic  condition  that 
the  vorticity  of  the  filaments  in  the  vortex  sheet  (the  sliding  vorticity) 
does  not  change  with  time. 

It  is  convenient  to  represent  these  conditons  non-dimensionally  by 

choosing  units  of  length  and  time  such  that  U=k=l,  so  these  parameters 

disappear  from  the  equations.  By  this  choice  the  growth  rate  of  the  wave 

amplitudes  will  be  obtained  non-dimensionally  in  units  of  kU. 

Turning  now  to  the  kinematic  condition  (l.l)  above  the  interface, 

we  represent  the  individual  change  as  the  sum  of  local  and  convective 

changes  in  the  symmetric  frame  and  express  the  vertical  velocity  by  the 

streamfunction  in  (1.2).  With  the  non-linear  convective  term  v”  7^ 

~  ~  s 

ignored,  the  condition  then  takes  the  form 


If  this  condition  is  satisfied  for  the  symmetric  wave,  the  kinematic  condition 
below  the  interface  is  automatically  satisfied.  One  of  the  kinematic 
conditions  in  (1.1)  is  replaced  by  the  condition  of  symmetry  in  (1.4),  so 
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it  Is  sufficient  to  consider  (1.6).  With  the  wave  elements  substituted  here 
from  (1.4),  the  exponential  factor  e^^  cemcels,  and  we  get 

A  +iA  =  iZ, 

which,  using  (1.5)  have  the  real  and  imaginary  parts 

• 

Ag  =  Y  =  A  sin  0 

(1.7)  (Kinematic  conditions). 

A  s  X  =  A  cos  6. 
s 

The  first  of  these  gives  two  equivalent  expressions  for  the  vertical  velocity 
at  the  interface  crest.  The  second  condition  states  the  fact,  already 
noted  earlier,  that  the  streamline  is  parallel  to  the  interface  at  the  nodes. 
The  ratio  of  the  two  kinematic  conditions  gives  the  growth  rate  of  the 
interface  amplitude  in  an  arbitrary  tilting  state  of  the  wave,  namely 

(1.8) 

The  growth  rate  of  the  wave  amplitude  is  obtained  from  the  time  derivative 

of  the  second  condition  combined  with  the  first,  namely 

•  •  • 

A  =  A  cos  9  -  9  A  sin  9  =  A  sin  9. 

6 

Applied  to  the  b-state  (9=90°)  this  combined  kinematic  condition  gives 
the  phase  velocity  9^=-l.  The  wave  moves  upwind  through  the  b-state  with 
the  speed  U  of  the  basic  flow.  In  the  b-state  the  propagation  of  the  wave 
is  governed  by  kinematic  conditions  only.  The  combined  condition  gives  the 
growth  rate  of  the  wave  amplitude  in  eui  arbitrary  tilting  state,  namely 
(1.9)  T^ln  A  =  (l  +  9)tan9. 

The  wave  is  neutral  in  the  non-tilting  states.  In  a  state  of  stationary 
phase  (9=0)  the  wave  and  the  interface  deformation  grow  at  the  same  rate, 
n»tan0  .  In  all  other  tilting  states  they  grow  at  different  rates. 

0 

We  now  tvirn  to  the  dynamic  boundary  condition  at  the  interface,  which 
states  that  the  surface  vorticity  of  the  vortex  sheet  filaments  is  con- 
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servative,  that  is 

^  (2U  +  Uj  -  Ug  )  s  0. 

With  the  individual  change  represented  as  the  sum  of  the  local  and  con¬ 
vective  changes  in  the  symmetric  frame  this  condition,  to  the  linear  approx¬ 
imation,  has  the  form 

(1.10)  -We)  =  “  (ui  +  Ug ).  (dynamic  condition). 

It  is  the  mathematical  formulation  of  the  principle  which  we  made  use  of 
in  the  qualitative  discussion  of  the  a-state  at  the  beginning  of  this 
section,  namely:  The  local  increase  of  vorticity  along  the  vortex  sheet 
comes  from  the  convergence  of  the  surface  vorticity.  This  interpretation 
of  the  dynamic  boundary  condition  was  given  by  Hoiland*  in  1942.  The 
principle  is  most  clearly  exliibited  in  the  a-state  (top  diagram  in  fig.  l) 
where  the  convergence  along  the  vortex  sheet  has  maximum  strength  at  t’le 
dovmwind  nodes. 

With  the  horizontal  velocities  in  (l.lO)  represented  by  the  stream- 
function  in  (1.2,3)  the  dynamic  condition  becomes 

•!,)  -  -t.). 

or  with  the  wave  elements  of  the  symmetric  wave  in  (1.4)  substituted 

(1.11)  Bind)  M  A  cos  9.  (Uk=l). 

This  is  the  dynamic  condition  (conservation  of  sliding  vorticity)  for  the 
symmetric  wave  in  the  vortex  sheet.  With  the  time  differentiation  carried 
out,  it  becomes 

•  • 

A  sin  9  -t-  9  A  cos  6  =  A  cos  6. 

Applied  to  the  a-state  it  gives  the  phase  velocity  of  the  wave  through  the 
a-state  as  9=1.  The  wave  moves  downwind  through  the  a-state  with  the 

*  E.  Hoiland:  The  Developed  Form  of  the  dynamic  Boundary  Condition  with 
Applications,  Arch,  for  Math,  og  Naturvid.  B.  XLVI.  Nr.  2,  1942. 
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speed  of  the  basic  flow,  so  it  has  at  that  moment  no  intrinsic  propagation 
through  the  fluid.  The  propagation  through  the  a-state  is  determined  by 
the  dynamic  condition  only,  namely  by  the  vorticity  concentrations  at  the 
downwind  nodes,  presicely  as  anticipated  by  qualitative  reasoning  earlier. 

The  dynamic  condition  gives  another  expression  for  the  growth  rate 
of  the  wave  amplitude,  namely 

•^InA  =  (l-9)cot0,  (Dynamic  condition) 

(1.12)  ^ 

•^InA  =  (l  +  9)tan9.  (Kinematic  condition). 

The  ratio  of  the  two  conditions  gives  the  phase  speed  in  an  arbitrary 
tilting  state,  namely 

(1.13)  9  =  cos®9-sin®9  =  cos26. 

This  formula  verifies  the  earlier  values  in  the  non-tilting  states,  and 
it  represents  the  phase  velocity  in  an  surbitrary  tilting  state  as  a  linear 
combination  of  the  values  in  the  non-tilting  states.  In  particular  it  shows 

that  the  wave  from  either  non-tilting  state  approaches  the  state  of  stationary 

« 

phase  9 -*9  *45°,  S-O. 

8 

With  the  value  of  0  substituted,  either  equation  (1.12)  gives  the 
explicit  value  of  the  growth  rate  of  the  wave  amplitude  in  8ui  arbitrary 
tilting  state,  namely 

(1.14)  -^InA  =  sin  29  =  -  >4 In  cos  20. 

By  integration  this  equation  gives  the  evolution  of  the  wave  from  an  arbitrary 
initial  state  with  the  phase  9^  and  the  amplitude  A^,  namely 

(1.15)  A®cos20  ■  A*  cos 20  *  A*9  . 

o  o 

A®  is  proportional  to  the  kinetic  wave  energy,  so  0A*  is  the  kinetic 
wave  energy  transport,  which  accordingly  is  constant  during  the  evolution 
of  the  wave.  As  thf  wave  slows  down,  the  kinetic  energy  of  the  wave  increases 
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a  proportional  amount.  If  the  wave  initially  has  a  phase  upwind  from  the 

a-state  less  than  45°  the  wave  moves  downwind  toward  the  a-state.  While 

approaching  the  a-state  it  speeds  up  and  decays.  It  moves  through  the 

a-state  with  maximum  speed  (the  speed  of  the  fluid)  and  minimum  amplitude. 

From  the  a-state  the  wave  moves  on  downwind  toward  the  state  of  stationary 

phase  9 -•  9  =  45°  and  grows  as  it  slows  down.  This  evolution,  which  is 
s 

expressed  mathematically  by  the  formula  in  (1.15) »  is  represented  graphically 

by  a  rectangular  hyperbola  in  a  polar  diagram  with  A, 9  as  the  polar  coordinate 

Each  point  on  this  hyperbola  represents  a  state  of  the  wave.  During  the 

evolution  the  "state-point"  moves  along  the  hyperbola  toward  increasing  0 

such  that  the  radius  vector  moves  with  a  constant  areal  velocity  0A®. 

If  the  wave  has  a  phase  less  than  45°  downwind  from  the  b-state  it  moves 

upwind  from  this  state,  and  the  evolution  is  represented  graphically  by  the 

conjugate  hyperbola  which, during  the  evolution,  is  traced  out  in  the  sense 

of  decreasing  9  through  the  b-state  and  on  toward  asymptotic  approach  to 

the  growing  state  of  stationary  phase. 

We  note  from  (I.I3)  that  the  v/ave  has  tv/o  states  of  stationary 

phase,  namely  the  downwind  state  of  0  =  45°  where  the  amplitude  grows  at  the 

s 

exponential  rate  1 «  kU,  and  the  symmetrically  located  upwind  state  (0  =  -  45°  ■ 

s 

where  the  wave  decays  at  the  same  rate.  These  two  states  of  stationary  phase 
are  known  as  the  normal  modes  of  the  vortex  sheet.  An  arbitrary  tilting 
state  of  the  wave  may  be  represented  as  the  resultant  of  the  two  normal 
modes  with  the  proper  amplitude  ratio.  As  one  mode  grows  and  the  other 
decays  the  resultant  tilting  wave  moves  on  toward  the  asymptotic  growing 
state.  An  arbitrary  tilting  state  of  the  wave  may  also  be  represented  as 
the  resultant  of  an  a-state  and  a  b-state  with  the  proper  amplitude  ratio. 

As  these  two  component  waves  move  on  toward  their  common  state  of  stationary 
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phase  in  accordance  with  the  evolution  in  (I.15),  their  resultant  has  the  evolu¬ 
tion  of  the  resul taint  wave  which  again  obeys  (1.1^).  The  pair  of  normsuL 
modes  (the  growing  and  decaying  states  of  stationary  phase)  are  thus  closely 
related  to  the  pair  of  non-tilting  states  of  the  symmetric  wave.  We  shall 
consider  this  relationship  further  in  the  next  section. 

The  evolution  in  (I.I5)  was  derived  from  lineeui'ized  theory  which  gives 
a  valid  approximation  only  as  long  as  the  amplitude  is  a  smadl  fraction 
of  the  wave  length  (kA^^l).  Let  us  consider  very  briefly  what  happens 
after  the  linear  approximation  brakes  down.  We  have  seen  (1.12)  that  the 
linear  theory  predicts  the  growth  rate  of  the  wave  amplitude  in  the  state 
of  stationary  phase  as  l  =  kU=2it/T,  where  T=I/U  is  the  period  during  vdiich 
the  fluid  moves  a  wave  length  in  the  symmetric  frame.  If  in  the  state  of 
stationary  phase  the  wave  at  some  initial  time  (t=0)  has  the  amplitude  A^, 
its  amplitude  at  the  time  t  has  grown  to  the  value  A =  A^exp(2nt/T).  The 
wave  amplitude  doubles  during  the  time  interval  t  =  0.11T,  while  the  fluid 
moves  a  tenth  of  a  wave  length,  so  the  motion  appears  to  be  highly  unstable. 
However,  this  exponential  growth  will  only  continue  as  long  as  the  amplitude 
is  very  much  smaller  than  the  wavelength,  when  the  linearized  theory  gives 
a  valid  approximation,  loiter  on,  as  the  amplitude  grows,  the  non-linear 
terms  become  important  and  can  no  longer  be  ignored.  The  non-linear  evolution 
must  be  studied  with  the  complete  non-linear  dynamic  equations,  mdiioh  must 
be  integrated  numerically  in  short  time  steps.  Such  integrations  were 
first  carried  out  by  Rosenhead*  in  1951.  The  integration  has  been  repeated 
more  recently  by  several  people  with  high  speed  electronic  computers.  These 
experiments  have  confirmed  Rosenhead' s  earlier  results.  Instead  of  the 
continuous  distribution  of  vortex  filaments  along  the  vortex  sheet,  Rosenhead 

*  The  Formation  of  Vorticies  from  a  Surface  of  Discontinuity.  Proo,  Roy. 

Soc.  1952,  Series  A,  Vol.  154,  p.  170-192.  -  — 
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carried  out  the  computations  with  twelve  vortex  filaments  unifoiraly 
distributed  per  wave  length.  His  results  are  shown  in  fig.  la.  The  concen¬ 
tration  of  the  vortex  filaments  toward  the  dovmwind  nodes,  which  was  predicted 
by  the  linear  theory,  is  clearly  reflected  by  Rosenhead's  discontinuous 
model  during  the  evolution  from  the  init’al  symmetric  state  at  the  top  to 
the  next  state  below  a  quarter  period  later.  However,  the  amplitude  grows 
at  a  much  slower  rate  than  the  linear  theory  predicts.  As  the  non-linear 
evolution  continues,  the  growth  of  the  amplitude  soon  stops,  while  the 
vorticity  becomes  increasingly  concentrated  on  the  steepening  downwind  slope c. 
The  vorticity  concentration  is  accompanied  by  a  rolling  up  of  the  vortex 
sheet  around  the  downwind  nodes,  a  feature  which  was  noted  already  by 
Helmholtz.  The  rolling  up  proceeds  with  increasing  rapidity  for  the 
filaments  near  the  nodal  points. 

There  can  be  little  doubt  that  a  vortex  sheet  which  is  subjected 
to  a  small  amplitude  initial  periodic  disturbance,  in  the  absence  of 
viscous  diffusion  of  the  vorticity,  behaves  essentially  in  the  way  Roienhead’ c 
computations  predict:  The  wave  crests  at  first  grow  almost  symmetrically 
at  the  exponential  rate  predicted  by  the  linear  theory  until  their  heights 
are  about  one  tenth  of  the  wave  length.  The  subsequent  non-linear  evolution 
is  asymmetric  with  gradual  concentration  cf  the  sliding  vorticity  towara 
the  nodal  points  on  the  steepening  dovmwind  slopes.  As  this  process  con¬ 
tinues,  the  vortex  sheet  rolls  up  around  the  downwind  nodes,  and  in  the 
end  the  vortex  sheet  appears  broken  up  into  periodica,  y  spaced  vortices 
with  most  of  the  vorticity  confined  within  circuleu:  regions  whose  diameter 
is  about  four  tenths  of  the  wave  length.  The  time  required  for  this  process 
from  an  initial  deformation  whose  amplitude  is  one  per  cent  of  the  wave 
length  to  the  final  vortex  is  roughly  the  time  it  takes  the  fluid  in  either 
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layer  to  move  a  wave  length. 

13ie  vortex  sheet  is  therefore  an  unstable  hydrodynamic al  system. 

The  smallest  disturbance,  whether  periodic  or  not,  will  result  In  an 
evolution  of  the  type  described  above  during  which  the  energy  of  the  distur¬ 
bance  grows  at  the  expense  of  the  energy  of  the  basic  flow.  This  type  of 
hydrodynamic  Instability  is  known  as  shearing  instability.  The  Instability 
will  be  partly  surpressed  if  the  fluid  in  the  layers  has  different  density 
with  the  lighter  fluid  on  top,  as  will  be  shown  In  the  next  section.  The 
vortex  sheet  between  fluid  layers  of  the  same  density  exhibits  the  shearing 
instability  in  its  pure  form. 

l^ere  is,  however,  a  stabilizing  effect  in  the  vortex  sheet  which  we 
have  not  considered,  and  idiich  cannot  be  removed,  namely  the  molecular  viscous 
diffusion  of  vorticity.  This  diffusion  process  is  going  on  all  the  time  in 
real  fluids  in  regions  where  there  are  vorticity  gradients.  This  process 
is  of  course  particularly  active  near  the  center  of  the  rolled  up  vortex 
sheet,  and  tends  to  convert  these  regions  into  space  vortices  with  continuous 
vorticity  distribution.  However,  the  diffusion  tends  to  modify  the  vortex 
sheet  itself  even  in  the  absence  of  any  disturbance.  This  modification  must 
be  taken  into  account  if  we  wish  to  predict  the  behavior  of  disturbances 
of  a  vortex  sheet  in  a  real  fluid.  This  problem  will  be  discussed  in  sections 
6  and  7  below. 

2.  Symmetric  waves  in  a  statically  stable  vortex  sheet.  -  Let  now 
the  unbounded  homogeneous  fluid  layers  with  the  relative  translation  ZU 
have  different  densities,  with  the  lighter  fluid  (pa  Pi  )  above  and  the 
heavier  fluid  (ps  Pg)  below  the  interface  which  separates  the  layers.  The 
static  stability  of  the  interface  is  measured  by  the  non-dimensional  parameter 
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(2.1)  a  =  (p»  -  Pi  )/(Pi  +  Pi  )  <1.  (Quasi-symmetry  condition;. 

This  system  is  not  symmetric,  for  the  lower  layer  has  greater  inertia. 

However,  if  the  density  difference  is  a  small  fraction  of  the  average  density 
as  indicated  in  (2.1),  the  asymmetry  of  the  inertia  can  be  ignored  to  a  good 
approximation.  The  justification  of  this  quasi-symmetry  approximation  is 
discussed  in  section  4. 

If  the  symmetry  approximation  is  justified,  a  symmetric  wave  with 
the  wave  elements  in  (1.4)  will  preserve  its  symmetry  at  all  times.  The 
evolutions  of  amplitudes  and  phase  aire  governed  by  the  two  kinematic  conditions 
in  (1.7)  and  a  dynamic  condition  which  describes  the  changes  of  the  sliding 
vorticity  of  the  vortex  filaments  in  the  interface.  These  changes  are  most 
clearly  visualized  in  the  a-state  of  the  wave  (top  diagram  in  fig.  l).  The 
convergence  along  the  interface  concentrates  the  sliding  vorticity  toward  the 
downwind  nodes,  as  in  the  case  of  the  statically  neutral  vortex  sheet. 

However,  the  fluid  filaments  at  these  nodes  are  now  light  at  the  top  and 
heavy  at  the  bottom,  and  they  are  tilting  downstream  with  the  slope  of  the 
interface,  so  the  upstream  side  is  heavier  than  the  downstream  side.  The 
gravitational  buoyancy  of  this  "over  weight"  will  tend  to  return  the  nodal 
filaments  to  their  equilibrium  level  positions,  and  thereby  induce  them  to 
rotate  in  the  sense  opposite  to  the  sliding  vorticity.  The  resultant 
vorticity  change  is  less  than  the  partial  change  from  either  of  these  causes. 

If  the  convergence  dominates,  the  wave  moves  downwind  through  the  a-state 
with  a  slower  speed  than  the  fluid.  The  wave  moves  upwind  through  the  b-state 
with  unit  speed.  Therefore  it  has  a  state  of  stationary  phase  less  than 
45°  downwind  from  the  a-state,  and  its  growth  rate  there,  n=  tan0  ,  is  less 
than  unity.  If  the  overweight  dominates,  the  wave  moves  upwind  through  the 
a-state  and  upwind  through  the  b-state.  It  has  no  states  of  stationary  phase. 
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It  moves  upwind  at  all  times  with  rhythmic  variations  of  phase  speed  and 
amplitude  between  the  extreme  values  in  the  non-tilting  states,  while  the 
energy  transport  in  the  symmetric  frame,  8 A’,  remains  constant. 

The  vorticity  change  from  the  overweight  must  be  proportional  to  the 
density  difference  and  the  fi'^ce  of  gravity,  that  is  gCp^  -Pi).  It  must 
therefore  be  proportional  to  the  non-dimensional  parameter 
(2,2)  u  =  sg/kU®  (Richardson  number), 

where  s  is  the  non-dimensional  static  stability  parameter  in  (2.1).  The 
parameter  a  is  called  the  Richardson  number  of  the  v/ave. 

To  find  the  proportionality  factor,  we  consider  the  a-state  of  the 
wave  for  vdiich  the  vorticity  change  from  convergence  is  just  balanced  by  the 
change  from  the  overweight,  so  the  resultant  vorticity  change  at  the  nodes 
is  zero.  This  wave  is  stationary  in  the  a-state.  It  has  a  stationary  neutred 
a-state  with  no  local  accelerations,  so  the  convective  accelerations  are 
balanced  by  the  acting  forces  (gravity  and  pressure  force).  But  the  con¬ 
vective  accelerations  are  not  changed  if  the  direction  of  the  flow  is  changed 

in  one  of  the  layers,  say  the  lower  layer  (see  top  diagram  in  fig.  1).  This 

new  field  is  also  steady.  But  this  is  the  field  of  a  gravity  wave  in  a 
statically  stable  interface  with  no  shear,  in  the  frame  wdiich  moves  with 
the  speed  U  of  the  wave.  The  speed  of  propagation  of  the  gravity  wave  in 
the  interface  between  unbounded  layers  is  given  by  U®  =  sg/ki  or  m,=  I« 

The  wave  mechanism  of  the  gravity  wave  is  very  simple:  Looking  at  the  a-state 

diagram  in  fig.  1  with  the  direction  of  the  flow  reversed  in  the  lower  layer 
and  inspecting  the  adjacent  streamline  channels  above  and  below  the  inter¬ 
face,  we  note  that  the  upper  flow  has  maximum  speed  over  the  crest  and  minimum 
speed  in  the  trough,  while  the  speed  below  the  Interface  changes  in  opposite 
rhythm.  This  means  that  the  Interface  filament  on  the  crest  has  maximum 


positive  vorticity  and  the  filament  in  the  trough  has  maximum  negative  vorticity. 
The  interface  is  one  of  the  streamlines,  and  the  filaments  move  downwind 
along  it  with  the  speed  U  of  the  wave.  On  the  journey  from  crest  to  trough 
the  positive  vorticity  on  the  crest  is  changed  to  the  negative  vorticity 
in  the  trough.  This  change  is  caused  by  the  overweight  in  the  downstream 
hill  between  crest  and  trough.  For  a  given  overweight  (density  difference) 
the  change  proceeds  at  just  the  right  rate  when  the  filament  moves  along 
the  streamline  with  the  speed  U®  =  sg/k. 

The  local  vorticity  change  along  a  statically  stable  vortex  sheet 
is  accordingly  given  by  the  formula 

■^(ui  -Ug)  =  -  (l-u)U'^(ui  +U9X  (Dynamic  condition). 

The  first  term  on  the  right  represents  the  vorticity  concentration  by  con¬ 
vergence  (see  1.10),  the  u-term  gives  the  vorticity  change  from  the  overweight. 
The  two  effects  act  in  opposite  sense,  and  the  ratio  of  their  magnitudes 
is  the  Richardson  number  of  the  wave  which  is  proportional  to  the  wave  length. 

The  (u=  l)-wave  has  no  local  vorticity  change.  Ihis  wave  has  a  stationary 
neutral  a-state  which  is  dynamically  equivalent  to  the  gravity  wave  we  obtain 
by  reversing  the  direction  of  the  flow  in  one  of  the  layers. 

The  dynamic  condition  has  been  anticipated  here  from  the  speed  of  the 
gravity  wave.  Because  of  its  importance  we  shall  now  derive  it  from  basic 
dynamic  principles.  Following  Hoiland's  procedure  (l.c.)  it  is  obtained 
by  integrating  the  equation  of  motion  along  a  closed  rectangular  curve 
with  horizontal  and  vertical  sides.  Its  diagonal  is  an  infinitesimal  segment 
of  the  interface.  The  pressure  is  continuous  across  the  interface  (the 
dynamic  condition)  so  the  pressure  integral  is  zero,  and  we  get 

j*  =  0. 
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The  rectangular  curve  has  only  four  infinitesimal  segments,  one  horizontal 
and  one  vertical  in  each  layer.  With  second  order  (non-linear)  terms  ig¬ 
nored  the  Integral  becomes 

(2.3)  ]^(piUi  -o,u»)dx  =  g(pa  -Oi)dZg. 

We  now  introduce  the  quasl-i^mmetry  approximation  by  assuming  that  the 
difference  of  inertia  between  the  layers  can  be  Ignored,  and  replace  the 
densities  Pi  and  in  the  Inertial  terms  on  the  left  by  the  average  density, 
)i(Pa-i-Pi).  ^e  dynamic  condition  for  the  quasi-symmetric  system  then  takes 
the  simple  form 

(2.4)  IJt  -u*)  =  2sg-j^,  (Dynamic  condition), 

where  a  is  the  stability  parameter  in  (2.1).  This  condition  is  intuitively 
clear:  The  sliding  vorticity  of  the  interface  filament  changes  at  a  rate 
which  is  proportional  to  the  stability  and  the  slope  of  the  interface.  The 
lower  half  of  the  filament  is  heavy  and  the  upper  half  is  light,  so  the  over¬ 
weight  tends  to  retvirn  the  filament  to  its  level  position.  With  the  individual 
changes  in  (2.4)  represented  as  the  sum  of  local  and  convective  changes,  the 
condition  (to  the  linear  approximation)  takes  the  form 

d  d 

"St  (ui  •  (Dynamic  condition). 

Here  the  local  change  of  the  sliding  vorticity  along  the  interface  appears 
as  the  resultant  of  two  effects,  namely :  (i)  the  concentration  of  the  sliding 

vorticity  of  the  basic  flow  by  convergence  of  the  vortex  filamPnts  along 
the  Interface,  and  (ii)  the  action  of  the  overweight  of  the  filaments  which 
tends  to  return  them  to  their  level  positions.  To  compare  the  relative 
magnitude  of  the  two  effects  we  evaluate  the  interface  slope  from  the  difference 
of  the  kinematic  conditions  in  (1.1),  which  gives 
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Wlth  this  value  of  the  slope  substituted,  the  dynamic  condition  takes  the 
form  idilch  was  anticipated  from  the  dynamics  of  the  gravity  waves,  namely 

(2.5)  =  -  (1 -»i)U'^(ui  +Uj).  (dynamic  condition) 

Having  derived  this  condition  from  basic  principles  t  it  now  Includes  the  theory 
of  the  gravity  wave  as  equivalent  to  the  a-state  of  the  (usl)-wave,  and 
predicts  the  speed  of  this  wave  as  s  sgA* 

With  the  wave  elements  of  the  symmetric  wave  In  (1.4)  substituted, 
the  dynamic  condition  for  this  wave  becomes 
d  *  * 

(2.6)  (A  sin  9)  =  A  sin  0  +  9  A  cos  9  =  (1  -  u)A  cos  0  . 

Applied  to  the  a-state  of  the  wave  (030)  It  shows  that  the  wave  In  that 
state  has  the  speed  6  =  1  -  u.  The  short  (u<l)-waves  move  downstream,  the 

ft 

longer  (u>l)-waves  move  upstream  through  the  a-state.  All  waves  move  up¬ 
stream  through  the  b-state  with  the  speed  U  of  the  flow,  governed  by  the 
kinematic  conditions  which  are  Independent  of  the  static  stability  across 
the  Interface.  So  the  short  waves  are  unstable  with  a  growing  state  of 
stationary  phase,  while  the  long  waves  are  stable,  moving  upstream  with 
rhythmic  variations  of  phase  speed  and  amplitude. 

The  details  of  the  evolution  are  obtained  by  comparing  the  two  equivalent 
expressions  for  the  growth  rate  of  the  wave  amplitude  from  the  kinematic 
and  dynamic  conditions, 

•j^lnA  =  (l-u-0)cot9,  (Dynamic  condition) 

(2.7) 

■^InA  =  (l+0)tan0.  (Kinematic  condition). 

The  ratio  of  these  give  the  phase  speed  in  an  arbitrary  state  as  a  linear 

combination  of  the  non-tilting  speeds 

•  •  • 

(2.8)  0  =  (1  -  u)cos®0  -  sin’s  =  8  cos’ 9  +  0.  sin’ 0  , 

a  b 
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The  short  ((i<l}-waveB  have  two  states  of  stationary  phase,  namely 

(2.9)  tan*e„  =  -6/0.  =  1  -  u-  (0=0). 

The  dynamic  and  kinematic  conditions  give  two  equivalent  values  of  the  growth 
rate  in  these  states,  namely 

=  (i-u)coteg  =  e^cote^ 

(2.10) 

=  “  =  tan  0g  =  -  9^  tan  9^ 

and  their  product  determines  the  growth  rate 

e  « 

(2.11)  n^  =  -9  9.  =  1 -u  =  -m®. 

a  D 

The  growing  state  (n>0)  has  the  positive  phase  9  downstream  from  the  a-state. 

8 

l^e  decaying  state  (n  <  0)  has  the  same  negative  phase  located  symmetrically 
upstream  from  the  a-state. 

The  long  (u>l)-waves  move  upstream  with  rhythmic  variations  of  speed 
between  the  extreme  values  in  the  non-tilting  states. 

The  evolution  of  the  wave  from  an  arbitrary  tilting  state  is  obtained 
by  substituting  the  value  of  0  in  (2.8)  into  either  one  of  the  basic  conditions 
(2.7),  which  gives 


d  d  * 

=  (2  -  n)sin  0  cos  0  =  -  In  0  . 

The  integral  of  this  from  an  arbitrary  initial  state  with  the  amplitude 

A  and  the  phase  9  is 
o  o 

(2.12)  A*0  =  A*C(l  -  u)cos®0  -  6in*6]  =  A®[(1  -  ji)co8“0  -sin*0  ]. 

O  0  0 

It  shows  that  the  energy  transport  of  the  wave  is  constant.  If  the  initial 

state  is  an  a-state  with  the  amplitude  A  s  A  ,  the  equation  becomes 

o  & 

y  a  ,  .a 


(2.13) 


(1  -  U  _  ^A  ^n  9^ 


1-  u. 


For  the  short  (u,<l)-wavos  this  is  a  hyperbola  in  a  (A,  0) -polar  diagram 


Fig  2  Evolution  of  short  unstoble  Helmholtz  wove 
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with  the  asymptotic  slope  tan  9  =  /l  -  u,  =  n.  In  the  evolution  from  the 

s 

a-state  (0=0)  the  area  between  the  hyperbola  and  the  asymptote  is  swept 
out  with  the  constant  areal  velocity 

•  u 

HA*e  =  0  =  (A^  tan  0  )n, 

a  a  a  s 

so  the  areal  velocity  is  proportional  to  the  asymptotic  growth  rate,  as 

indicated  in  fig.  2.  The  evolution  from  an  initial  b-state  (0=90°)  is 

represented  by  the  corresponding  conjugate  hyperbola.  The  evolution  from 

an  arbitrary  initial  state  (A  ,9  )  is  obtained  by  choosing  the  corresponding 

o  o 

point  on  the  hyperbola  as  the  initial  point. 

For  the  long  (u,>l)-wave6  equation  (2.13)  is  an  ellipse  with  the 

amplitude  ratio  A^/A^  =  /u -1.  =  m,  which  is  traced  out  in  the  upstream  sense 

of  decreasing  9  during  the  evolution  of  the  wave.  The  area  of  the  ellipse 

is  swept  out  at  the  constant  rate 

>40A2  =  A»  =  )4e.A»  =  ramA  A.  , 

a  a  b  b  a  b 

so  the  entire  area  of  the  ellipse,  uA^A^,  is  swept  over  during  the  time 
T=2n/m.  This  is  the  time  during  which  the  wave  moves  two  full  wavelengths 
while  the  interface  performs  a  complete  standing  oscillation.  So  T  is  the 
period  and  m  in  (2.11)  is  the  frequency  of  these  stable  oscillations.  We 
note  that  the  asymptotic  growth  rate  of  the  vmstable  waves  and  the  frequency 
of  the  stable  waves  are  both  equal  to  the  geometric  mean  of  the  absolute  values 
of  the  phase  velocities  in  the  non-tilting  states,  measured  in  the  xmits 
of  U  and  k. 

All  waves  move  upwind  through  the  b-state  with  unit  speed.  The  stable 
(>i<2)-wave8  move  slower  through  the  a-state  so  they  decay  during  the  passage 
from  the  a-state  to  the  b-state.  The  longer  (ji>2)-waves  move  with  a  speed 
greater  than  unity  through  the  a-state,  so  these  waves  grow  from  the  a-state 
to  the  b-state.  The  (tA=2)-wave  is  a  neutral  wave  which  moves  upwind  with 
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conatant  unit  speed  and  constant  amplitude. 

3.  Relation  between  the  non-tilting  states  and  the  normal  modes  of 
the  symmetric  wave.  -  As  mentioned  in  section  1,  the  states  of  stationary 
phase  of  the  xinstable  symmetric  waves  are  called  the  normal  modes  because 
the  wave  grows  at  a  constant  exponential  rate.  The  statically  stable 
quaei-symmetric  vortex  sheet  between  unbounded  layers,  defined  by  the  parameters 
(U,8),  has  a  pair  of  unstable  normal  modes  of  this  type  for  every  wave  length 
shorter  than  the  stationary  a-wave.  For  longer  wave  lengths  the  normal  modes 
are  a  pair  of  neutral  waves  with  the  interface  deformation  propagating  to 
the  right  or  to  the  left  with  constant  amplitude  and  the  constant  non- 
dimensional  phase  speed  m  in  (2.11)  relative  to  the  symmetric  frame.  In 
the  frame  which  moves  with  the  mode  the  interface  is  a  stationary  streamline. 

The  speeds  of  the  flow  in  the  layers  relative  to  this  stationary  frame  have 
the  ratio  (l-m)/(l  +  m),  and  the  amplitudes  of  the  upper  and  lower  part  of 
the  wave  have  evidently  the  same  ratio.  So  the  stable  modes  are  not  symmetric. 
However,  the  resultant  of  such  a  pair  with  the  same  amplitude  is  the  stable 
symmetric  wave  with  the  standing  oscillation  of  the  interface  and  the  upstream 
propagation  of  the  wave  as  described  at  the  end  of  section  2.  The  properties 
of  the  long  stable  (u>l)-modeB  are  therefore  known  implicitly  from  the  stable 
oscillations. 

The  boundary  conditions  (2.7)  for  the  symmetric  wave  gave  the  relations 
(2.10)  between  the  phase  velocities  in  the  non-tilting  states  and  the  growth 

rate  and  phase  of  the  unstable  modes,  namely 

•  • 

(3.1)  0  =  ntan9  ,  0.  =  -ncot0  ,  (n>0) 

a  B  b  s 

where  0  denotes  the  phase  of  the  growing  mode  measured  positive  downwind  from 

O 

the  a-state.  These  formulas  are  characteristic  of  all  types  of  symmetric 
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wavee,  regardless  of  the  physical  mechanism  which  governs  their  propagation 
and  growth.  To  show  thls«  consider  an  arbltreLry  pair  of  unstable  modes  with 
the  growth  rates  ±nand  the  phase  ^  each  of  these  modes  taken 

separately  the  upper  field  has  the  evolution 

^  coe(x  -  0g), 

(3.2  _  (n>0) 

ti  "  =  /ae"  cos(x+  0  ), 

8 

where  an  arbitrary  amplitude  factor  is  left  out  in  each  field.  The 

resultant  of  these  modes  with  equal  initial  amplitudes  has  the  upper  field 

(3*3)  =  chntcos0  cosx  +  sh  nt  sin  0  sinx  =  Acos(x-0). 

s  s 

This  is  a  symmetric  wave  starting  from  an  initial  a-state  (0=t=O)  with 

the  amplitude  A  =  cos  0  ,  and  the  phase  of  the  wave  has  the  evolution 
&  8 

(3.^)  tan  0  =  tan  0  th  nt. 

s 

The  wave  moves  toward  the  asymptotic  state  0-»0  as  t-»*,  with  the  phase 

8 

speed 

e 

(3.5)  9  =  n  tan  9  (cos9/chnt)®  =  n  tan  9  (A  /A)®, 

s  s  a 

e 

in  other  words  with  the  constant  energy  transport  0  A*  =  n  tan  0^  A® .  With 
ch*®nt =  1  -  th®nt  eliminated,  using  (5.^)»  the  phase  speed  formula  becomes 

•  e  • 

(3*6)  0  =  ntan0  cos®  0  -  n  cot  0  sin®0  =  0  cos®0  +  0.  8in*0. 

s  s  a  ^ 

The  result  to  the  right  is  obtedned  by  applying  the  formula  to  the  non-tilting 
states  (0=0,90°).  But  these  values  of  the  phsise  speeds  in  the  non-tilting 
states  are  the  values  in  (3»1)  which  were  obtained  by  examining  the  symmetric 
wave  as  an  entity  in  its  evolution  from  an  arbitrary  state.  In  other  words 
the  symmetric  wave  heis  the  same  evolution  whether  it  is  regarded  as  an 
entity  or  as  a  resultant  of  normal  modes.  Hiis  is  a  consequence  of  the  linear 
theory.  If  several  component  waves  satisfy  the  linearized  boundary  conditions, 
their  resultant  will  also  satisfy  these  conditions  at  all  times.  We  may 
therefore  use  the  spectrum  of  normal  modes  with  their  simple  exponential 
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growth  rates  to  exaalne  the  early  linear  stage  in  the  evolution  of  arbitrary 
non-periodic  local  Initial  disturbances  by  representing  these  local  distur¬ 
bances  as  the  resultant  of  the  entire  spectrum  of  symmetric  waves  with  their 
amplitudes  and  phase  properly  adjusted  to  fit  the  local  disturbance  at 
the  initial  moment,  using  the  general  method  of  Fourier.  Ihe  evolution 
of  the  local  disturbance,  for  as  long  as  the  linear  approximation  is  valid, 
is  then  obtained  as  the  resultant  of  the  separate  evolutions  of  the  component 
symmetric  waves,  which  is  known  from  section  2. 

If  the  growth  rate  and  stationeury  phase  of  the  growing  mode  are  known, 
the  phase  velocities  in  the  non-tilting  states  are  given  by  (3»l)-  Con¬ 
versely  if  the  phase  velocities  in  the  non-tilting  states  are  known  the 
growth  rate  and  phase  of  the  modes  are  given  by  the  corresponding  formulas 

(3*7)  n*  =  -0^9^  =  -m*;  ”  "  ®a^®b  ’ 

These  formulas  define  both  the  growing  and  the  decaying  mode.  They  have 
the  same  positive  and  negative  growth  rate  emd  the  same  stationary  phase 
downstream  and  upstream  from  the  a-state.  inie  sign  of  the  phase  of  the 

t 

growing  mode  (n>0)  is  given  by  the  sign  of  6  . 

El 

All  the  characteristics  of  the  symmetric  wave  may  then  be  obtained 

either  from  the  growth  rate  and  phase  of  the  modes  (n,0  )  or  from  the  phase 

8 

t  • 

velocities  in  the  non-tilting  states  of  the  wave  (0J^*0^,)•  Either  pair 
of  parameters  may  be  derived  independently  of  the  other  pair,  which  then 
in  turn  are  obtained  from  the  interrelated  formulas  in  (3.1)  or  (3.7). 

The  question  of  which  pair  should  be  selected  for  primary  derivation  depends 
upon  the  circiunstances,  and  sometimes  it  is  a  matter  of  taste.  The  more 
mathematically  inclined  student  would  probably  in  moat  cases  prefer  to 

determine  the  normal  mode  pair  (n,0  )  first.  The  phase  velocities  in  the 

•  • 

non-tilting  states  (0_,0.  )  on  the  other  hand  are  often  accessible  by  a 

8  D 
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sitnple  physical  consideration,  and  may  sometimes  be  derived  intuitively 
from  fundamental  physical  principles  with  very  little  mathematical  labor. 

However,  both  sets  of  parameters  are  useful  tools  for  the  physical 
analysis  of  the  wave  mechanism:  If  the  wave  moves  in  opposite  directions 
through  the  non-tilting  states  it  starts  growing  when  it  leaves  these 
states,  and  it  approaches  from  either  side  the  growing  state  of  stationary 
phase  whose  growth  rate  and  phase  are  uniquely  determined  by  the  non-tilting 
phase  speeds.  If  the  wave  moves  in  the  same  direction  through  the  non¬ 
tilting  states,  it  is  a  progi-essive  wave  in  that  direction  with  rhythmic 

variations  of  speed  and  amplitude  in  such  a  vay  that  the  energ;/-  transport 

• 

0A»  remains  constant.  Meanwhile  the  interface  performs  standing  oscillations, 
with  the  frequency  equal  to  the  geometric  mean  of  the  non-tilting  phase 
speeds.  Both  growth  rate  and  stationary  phase  are  imaginary  in  the  stable 
wave.  But  the  formula  for  the  field  in  (5.5)  is  real  and  describes  the 
evolution  of  the  stable  wave  from  an  initial  a-state. 

4.  Justification  of  the  symmetry  approximation.  -  We  shall  now 
examine  the  magnitude  of  the  error  we  make  in  the  quasi-symmetric  wave 
theory  by  ignoring  the  kinematic  effect  of  the  density  difference  in  the 
dynamic  boxmdary  condition.  To  determine  the  error  we  must  find  the  behavior 
of  the  waves  in  a  system  with  arbitrary  densities  in  the  layers,  and  examine 
how  much  it  departs  from  the  behavior  of  the  quasi-symmetric  waves  when 
the  density  difference  is  small. 

To  make  the  analysis  a  little  more  general  we  shall  let  the  layers 
be  bounded  by  rigid  horizontal  planes  at  equal  distances  above  and  below 
the  interface,  so  both  layers  have  the  same  finite  depth  h.  This  system 
has  a  new  length  parameter  h,  besides  the  wave  number k,  but  the  theory  for  the 
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waves  in  the  bounded  system  must  reduce  to  the  wave  theory  for  unbounded 
layers  for  the  waves  which  are  very  much  shorter  than  the  depth  of  the 
layers  (kh»l).  In  fact,  we  see  from  (1.3)  that  the  amplitude  of  the  wave 
in  the  imbounded  system  decreases  exponentially  with  the  distance  from  the 
interface  and  is  reduced  to  about  4  per  cent  of  the  central  value  at  a 
distance  of  half  a  wave  length  from  the  interface.  If  the  rigid  boundaries 
are  further  away  they  have  very  little  effect  on  the  waves.  Hov/ever,  if 
they  are  closer,  the  streamfunction  of  the  wave  must  be  such  that  the  vertical 
motion  dt/dx  is  zero  at  the  boundaries.  Since  the  streamfunction  is 
periodic  in  the  horizontal  direction,  this  means  that  it  is  zero  at  the 
outer  rigid  boundaries.  The  motion  is  of  course  still  Laplacean,  so  the 
streamfunction  must  be  a  linear  combination  of  the  basic  solutions  in  (1.3) 
which  is  zero  at  the  height  h  above  and  at  the  depth  h  below  the  inter¬ 
face.  In  other  words,  the  streamfunction  of  the  wave  in  the  bounded 
system  must  have  the  form 

t  =  ti  sh  k(h  -  z)/shkh,  (Upper  layer) 

(4.1) 

If  =  ^  sh  k(h  +  z)/8hkh,  (Lower  layer) 

where^  as  before,  and  ^  are  the  values  in  the  two  layers  at  the  interface 
level.  The  corresponding  horizontal  velocities  au-e 

u  =  k  chk(h  -  z)/Bh  kh,  (Upper  leiyer) 

(4.2) 

-u  =  k  chk(h  +  z)/sh  kh.  (Lower  layer) 

We  see  that  the  wave  has  no  vertical  wind-shear  near  the  rigid  boundaries. 
This  may  be  anticipated  from  the  fact  that  the  irrotii  ^onal  flow  is  here 
parallel  to  the  boundaries,  so  both  the  curvature  and  shear  are  zero.  We 
shall  consider  the  consequences  of  this  circumstance  later. 

If  the  lower  fluid  is  much  heavier  than  the  upper  fluid,  the  inertia 
of  the  system  is  highly  asymmetric  with  respect  to  the  interface.  It  is 
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therefore  no  reason  to  expect  that  an  initially  symmetric  wave  with  the 
wave  elements  in  (1.4)  will  remain  symmetric  during  its  evolution.  Rather 
the  contrary  must  be  expected.  We  shall  therefore  investigate  whether  this 
asymmetric  system  has  normal  modes  for  which  the  interface  deformation 
has  the  evolution 

(4.3)  z  =  A  X  Real  part  of  exp  i (lex  +  rat), 

s  so 

The  frequency  equation,  m=m(k)  for  the  mode  is  found  by  applying  the 
two  kinematic  conditions  (1.1 )  and  the  dynamic  condition  (2.3)  to  the  interfac‘‘ 
deformation.  If  for  a  given  value  of  k  the  frequency  equation  gives  a  real 
positive  value  for  m,  the  mode  is  a  stable  neutral  wave  which  moves  in  the 
symmetric  frame  toward  negative  x  with  the  phase  speed  m/k..  If  the  frequency 
equation  gives  an  imaginary  frequency  m=  -  in,  the  mode  is  an  unstable  wave 
whose  amplitude  grows  at  the  exponential  rate  n. 

We  shsdl  again  represent  the  boundary  conditions  non-dimensionally 
by  using  units  of  length  and  time  such  that  k=U=l.  Ilie  local  change  in 
the  symmetric  frame  of  the  deformation  in  (4.3)  is  then 


dz 

ir=  ^“*s 


m  ^^8 

=  inr 


» 


so  the  kinematic  conditions  (1.1)  for  the  mode  eu'e 
(4.4)  z  =  ti/(ii>  +  l)  =  ^^/(ra-l). 

s 

The  dynamic  condition  in  (2.3)  is 

]5^(PiUi  -  OaUg  )  =  g(pa  -  Pi  )  . 

From  (4.2)  the  horizontal  accelerations  at  the  interface  level  are 


^=cthkh^;  .^=cthkh^;  (k=l) 


or  with  the  streamfunction  values  of  the  mode  substituted  from  (4.4) 
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Tv,  Dz  dz 

=  (m  +  1 )  cth  kh  =  (m  +  1  )*  cth  kh  • 

Du. 

“  =  (m  - 1)  cth  kh  =  (m  -  l)®  cth  kh  • 

With  these  substituted  the  dynamic  condition  becomes 
Pi  (m  +  1)*  +  p*  (m  -  1)*  =  g(p(  -  0i  )thkh. 

This  is  the  frequency  equation  for  the  normal  modes  (4.3)  in  the  statically 
stable  vortex  sheet  between  bounded  layers.  The  Richardson  number  of  the 
waves  in  the  bounded  system  is  defined  as 

(4.5)  '1  =  (sgAl^)th:<;h. 

It  reduces  to  the  value  in  the  unbounded  system  when  the  wave  is  very  much 
shorter  than  the  depth  of  the  layers.  V^ith  this  value  of  u  substituted  the 
frequency  equation  becomes 

m*  -  2sm  =  u  -  1 

or,  with  s*  added  on  both  sides, 

(4.6)  (ra-s)»  =  u- (l-5>). 

If  the  density  difference  is  a  small  fraction  of  the  mean  density  (s^l)  the 
frequencies  of  the  normal  modes  are  obtained  with  good  approximation  by 
ignoring  s  in  (4.6),  giving  m*  =  u-1  =  -n*,  which  is  the  value  vie  obtained 
in  (2.11)  by  introducing  the  quasi -symmetry  approximation  in  the  dynaunic 
boundary  condition.  We  can  now  test  the  approximation  quantitatively  by 

I 

comparing  the  correct  value  of  the  frequency  in  (4.6)  with  the  approximate 
value  for  any  given  value  of  s.  Atmospheric  inversions  seldom  exceed  a 
temperature  jump  of  10®C.  Ihe  corresponding  s-value  is  less  than  0.02,  so 
practically  no  error  is  made  by  putting  1  -  s*  =  1  in  (4.6).  The  quasi- 
symmetric  vadue  of  the  frequency  differs  from  the  correct  value  by  an  amount 
roughly  equeil  to  s.  The  quasi-symmetric  theory  is  good  to  this  approximation. 

When  the  density  difference  between  the  layers  is  appreciable  (as  for 
example  the  interface  between  the  ocean  amd  the  atmosphere),  the  dynamic 
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asymmetry  across  the  Interface  can  of  course  not  be  ignored.  In  the  limiting 
case  with  zero  density  in  the  upper  layer  (s  =  1)  the  system  becomes  one  single 
layer  with  a  free  surface,  and  the  frequencies  of  the  normal  modes  in  this 
layer  are  given  by 

(m-l)*  =  jA  =  (gAU*)thkh. 

Ihe  normal  modes  (4.3)  are  here  two  families  of  neutral  gravity  waves  which 
move  in  "the  symmetric  frame"  toward  negative  x  with  the  non-dimensional 
phase  velocities  m  =  l±/u.  The  fluid  in  the  layer  moves  in  the  symmetric 
frame  toward  negative  x  with  the  speed  U  =  1.  Relative  to  the  fluid  the  modes 
move  of  course  symmetrically  in  opposite  directions  with  the  same  speed, 

U®  =  (gA)thkh. 

Vdien  the  densities  in  the  layers  have  arbitrairy  values  with  the  heavier 
fluid  in  the  lower  layer,  the  static  stability  s  has  an  intermediate  value 
between  the  limiting  values  zero  and  one,  and  the  normal  modes  of  the  system 
have  the  frequencies  in  (4.6).  The  long  modes,  u>l-s*,  have  real  frequen¬ 
cies,  whose  values  may  be  represented  graphically  in  a  non-dimensional  m,u- 
diagram  (see  fig.  4)  by  a  parabola  whose  vertex  has  the  coordinates 

m  =  s, 

(vertex  of  dispersion  parabola) 

The  parabola  has  the  saune  shape  for  all  values  of  s.  As  s  is  changed  con¬ 
tinuously  from  zero  to  one,  the  vertex  moves  along  a  similar  inverted  par¬ 
abola,  m*  =  1-u  ,  which  connects  the  vertex  points  of  the  limiting  vertex 
parabolas  (s  =  0,1).  With  the  coordinates  of  the  vertex  substituted,  the 
frequency  equation  (4.6)  becomes 

a 

(4.7)  (m-m)  =  n-u. 

6  8 

The  long  (u>|ji  )-modes  in  (4.3)  are  neutral  waves  which  propagate  in  the 
symmetric  frame  with  the  non-dimensional  phase  velocities  m  in  (4.7).  Let 
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U8  now  use  another  frame  of  reference  which  moves  toward  negative  x  relative 
to  the  symmetric  frame  with  the  speed  U  =  m  U  s  sU.  In  this  frame  the  non- 
dimensional  phase  velocity  of  the  mode  is  m  =  m  -  m  •  euid  the  dispersion 
equation  for  the  mode  has  the  simplest  possible  foi'm,  namely 
(4.8)  DOo  =  u-  n_* 

B 

The  frame  moves  relative  to  the  symmetric  frame  with  the  speed  U  such  that 

s 

(Pa  +  Pi  =  Pa^-  PjU.  It  is  the  frame  in  which  the  center  of  mass  of  the 
s 

system  is  stationary.  The  long  stable  neutral  modes  ((*>(i_)  move  in  opposite 

s 

directions  relative  to  the  center  of  mass  of  the  system  with  the  same  phase 

velocity  (n-u  )^.  The  short  unstable  (ti<u  )-modes  have  stationeiry  phase 
8  8 

relative  to  the  center  of  mass  and  they  grow  or  decay  at  the  exponential 

rate  (u  -  The  treuisitional  (u  =  ia  )-mode  is  a  neutral  wave  which  is 

8  8 

stationary  relative  to  the  center  of  mass. 

All  the  dispersion  (or  frequency)  parabolas  in  fig.  4  have  one  point 
in  common«with  the  coordinates  m  =  0,  u  =  !•  This  means  that,  regardless  of 
the  density  difference  between  the  layers,  one  of  the  stable  modes  of  the 
system  is  stationeiry  in  the  symmetric  frame,  namely  the  mode  for  vdiich 
U®  =  (sg/k)thkh.  If  for  this  stationary  mode  the  direction  of  the  flow  is 
reversed  in  one  of  the  layers,  the  acceleration  remains  unchanged  everywhere. 
So  the  new  situation  is  edso  stationary.  But  this  is  the  field  of  a  gravity 
wave  in  the  frame  which  moves  with  the  speed  U  of  the  wave  relative  to  the 
fluid  in  the  layers.  In  other  words,  the  gravity  wave  propagates  through  the 
fluid  with  the  speed  U  which  makes  its  Richardson  number  ^  in  (4.3)  equal  to 
unity.  This  speed  gives  the  overweight  the  right  time  to  reverse  the  vorti- 
city  of  the  interface  filaments  diu-lng  their  passage  from  crest  to  trou^. 

If  the  wave  is  very  much  shorter  than  the  depth  of  the  layers  (kh>  1), 
the  hyperbolic  factors  in  (4.1  and  4.2)  are  very  neeurly  equal  to  the  simple 
exponential  functions,  exp(7kz),  and  the  Richardson  number  of  the  wave  in 
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(4.3)  has  very  nearly  the  value  sg/k^  of  the  unbounded  system.  Ihe  boundaries 
have  no  appreciable  effect  on  the  waves  which  are  much  shorter  than  the  depth 
of  the  layers.  13ie  short  wave  approximation  with  the  boundary  effects  Ignored 
Is  In  fact  quite  good  If  the  depth  of  the  layer  Is  more  than  half  a  wave  length. 
For  waves  which  are  very  much  longer  than  the  depth  of  the  layers  (kh«  l)  we 
obtain  another  long-wave  approximation  by  Ignoring  squares  and  higher  powers 
of  kh  In  the  general  formulas.  It  turns  out  that  the  same  long-wave  approxi¬ 
mation  la  obtained  when  we  Ignore  the  vertical  acceleration  In  the  vertical 
dynamic  equation,  and  take  the  pressure  distribution  along  the  verticauL  to  be 
hydrostatic. 

5.  The  quasi-static  approximation.  -  Let  us  first  consider  the  approxi¬ 
mation  as  a  long  wave  approximation  of  the  general  theory  in  seJtion  4.  If 
kh<  1  we  may  develop  the  hyperbolic  functions  (4.1  and  4.2)  in  series  and 
ignore  higher  powers  of  kh.  The  streamf unction  in  the  upper  layer  to  the 
order  of  (kh)*  is 

(5.1)  ♦  =  ti(l-f)Cl  +  ^f(f-2)(kh)*]  (z>0). 

To  the  same  order  of  approximation  the  upper  horizontal  velocity  field  is 

(5.2)  hu  =  ♦![!+ ()6  p- -  I  +  >i)(kh)*]  (z>0). 

The  Richardson  number  of  the  wave,  to  the  same  approximation,  is 

(5.3)  =  (sgh/U®)[l  -  yi(kh)*]. 

If  the  wave  is  much  longer  than  h,  the  factors  in  the  brackets  in  these  formulas 
differ  from  unity  by  an  amount  less  than  (kh)*.  The  long  wave  approximation  is 
obtained  by  taking  these  factors  equal  to  one.  Tbe  error  in  this  approximation 
is  only  a  fraction  of  (kh)*. 

We  note  that  the  long-wave  approximation  gives  a  motion  with  no  vertical 
shear.  Since  the  motion  is  irrotational  it  must  have  enough  shear  to  compensate 


-32- 


for  the  curvature  of  the  streamlines.  But  the  streamlines  have  very  slight 
curvature  in  the  very  long  waves,  and  the  shear  is  therefore  so  sli^t  in  the 
long  waves  that  it  may  be  ignored  without  significant  error. 

!Die  long  wave  approximation  gives  the  vertical  velocity  at  an  arbitrary 
level  in  the  upper  layer  as 

(5.4)  ^ 

We  note  that  this  approximation  satisfies  the  solenoidal  condition  2,‘Z,®  ® 
the  incompressible  fluid.  On  the  other  haind,  the  lamellar  condition  7 X  v  =  0 
is  slightly  violated,  since  the  v/ind  sheai‘  of  the  slightly  curved  flov/  is 
ignored.  The  vertical  velocity  in  (5.4)  increases  linearly  from  zero  at  the 
upper  rigid  boundary  to  its  maximum  value  at  the  interface  level,  namely 

(5.5)  .5  .his.. 

e«  Dt  Dt  • 

In  the  absence  of  shear  the  vertical  fluid  column  remains  vertical,  and  its 
height  hi  decreases  at  the  rate  of  the  vertical  motion  at  its  base,  which  in 
turn  is  equal  to  the  horizontal  divergence  in  the  column. 

Let  us  now  compare  this  long  wave  approximation  with  the  quasi-static 
approximation  which  is  made  when  the  dynamic  effect  of  the  vertical  accelera¬ 
tion  on  the  pressure  is  ignored,  so  the  pressure  distribution  ailong  the  vertical 
is  taken  to  be  hydrostatic.  Ihe  horizontal  and  vertical  dynamic  eqiiations  for 
quasi-static  flow  are  therefore 

PlV^t  =  -  7p, 

(5*6)  (Quasi-static  equations) 

pg  =  -  8p/dz. 

Let  us  apply  these  equations  to  an  arbitrary  motion  of  a  homogenous  fluid 
layer.  Since  the  pressure  distribution  along  the  vertical  is  hydrostatic,  the 
vertical  thickness  of  the  isobaric  layers  is  the  same  everywhere  in  a  homo¬ 
genous  layer.  All  the  isobaric  surfaces  have  the  same  shape  and  may  be  gene¬ 
rated  by  displacing  one  of  them  along  the  vertlcad.  Ihis  means  that  the 
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horizontal  pressure  gradient  has  no  variation  with  height,  and  therefore  the 
horizontal  acceleration  has  no  variation  with  height  in  a  quasi-static  motion 
of  a  homogenous  fluid  layer.  If  the  motion  has  no  vertical  shear  at  a 
certain  time,  no  shear  can  develop  later.  Ihe  vertical  fluid  columns  will 
remain  vertical  at  all  times.  We  have  seen  that  this  condition  is  very  nearly 
satisfied  in  waves  which  are  much  longer  than  the  depth  of  the  layer,  so  the 
motion  in  these  waves  obeys  the  quasi-static  approximation  with  good  accuracy. 
The  kinematic  conditions  in  (5*5)  which  were  derived  from  the  long  wave 
app*’oximation  may  also  be  derived  from  the  quasi-static  theory,  because  they 
are  the  consequence  of  the  condition  that  the  motion  has  no  vertical  shear. 

In  other  words,  the  long  wave  approximation  is  equivalent  to  the  quasi-static 
approximation,  and  vice  versa.  It  may  be  noted  from  (5 •5)  that  the  ratio 
between  the  amplitudes  of  the  vertical  and  the  horizontal  components  of  the 
motion  is  equal  to  kh,  that  is  the  ratio  between  the  vertical  and  horizontal 
scale  of  the  motion.  Quite  generally  then,  the  quasi-static  equations  (5.6) 
may  be  used  as  a  good  approximation  for  the  description  of  the  motion  if  the 
ratio  between  the  horizontal  and  vertical  scale  of  the  motion  is  smadl.  The 
percentual  error  we  make  by  using  these  equations  has  the  order  of  the 
square  of  this  ratio.  We  shall  show  later  in  Chapter  II  that  the  same  con¬ 
dition  for  the  quasi-static  approximation  applies  to  compressible  flow  with 
a  continuous  variation  of  the  density  along  the  vertical. 

The  Bichardson  number  has  the  constant  value  sgh/U”  for  all  waves  which 
are  long  enough  to  obey  the  quasi-static  approximation  with  good  accuracy. 

This  means  that  the  long  quasi-static  gravity  waves  in  the  resting  double 
layer  all  move  with  the  same  phase  speed  thi'ough  the  fluid,  where  =  sgh. 
The  quasi-static  waves  are  non-dispersive.  This  speed  exceeds  the  correct 
value  by  about  0.7(kh)’  per  cent.  The  phase  speed  of  the  long  quasi-static 
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gravlty  waves  may  be  taken  as  a  convenient  measure  for  the  static  stability 
of  the  resting  double  layer.  When  the  layers  have  the  relative  translation 
2U,  the  Richardson  number  for  the  corresponding  quasi-static  Helmholtz  waves  is 
(5*6)  u  =  sgh/U®  =  (Cj^/U)®.  (Quasi-static) 

It  is  a  non-dimensional  constant  which  includes  all  the  physical  parameters 
of  the  system,  and  it  measures  the  ratio  between  the  static  stability  and  the 
shear  across  the  interface.  In  the  quasi-symmetric  systems  (s<  1),  the 
Richardson  number  defines  the  dynamic  properties  of  the  quasi-static  waves: 

If  the  static  stability  dominates  (u>l)  the  quasi-static  waves  are  stable 
and  have  the  frequency  ra  =  / u -  1 .  If  the  shear  dominates  (ii<l)  the  quasi¬ 
static  waves  are  unstable  with  the  growth  rate  n  =  /I  -  (i.  In  the  systems 
with  arbitrary  densities  in  the  layers  the  transition  from  stable  to  unstable 
quasi-static  waves  occurs  at  a  Richardson  number  less  than  unity,  namely 
u^=  l-s». 

6.  Lateral  diffusion  of  vorticity  from  a  shear  layer.  -  In  order  to 
see  how  the  vorticity  diffusion  effects  the  dynamic  behavior  of  a  vortex 
sheet,  we  shall  examine  how  the  undisturbed  vortex  sheet  is  modified  by  the 
diffusion  process. 

Let  us  begin  by  considering  a  more  general  system  with  the  vorticity 
concentrated  in  a  central  layer  between  unbounded  irrotational  layers.  Uie 
boundaries  between  the  central  layer  and  the  outer  layers  are  not  sharp. 

Hie  vorticity  is  supposed  to  decrease  symmetrically  in  both  directions  from 
a  meiximum  value  at  the  center  of  the  shear  layer  and  approach  zero  asymp¬ 
totically  with  increasing  distance  from  the  center  level  (see  fig.  6). 

The  diffusion  process  is  particularly  simple  if  the  vorticity  distri¬ 
bution  at  some  initial  moment  is  given  by  the  Oaussean  error  function 


Fig.  6  Symmetric  vorlicity  distribution  ocrost  o  central  layer 
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(6.1)  <1  =  <lo  ®xp[  -  n(z/d)*  ]  =  dU/dz. 

In  the  synunetrlc  frame  which  moves  with  the  fluid  at  the  central  level  the 
asymptotic  maximum  speed  In  the  outer  layers  Is 

Uq  =  /  q  dz  =  q©  J 
o  o 

With  this  value  of  substituted,  the  Oaussean  vortlclty  distribution  In 

(6.1)  becomes 

(6.2)  q=  (2Uo/d)  expC- n(z/d)*]. 

The  parameter  d  may  be  regarded  as  the  depth  of  the  central  layer  In  an 
"eqtilvalent"  three  layer  system  with  the  constant  central  level  vortlclty  in 
a  shear  layer  of  depth  d  between  imbounded  Irrotational  outer  layers  which 
have  the  relative  translation  ZUq, 

If  the  vortlclty  distribution  in  the  layer  is  Qaussean  at  any  one  time, 
it  will  remain  Ckiussean  at  all  times,  while  the  depth  grows  with  time.  This 
is  readily  seen  by  substitution  of  q  from  (6.2)  Into  the  Navler-Stokes  vortl- 
city  equation  for  two 'dimensional  shesir  flow,  namely 

(6.3)  q  =  v7®q  =  pq".  (p  =  kinematic  viscosity) 
Here  the  accent  denotes  hel^t  differentiation  and  the  dot  denotes  local  time 
differentiation  in  the  symmetric  frame.  The  logarithmic  fora  of  (6.2)  is 

In  q  =  ln(2Uo)-ln  d-n(z/d)*. 

The  logarithmic  height  derivative  is 

q'/q  =  -2nz/d*. 

The  second  hel^t  derivative  is 


The  logarithmic  time  derivative  is 
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With  these  substituted  the  vorticity  equation  (6.3)  gives 
(6.4)  d  =  2TTv/d, 

which  has  the  integral 

d*  -  dj  =  4iiv(t  -  to). 

The  rate  of  growth  of  the  depth  of  the  Gaussean  shear  layer  is  proportioned 
to  the  inverse  depth  of  the  layer.  We  may  choose  the  initial  sheeur  layer  as 
thin  as  we  like.  In  the  limit,  if  v/e  choose  do-*0,  we  have  initially  a  vortex 
sheet,  and  we  see  that  the  vortex  sheet  immediately  is  changed  by  the  viscous 
vorticity  diffusion  into  a  shear  layer  of  finite  depth,  and  this  layer  grows 
in  thickness  at  the  rate  in  (6.4).  It  is  clear  then  that  the  dynamic  theory 
for  the  instability  of  the  vortex  sheet  in  section  1  (which  is  based  on  the 
assumption  that  the  vortex  sheet  retains  its  identity  at  edl  times)  is  not 
very  satisfactory.  To  improve  on  the  theory  we  must  examine  the  instability 
of  a  shear  layer  of  finite  depth. 

7.  The  instability  of  a  constant  shear  layer.  -  The  Gaussean  shear 
layer  which  develops  by  the  vorticity  diffusion  from  a  vortex  sheet  is  not 
easy  to  discuss  mathematically.  We  shall  therefore  in  this  section  examine 
a  shear  layer  of  depth  d  with  the  constant  wind  shear, 

q  =  2U/d, 

between  unbounded  irrotational  outer  leiyers  which  have  the  relative  trans¬ 
lation  2U.  All  three  layers  have  the  same  constant  density.  We  shall 
ignore  the  vorticity  diffusion  at  the  boundaries  of  the  sheeu:  layer. 

Let  the  upper  interface  at  some  initiad  time  be  given  a  small  ampli¬ 
tude  sinusoidal  deformation, 

(7.1)  Zi  =  A  COB  kx,  (kA  <1) 

s  s 

while  for  the  moment  the  lower  interface  is  undeformed.  Since  the  vorticity 
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is  conserved  in  this  system  (see  60)f  this  deformation  meeuis  that  the  con¬ 
stant  vorticity  in  the  shear  layer  has  been  added  in  the  sinusoidal  segments 
above  the  interface  level,  and  the  same  vorticity  is  removed  from  the  segments 
below  the  interface  level. 

The  field  associated  with  the  sinusoidal  deformation  of  the  upper  inter¬ 
face  may  be  represented  as  the  resultant  of  two  component  fields.  One  compcaieni 
is  the  undisturbed  shear  layer  with  the  undeformed  interface.  The  second 
component,  the  ^-field,  is  irrotational  everyv;here  except  in  the  areas  of  adc'ed 
vorticities  between  the  interface  level  and  the  deformed  interface.  It  is  easy 
to  see  what  the  general  appearance  of  this  periodic  field  must  be.  It  has  a 
slight  asymmetry  with  reference  to  the  interface  level  since  all  the  positive 
vorticity  is  above  and  all  the  negative  vorticity  is  belov/  the  interface  level. 
However,  if  the  interface  aunplitude  is  a  small  fraction  of  the  vra.ve  length, 
this  asymmetry  may  be  ignored  (the  linear  approximation).  The  field  is  then 
very  nearly  symmetric  with  reference  to  the  interface  level  with  sinusoidal 
distribution  of  vortex  filaments  along  the  interface  level  and  zero  vorticity 
elsewhere. 

All  the  vortex  filaments  may  be  grouped  in  pairs  distributed  symmetri¬ 
cally  with  reference  to  any  one  of  the  nodes.  The  vertical  line  through  the 
node  is  the  central  straight  streamline  of  each  pair,  so  it  is  a  streamline 
of  the  resultant  periodic  ;v-field. 

It  is  readily  seen  that  this  periodic  field  will  propagate  the  defor¬ 
mation  upwind  relative  to  the  fluid  at  the  interface.  The  wave  speed  is 
obtained  by  application  of  Stokes's  theorem  to  a  small  rectangular  section 
of  the  sinusoidal  segment  of  added  vorticity.  The  area  of  this  section  is 
Zidx.  The  contribution  to  the  circulation  round  the  segment  along  the  verti¬ 
cal  sides,  ZjW,  are  sraauLl  of  the  second  order.  The  circixlation  round  the 
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•TMi  sxdx  to  the  first  order  (linear)  approximation  is 

(7>2)  (u'*'-u~)dx  >  qsidx  B  (2U/d}zidx.  (The  dynamic  condition) 

Here  denotes  the  tangential  component  of  the  Laplacean  v-field  above  the 
sinusoidal  segment  of  added  vorticity  and  u'  denotes  the  tangential  component 
below.  This  is  the  dynamic  condition  for  the  deformed  interface. 

Let  the  v-field  be  represented  by  a  streamfunctlon  (as  in  1.2),  whose 
value  at  the  interface  level  (z  ■  0)  be  denoted  by  •  Above  and  below  the 
interface  the  streamfunction  (as  in  1.3)  has  the  values 

The  tangential  wind  shear  across  the  segment  of  added  vorticity  is  then 

u*  -  u"  =  2k'^  . 

With  this  value  substituted  the  dynamic  condition  (7.2)  takes  the  simple  form 
(7.3)  n  =  (U/h)zi  (K  =  kd) 

This  field  propagates  the  deformation  upwind  relative  to  the  fluid  with  a 
certain  phase  velocity  Cg.  In  a  frame  of  reference  which  moves  with  this 
wave  the  wave  is  stationary  and  the  deformed  interface  is  one  of  the  stream¬ 
lines,  having  the  slope 

or  kinematic  condition) 

This  is  the  kinematic  condition  at  the  interface.  The  kinematic  and  dynamic 
conditions  combined  give  the  intrinsic  upwind  phase  velocity  of  the  wave  at 
the  moment  when  only  the  upper  interface  is  deformed,  namely 
(7.V)  Cg  .  U/a,  Cg/U  -  V2xd. 

This  wave  is  called  a  Rayleigh  wave  aifter  Lord  Rayleigh  who  discovered  it  in 

1880. • 

nie  v-field  associated  with  the  added  vorticity  at  the  upper  deformed 
interface  gives  a  vertical  motion  to  the  fluid  particles  at  the  lower  initially 

4  '  '  ■ ' 

Lord  Rayleigh:  On  the  Stability  and  Instability  of  Certain  Fluid  Motions. 
Scientific  Papers.  Article  66,  Vol.  I,  p.  474. 
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level  interface.  In  the  next  Inetant  the  lower  Interface  will  therefore  have 
a  sinusoidal  deformation  with  crests  and  troughs  in  phase  with  the  upper 
Interface  nodes.  Ibis  lower  deformation  in  turn  is  associated  with  a  periodic 
component  field  of  the  type  (7.3)  centered  at  the  lower  Interface.  The  evolu¬ 
tion  of  the  wave  is  therefore  governed  by  the  resultant  of  the  two  periodic 
fields  from  the  added  vorticitles  at  both  interfaces*  and  this  evolution  is 
not  simple. 

However,  if  at  the  initial  moment  both  interfaces  are  given  equal 
sinusoidal  deformations  with  the  same  amplitude  and  phase  we  have  a  symmetric 
wave.  Because  of  the  kinematic  and  dynamic  symmetry  of  the  shear  layer  this 
wave  will  remain  symmetric  at  all  times.  Similar  to  the  Helmholtz  wave  in 
the  vortex  sheet  (see  fig.  l)  we  call  the  non-tilting  state  of  the  symmetric 
Bayleigh  wave  with  the  interfaces  deformed  in  phase  and  the  upper  and  lower 
component  fields  in  opposite  phase  the  a-state  of  the  wave.  The  second 
non-tilting  state  vdth  the  interfaces  deformed  in  opposite  phase  and  the 
component  field  in  phase  is  called  the  b-state.  (See  the  wave  diagrams  in 
fig.  7.) 

In  the  a-state  the  upper  and  lower  component  fields  have  opposite 
phase.  The  straight  vertical  streamlines  of  each  component  (pansing  through 
the  interface  nodes)  coincide  but  they  have  opposite  sense.  The  upper  com¬ 
ponent,  if  operating  alone,  would  propeigate  the  upper  deformation  intrinsically 
upwind  through  the  fluid  with  the  speed  U/h,  The  lower  component  field  tends 
to  oppose  this  propagation.  However,  since  the  intensity  of  the  field  decreases 
exponentially  with  the  distance  from  the  Interface,  it  is  reduced  by  the  factor, 
(7.5)  or  =  e"*'*  (K  =  kd) 

at  the  upper  interface  level.  Its  contribution  to  the  propagation  of  the  upper 
wave  is  reduced  by  the  same  factor.  The  upper  deformation  therefore  propagates 
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Intrlnslcally  upwind  through  the  a-state  with  the  speed 

(7.6)  C,  =  U(l-a)/H. 

But  the  fluid  at  the  Interface  moves  downwind  with  the  speed  U  In  the  symmetric 
frame,  so  the  upper  wave  moves  downwind  through  the  a-state  with  the  speed 
U  -  C,  in  the  symmetric  frame.  Because  of  the  symmetry  the  lower  wave  heis  the 
same  downwind  speed. 

In  the  b-state  with  the  interfaces  deformed  in  opposite  phase  (see  the 
lower  wave  diagram  in  fig.  7)  the  upper  and  lower  component  fields  have  the 
same  phase.  They  propagate  the  deformations  intrinsically  upwind  through  the 
fluid  with  the  speed 

(7.7)  Ck  =  U(l  +  a)/H. 

In  the  symmetric  frame  the  symmetric  wave  moves  downwind  through  the  b-state 
with  the  speed  U-C^.  The  arithmetic  mean  of  the  intrinsic  phase  velocities 
in  the  non-tilting  states  is 

(7.8)  )6(C, +Cb)  =  U/k  =  Co. 

It  is  readily  seen  that  the  wave  has  this  intrinsic  speed  in  a  state  half  way 
between  the  non-tilting  states. 

The  phase  velocities  in  the  non-tilting  states  are  represented  graphi- 
caG.ly  by  the  C, -line  and  the  C^-line  in  the  diagram  to  the  left  in  fig.  7. 

The  Co-line  is  the  straight  broken  line  through  the  origin  and  the  point 
(2nd,  U).  It  represents  the  contribution  to  the  propagation  of  the  deforma¬ 
tion  from  its  own  local  field.  This  contribution  is  proportional  to  the  wave 
length.  The  C, -line  and  the  C^-line  are  located  symmetrically  on  either  side 
of  this  line.  For  waves  shorter  than  the  double  depth  of  the  shear  layer 
(h>3)  the  exponential  term  Qr=  e”^ < 0.04,  so  for  these  short  waves  the  C, - 
and  C* -lines  are  practically  coincident  with  the  central  Cg-iine.  The  reason 
is  that  the  intensities  of  the  component  fields  are  insignificant  at  the  level 
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of  the  other  Interface  In  the  very  short  wave,  so  each  component  wave  behaves 
as  If  the  other  were  absent.  As  the  wave  length  is  Increased  the  action  from 
the  other  component  field  becomes  increasingly  important.  The  Ca-line  departs 
more  and  more  from  the  central  Cg-line,  and  the  Cb-line  departs  the  same 
distance  to  the  other  side.  For  the  (H  =  l)-wave  the  departure  is  U/e.  Ulti¬ 
mately  for  waves  which  are  very  much  longer  than  the  depth  of  the  shear  layer 
(k<1)  we  may  develop  the  exponential  term  in  series  and  ignore  the  higher 
terms 

C,  sld-e"**)  -U(1->4k).  (forL»2iid). 

Hie  C. -line  has  the  long  wave  asymptote  C,  s  U.  Using  this  result  in  (7*8)  we 
find  that  the  Cb-line  has  the  asymptote  Cb -•  U(2/h  -  l).  It  meets  the  C, - 
asymptote  at  hsI  and  has  twice  the  slope  of  the  CQ-line. 

Since  C,  <U  for  all  wave  lengths,  all  waves  move  downwind  through  the 
a-state  in  the  symmetric  frame.  As  soon  as  the  wave  leaves  the  a-state  the 
lower  field  begins  to  augment  the  amplitude  of  the  upper  deformation  and 
vice  versa.  The  wave  starts  growing  when  it  leaves  the  a-state.  In  the 
a-state  the  shear  layer  as  a  whole  has  a  sinusoidad  deformation  with  uniform 
depth  everywhere.  As  the  wave  moves  downwind  from  the  a-state  the  shear  layer 
becomes  periodically  inflated  and  deflated  with  the  inflated  parts  centered 
downwind  from  the  interface  crests.  In  these  respects  the  shear  layer  and 
the  vortex  sheet  have  a  similar  behavior.  In  both  the  symmetric  wave  moves 
downwind  through  the  a-state  while  the  vorticity  bejomes  concentrated  towau*ds 
the  nodes  downwind  from  the  crests. 

The  wave  of  Cb  =  U  has  a  stationary  neutral  b-state.  This  wave  has 
the  wave  number 

-H 

(7.9)  Cb  =  U:  K  =l  +  e  ®  =  1.2785*  **  (L=4.9d). 

S  8 

It  is  the  wave  length  for  which  the  Cb-line  meets  the  U-line  in  fig.  7. 
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Waves  shorter  than  the  stationary  neutral  b-wave  move  downwind  through 
both  non-tilting  states,  so  they  move  progressively  downwind  at  all  times. 

They  pass  through  the  a-state  with  maximum  speed  and  minimum  amplitude.  They 
grow  and  slow  down  after  leaving  the  a-state  and  pass  through  the  b-state 
with  minimum  speed  and  maximum  amplitude.  They  are  stable  waves  with  rhythmic 
variations  of  speed  and  energy  between  the  extreme  values  in  the  non-tilting 
states. 

The  longer  (H<H^)-waves  move  downwind  through  the  a-state  and  upwind 
through  the  b-state.  From  either  non-tilting  state  they  approach  asymptotically 
a  state  of  stationary  phase  with  a  downv/ind  tilt  from  the  a-state,  as  shown  in 
fig.  7  for  the  wave  k=0.8.  As  this  asymptotic  state  of  stationary  phase  is 
approached  the  wave  keeps  on  growing,  ultimately  at  the  rate  which  is  character¬ 
istic  for  the  state  of  starionary  phase.  This  growth  rate  must  be  proportional 
to  the  instantaneous  strength  of  the  component  fields,  and  hence  to  the  instant¬ 
aneous  interface  amplitude  A  ,  Denoting  the  proportionality  factor  by  n,  the 

s 

instantaneous  growth  rate  is  dA g/dt  =  nA^.  In  the  asymptotic  state  of  stationary 
phase  the  long  symmetric  waves  in  the  shear  layer  grow  at  a  constant  exponential 
rate. 

The  evolution  of  the  symmetric  wave  from  an  arbitrary  state  is  obtained 
by  applying  the  boundary  conditions  to  such  a  state.  The  deformations  of  the 
interfaces  in  an  arbitrary  symmetric  state  are 

zi  =  A  cos  (kx-<7)  =  real  part  of  (  ,  e"'^), 

C7.10)  ® 

Zg  =  A  cos  (kx  +  a)  =  real  peirt  of  (^‘^e  ). 

s 

The  phase  cr  is  measured  positive  downstream  from  the  a-state.  The  complex  wave 
parameter  C  has  the  value 

(7.11)  C  =  A  e"^®  =  5  -  iT). 

6 

The  deformation  Zj  of  the  upper  interface  is  associated  with  a  Laplacean 
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field  whose  streamf unction  has  the  interface  level  value  ii  =  (U/h)zi  ,  and 
whose  Intensity  at  the  lower  interface  is  reduced  by  the  factor  orse**’^.  The 
kinematic  conditions  are  that  the  interface  particles  are  moved  Jointly  by 
the  two  component  Laplacean  fields  associated  with  the  Interface  deformations. 
The  vertical  velocity  at  the  upper  interface  in  non-dimensional  units  (k=U=l) 
is 


With  the  deformations  of  the  symmetric  wave  in  (7.10)  substituted,  its  ampli¬ 
tude  becomes 


C-HC  =  i(C-aC*)/K. 

The  real  and  Imaginary  parts  of  this  equation  are  the  evolution  equations  for 
the  symmetric  wave.  They  are 


(7.12) 


(l-n,)§, 


n,  =  (1  -  a)/K 


5  =  -  (1  -  nb  )T).  Hb  =  (1  +  cr)/H 

The  paramaters  n,  and  Ub  are  the  non-dimensional  intrinsic  upstream  phase 
velocities  in  the  non-tilting  states  of  the  wave.  These  equations  are  equiva¬ 
lent  to  the  corresponding  equations  for  the  symmetric  Helmholtz  waves  in  (2.7). 
They  have  the  same  form  when  the  values  of  CtT)  are  substituted  from  (7.11)t 
namely 


•rr  In  A  =  (1  -  n,  -  &)  cot  a, 

(7.13)  ”  ® 

In  A^  =  (nb  - 1  *b)  tana. 

They  may  be  Integrated  in  precisely  the  same  manner,  to  find  the  evolution 
from  an  arbitrary  state  of  the  wave.  In  particular,  if  the  growth  rate  in 
the  states  of  stationary  phsise  (fr=0)  is  denoted  by  n  we  get 

n  tan  a  =  1  -  n.  s  &  , 

(7.14)  ® 

-  n  cot  a  =  1  -  nb  *  ftb » 
s 

precisely  as  in  (3.1)  for  the  Helmholtz  waves.  With  the  values  of  n«  and  nb 
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substltuted  from  (7*12 )  the  growth  rate  n,  with  the  proper  dimensions  kU 
restored,  becomes 

(7.15)  (iihAU)’  =  (nd/U)®  =  =  e"^*^  -  (1  -  h  )® . 

It  is  zero  for  H=l±e”**,  that  is  for  the  stationary  b-wave  in  (7*9)  and 
for  an  infinitely  long  wave  (h-»0).  Between  these  limits  n  in  (7.15)  is  real, 
and  the  waves  are  unstable.  Ihe  mode  of  maximum  growth  rate  (dn/dK=0)  has 
the  wave  number 

(7.16)  =  l-exp(-2K  )  =  0.8.  (L  =7.9d) 

m  m  m 

It  is  about  eight  times  longer  than  the  depth  of  the  shear  layer.  From  (7.15) 
its  growth  rate  is 

(7.17)  II  =  0.4  U/d  =  >4  Uk  .  (L=7.9<i) 

TD  '  m  m 

It  is  very  nearly  half  the  growth  rate  of  the  same  wave  if  the  depth  of  the 
shear  layer  shrinks  to  zero  and  becomes  a  vortex  sheet.  The  very  long  modes 
(to  the  order  of  r)  have  the  growth  rate  n  =  kU,  the  same  eis  in  a  vortex  sheet. 
From  (7»14)  the  stationary  phase  is  given  by 

tan^o  =  4^  =  or  cos  2a  =  (l-H)e’'. 

The  K  -wave  has  a  stationary  b-state  (a  =90°).  The  growing  (R=l)-fflode  is 

6  5 

45°  downstream  from  the  a-state.  The  mode  of  maximum  growth  rate  (r  =  0.8) 

m 

has  the  phase  of  cos2o  =  /O.2.  It  is  31° 45'  downstream  from  the  a-state. 

Hiis  wave  is  shown  in  fig.  7« 

Waves  which  are  very  much  longer  than  the  depth  of  the  shear  layer 

move  very  slowly  downstream  through  the  a-state  and  very  fast  upstream  through 

the  b-state,  so  their  stationary  phase  is  near  the  a-state.  To  the  order  of 

R  the  stationary  phase  of  the  very  long  waves  is  a  =)4r.  The  upper  crest 

8 

is  displaced  the  distance  d  downstream  from  the  lower  crest,  and  the  growth 
rate  (to  the  same  order)  is  kl).  The  resultant  streamfunction  at  the  upper 
interface  level,  with  higher  powers  of  r  ignored,  is 
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ti  =  A  [cos  x-a  cob(.x  +  k)  ]/h  =  /S  a  cos  (x  -  45® ) . 

8  8 

!nie  streamline  crest  is  45®  downstream  from  the  interface  crest,  the  same  as 

in  the  growing  mode  of  the  vortex  sheet.  In  the  limit  as  d-»0  the  value  of  the 

streamfunction  and  the  growth  rate  kU  become  exact  for  sdl  wave  lengths.  The 

shear  layer  behaves  approximately  as  a  vortex  sheet  for  waves  which  are  very 

much  longer  than  the  depth  of  the  layer. 

Let  us  now  consider  the  justification  of  the  linear  approximation 

which  has  been  used  in  the  above  theory.  It  is  based  on  the  assumption  that 

the  wave  amplitude  is  much  smaller  than  the  wave  length  (kA  <s^l).  In  the 

s 

short  stable  waves  the  amplitude  cannot  grow  beyond  the  b-state  amplitude. 

If  kAg^**^!  the  linear  approximation  is  good  during  the  entire  oscillation  of 

the  short  stable  wave,  and  the  linear  theory  is  valid  for  these  oscillations. 

On  the  other  hand,  a  long  unstable  wave  should,  according  to  the  linear 

theory,  from  any  state  except  the  decaying  state  of  stationary  phase  (the 

decaying  modei  move  toward  the  growing  state  of  stationary  phase  (the  growing 

mode)  where  ultimately  it  grows  at  an  exponential  rate.  After  a  certain  tim. 

the  condition,  kA  «1,  is  no  longer  satisfied,  and  the  linear  theory  is  no 
8 

longer  valid.  The  subsequent  non-linear  evolution  of  the  unstable  waves  is 
a  difficult  mathematical  problem,  but  its  main  characteristics  may  be  anti¬ 
cipated  by  qualitative  physiced  reasoning. 

Let  us  consider  the  wave  of  maximum  growth  rate  shown  in  fig.  7.  As 
this  wave  leaves  the  a-state,  the  shear  layer  becomes  periodically  inflated 
and  deflated  with  the  inflated  parts  centered  downwind  from  the  interface 
crests.  As  the  wave  approaches  the  state  of  stationary  phase,  and  the  inter¬ 
face  keeps  on  growing  at  an  exponential  rate  (according  to  the  lineary  theory), 
the  kinematics  of  the  wave  place  an  upper  limit  on  the  linear  evolution. 


namely  the  moment  when  the  deflated  parts  shrink  to  zero  in  the  middle.  Let 
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the  wave  of  maxlnum  growth  rate  be  started  from  an  a-state  with  an  amplitude 

which  is  one  per  cent  of  the  wave  length  (A  sO.OlL  ,  or  k  A  =  0.06).  Ihe 

sa  m  D  sa 

linear  theory  is  then  initially  quite  good.  If  the  lineeu?  evolution  is 

extrapolated  up  to  the  kinematic  limit,  the  wave  will  have  reached  very 

nearly  to  the  state  of  stationary  phase  at  that  tine,  so  the  kinematically 

limiting  interface  amplitude  A  is  given  by  the  formula,  A  sin  a  =  H  d,  and 

sn  601  n 

we  have 

(7.18)  A  k  =:)4K/8ino  =0.76.  (k  =0.8,  a  =31®45') 

681  81  n  to  81  81 

This  value  does  not  satisfy  the  condition  A  k<l  for  the  linear  approximation 

6 

Near  the  kinematic  limit  the  evolution  of  the  wave  is  certainly  no  longer 
linear.  The  non-linear  terns  which  have  been  ignored  are  already  becoming 
increasingly  important.  To  get  a  rough  idea  of  the  nature  of  this  non-linear 
evolution  we  shall  nevertheless  extrapolate  the  linear  evolution  uj.  to  the 
time  of  the  kinematic  limit  and  then  use  a  different  kind  of  reasoning  for 
the  subsequent  evolution. 

Let  us  first  determine  the  time  which  the  wave  needs  to  grow  from  the 

initial  a-state  of  A^^=0,01L^to  the  kinematically  limiting  amplitude 

A  k  =0.76,  or  A  _/A  =12.  This  time  (see  3.3)  is  given  by 

sm  m  snr  sa 

ch  n  t  =  12  cos  a  =  10.2  =  ch  3,  or  n  t  =  3. 
m  m  ’  m 


From  (7.17) 


n  =  )4Dk  =  w/T  , 

m  m  m’ 


(T^=L  /U). 
m  m 


Ihus,  from  an  initial  a-state  with  the  amplitude  one  per  cent  of  the  wave 
length,  the  wave  of  maximum  growth  rate  reaches  the  kinematic  limit  in  the 
state  of  stationary  phase  in  the  time  t=  (3/w)T  ,  that  is,  roughly  the  time 

ffl 

during  which  the  outer  fluid  moves  a  wave  length  in  the  symmetric  frame. 

The  subsequent  non-linear  evolution  may  be  estimated  as  follows: 


!Ihe  entire  field  of  motion  nay  be  regarded  as  the  resultant  of  all  the  partial 
fields  associated  with  the  vorticity  in  the  separated  regions  of  rotating 
fluid.  Let  us  first  consider  one  of  these  partial  fields  alone  and  estimate 
how  the  fluid  would  move  if  all  the  other  fields  were  absent.  We  represent 
this  partial  field  as  the  resultant  of  two  component  fields,  namely,  one 
component  field  coning  from  the  largest  ellipse  which  can  be  inscribed  in  the 
separated  rotational  region,  and  the  other  component  coming  from  the  remaining 
vorticity  within  the  two  asymmetrically  located  arms  which  connect  the  ellipse 
with  the  points  of  separation.  It  is  clear  that  the  major  axis  of  the  ellipse 
is  parallel  to  the  interface  tangents  at  the  nodes,  so  from  (7.18)  its  incli¬ 
nation  is  given  by 

tan  0  *  A  k  =  0.76.  (e  =  37®15'). 

sm  m 

The  minor  axis  is  one  half  of  the  normsd  distance  between  these  tangents, 
that  is 

b  =  J4d  cos  9  +  (a _/k  )  sin  9  =  Hd  cos  0  (1  +  2o  A  k _/k  )  =  0.82  d, 
nr  m  m  sm  nr  m 

or  b  =  0.82d.  The  analytical  determination  of  the  major  axis  is  not  quite 

so  simple,  but  a  careful  graphical  determination  gives  the  value  a  =  2.2d. 

The  elliptic  vortex  covers  the  nurea  nab  =  5»7d®.  The  entire  separated 

rotational  region  includes  the  fluid  within  one  wave  length  of  the  shear 

layer  and  hence  covers  the  area  L  d=  7.9d®.  The  elliptic  vortex  has  72 

m 

per  cent  and  the  connecting  arms  28  per  cent  of  the  vorticity  of  the  partial 
field. 

The  field  of  an  Isolated  elliptic  vortex  in  an  unboimded  irrotational 
surrounding  fluid  was  derived  by  Kirchhoff*.  The  elliptic  boundary  has  a 
solid  rotation  in  the  sense  of  the  vorticity  inside  with  the  anguleu:  speed 

n=  qab/(a  +  b)*.  (elliptic  vortex) 

*  Kirchhoff,  Mechanics  XX;  Lamb,  Hydrodynamics,  1932,  p.  232. 


nie  inscribed  elliptic  vortex  in  the  separated  region  of  the  shear  layer 
would  at  the  moment  of  separation  rotate  with  the  angular  speed 

Q  =  0.2  q=  0.4n/d=  n_ 

n 

if  all  the  other  vorticity  were  absent. 

Let  us  for  the  moment  ignore  the  component  field  coming  from  the 

connecting  arms,  and  consider  how  the  partial  fields  from  the  other  inflated 

separated  regions  would  move  the  elliptic  vortex.  Each  of  these  have  the 

vortex  strength  k  =  2UL  ,  and  their  action  in  the  region  of  the  elliptic 

m 

vortex  is  roughly  the  same  as  that  of  vortex  filaments  of  the  same  strength 
placed  at  their  respective  centers.  Kiese  filaments  may  be  grouped  symmetri¬ 
cally  in  pairs  whose  singular  point  coincides  with  the  center  of  the  elliptic 
vortex.  Near  this  center  all  these  partial  fields  are  pure  deformations  whose 
principal  axes  are  inclined  45®  toward  the  central  level.  This  resultant 
deformation  field  is  stationary.  At  the  moment  of  separation  the  major  axis 
of  the  elliptic  vortex  is  inclined  37°  toward  the  level.  A  short  time  later 
it  has  turned  into  coincidence  with  the  axis  of  inflow  of  the  deformation 
field.  The  deformation  field  has  therefore  at  first  very  little  direct 
influence  on  the  rotation  of  the  ellipse.  However,  it  deforms  it  into  a 
more  circular  shape,  and  thereby  causes  its  rotation  to  speed  up.  The  defor¬ 
mation  continues  until  the  major  axis  of  the  ellipse  is  vertical  and  bisects 
the  angle  between  the  stationary  axes  of  deformation.  At  that  moment  the 
ratio  between  the  major  and  minor  axis  is  a  minimum  and  the  ellipse  has  the 
fastest  rotation  from  Its  own  field.  The  deformation  now  also  contributes 
to  the  rotation  of  the  ellipse.  This  process  can  be  examined  more  quanti¬ 
tatively  ; 

The  two  members  of  the  pair  of  vortex  filaments  have  the  distance 

NL  from  the  center  of  the  ellipse.  The  field  of  this  pair  has  the  stream- 
n 
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t  =  (2UA„)  In  (BiB,). 

ni 

Using  polar  coordinates  (11,9)  from  the  center  of  the  ellipse,  we  have 

(R,  )»  =  R*  +  (NL  )*  -  2(HNL  )»  cos  20  . 

m  m 

In  the  region  of  the  ellipse  we  may  ignore  higher  powers  of  (I(/NL^)’,  and  write 

R,I^  =  (NL  )*  -R8  cos2e. 

D 

Within  this  region  the  streamf unction  has  the  ascendent 

7t  =  (4uVk  )[ein20  (R7e)-cos2e7R]/(NL_)*. 

The  vector  in  the  brackets  is  a  unit  vector  so  the  speed  of  the  field  at  the 

distance  R  from  the  center  has  the  constant  value 

V  =  2Rn  (eO-3  (n^  =  HUk„). 

m  mm 

This  is  the  speed  of  deformation  of  the  pair.  All  the  partial  fields 

have  the  same  orientation,  so  the  resultant  deformation  has  the  speed 

V  =  2Rn^S  (wN)-*  =  Rn^3  * 

At  the  moment  when  the  ellipse  is  oriented  along  the  axes  of  deformation,  its 

axes  are  changed  at  the  instantaneous  rates 

•  • 

b/b  =  -  a/a  =  n ^3, 


and  its  rotation  changes  at  the  rate 


Substituting  here, 
tion,  namely 


0  2  a-b 

n  =  J  TTT-  “m- 

for  a  rough  estimate,  the  values  at  the  time  of  separa- 


a  =  2.2  d,  b  s  0.82  d,  (L  =  n  , 

m 


*  H.  B.  Dwight:  Tables  of  Integrals  and  Other  Mathematical  Data,  Formula 
48.2.  MacMillan  Co.,  N.Y.,  1957. 
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the  ellipse  at  that  moment  heis  the  angular  acceleration  0  B0.3n^.  Let  us 

assume,  again  for  a  rough  estimate,  that  the  acceleration  has  this  constant 

value  irtille  the  ellipse  turns  from  Its  Inclined  position  at  the  moment  of 

separation  (9  =  37°15')  to  the  vertical  position  (6=90^),  and  evaluate  the 

corresponding  time  Interval  t.  It  Is  given  by 

0.15(n  t)*+n  t  =  (52.75A80)«  =  0.92, 
in  D 

or  t-0.82/n  =0.26T_.  If  we  Ignore  the  acceleration  and  let  the  ellipse 
n  in 

rotate  with  the  constant  Initial  angular  speed  n  It  would  turn  the  same 

n 

angle  In  the  time  Interval  0.29 T  . 

zn 

It  remains  to  consider  the  component  field  from  the  vortlclty  In  the 

two  arms  which  connect  the  ellipse  with  the  points  of  separation.  We  recall 

that  these  include  about  28  per  cent,  that  is,  each  aomi  about  l4  per  cent  of 

the  total  vortlclty  in  the  separated  inflated  region.  This  component  field 

Is  roughly  the  same  as  that  of  a  pair  of  vortex  filaments  of  the  same  strength 

placed  at  the  centers  of  vortlclty  of  the  arms.  At  the  moment  of  separation 

these  filaments  have  roughly  the  distance  2.4d  =  0.3L  from  the  center  of  the 

n 

ellipse,  and  their  center  line  is  Inclined  about  22HP  toward  the  center  level. 
In  the  region  Inside  the  ellipse  the  field  of  this  pair  is  approximately  a 
pure  deformation  with  axes  45^  inclined  toward  the  center  line.  At  the 
distance  S  from  the  center  of  the  ellipse  the  speed  of  this  deformation  field 

is 

V  =  0.l4*2Ka^(0.3it)-®  w  Hn 

This  field  has  about  the  same  strength  in  the  region  of  the  ellipse  as  the 
resultant  deformation  from  all  the  other  separated  vortlclty  concentrations. 
The  orientation  of  this  field  is  such  that  it  tends  to  slow  down  the  rotation 
of  the  ellipse.  Furthermore,  whereas  the  deformation  field  from  the  other 
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vortlclty  concentrations  is  quasl-stationary,  this  latter  field  from  the 
connecting  arms  follows  the  ellipse  In  its  rotation.  The  two  deformation 
fields  effect  the  rotation  of  the  ellipse  in  opposite  sense,  and  since  both 
effects  are  rather  small,  every  elliptic  vortex  rotates  roughly  as  if  all 
the  vorticlty  outside  the  ellipse  were  absent.  It  turns  with  the  angular 
velocity  n  from  the  inclined  position  at  the  time  of  separation  to  the 

ID 

vertical  position  during  the  time  the  outer  air  moves  about  one  quarter  of 
a  wave  length.  The  ellipse  is  deformed  toward  a  more  roiinded  shape  during 
this  transition.  We  note  that  these  conclusions  resemble  the  results  from 
Bosenhead's  computations  for  the  vortex  sheet. 

8.  The  statically  stable  shear  layer.  -  We  shall  now  examine  a  layer 
with  constant  shear  between  unbounded  irrotational  outer  layers  of  different 
densities  with  the  lighter  fluid  (p=  Pi )  above  and  the  heavier  fluid  (p=  Oa) 
below  the  shear  layer.  The  density  distribution  in  the  shear  layer  may  be 
assigned  in  different  ways.  In  this  section  we  consider  the  theoretically 
simplest  model  where  the  shear  layer  is  homogenous  with  a  density  which  is 
the  arithmetic  mean  of  the  outer  layer  densities.  In  this  model  the 
boundaries  of  the  shesu:  layer  have  the  same  static  stability 

(8.1)  )4s  =  (p,  -  Pi  )/(p,  +  Pi ). 

This  system  with  arbitrary  densities  in  the  outer  layers  was  investigated  by 
Q.  I.  Taylor*  in  1916.  The  quasi-symmetric  system  (s<  1)  was  investigated 
later  by  Goldstein**.  The  quasi-symmetric  shear  layer  is  governed  by  a 
Richardson  number  which  has  the  value 

(8.2)  4  =  )48gd/U*. 

*  Taylor,  Q.  I.,  19315  Effect  of  Variation  in  Density  on  the  Stability  of 
Superposed  Streams  of  Fluids.  Proc.  Roy.  Soc.  A  132,  pp.  499-523. 

Goldstein,  S.,  1931:  On  the  Stability  of  Superposed  Streams  of  Fluids  of 
Different  Densities.  Proc.  Roy.  Soc.  A  132,  pp.  524-548. 
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We  shall  examine  the  behavior  of  symmetric  waves  in  this  quasi- 
symmetric  shear  layer.  Similar  to  the  statically  neutreil  shear  layer,  the 
deformations  of  the  botmdaries  generate  Laplacean  fields  associated  with 
the  displaced  vorticlty  in  the  shear  layer.  !nie  Interface  level  value  of 
the  component  field  associated  with  the  upper  deformation  Z]  is  (U/h)zi . 
However,  the  static  stability  across  the  Interface  now  gives  rise  to  sliding 
vorticlty  along  the  deformed  interface  whose  change  is  governed  by  the  dynamic 
condition  in  (2.4),  namely. 


Let  ^1  denote  the  interface  level  value  of  the  streamfunction  of  the  Laplacean 
field  associated  with  this  sliding  vorticlty.  In  terms  of  the  streamfunction 
the  dynamic  condition  has  the  form 

(8*3)  “  (HsgA)^^  • 

Let  us  find  out  under  what  conditions  this  sliding  vorticlty  wave  is  stationaiy 
and  neutral  in  the  symmetric  frame.  Ihe  individual  change  is  then  equal  to 
the  convective  change,  so  the  dynamic  condition  becomes  tj  =  OJsgAUfej.  The 
sliding  vorticlty  wave  is  in  phase  with  the  interface  deformation  with  maximxim 
value  at  the  crest.  Hie  field  strength  is  such  that  the  sliding  vorticity  of 
the  interface  filament  is  exactly  reversed  by  the  overweight  while  the  fila¬ 
ment  moves  from  crest  to  trough  in  the  symmetric  frame.  This  is  the  wave 
mechanism  of  the  neutral  gravity  wave,  so  we  shall  call  this  part  of  the 
symmetric  wave  a  gravity  wave.  Hie  second  part,  which  comes  from  the  dis¬ 
placed  vorticity  of  the  shear  layer,  has  the  field  strength  (D/k)zi  of  the 
Bayleigh  wave.  We  shall  call  this  part  of  the  symmetric  wave  a  shear  wave. 

The  gravity  wave  is  momentarily  stationary  and  neutral  in  the  symmetric 
frame  if  it  is  in  phase  with  the  shear  wave  with  an  amplitude  which  is  u 


Fig.  8  Stability  diagrom  for  a  strotified  sheer  layer  (after  Goldstein) 
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tlaea  the  ^tar  wave  anplltude.  Hxia  condition  may  be  aatiafied  In  an  arbi¬ 
trary  tilting  atate  of  the  aymmetrio  wave.  However,  if  the  entire  aynmetrlc 
wave  ahall  be  atationary  and  neutral,  thla  condition  muat  be  obeyed  by  both 
the  gravity  wave  and  the  ahear  wave  almultaneoualy,  which  meana  that  the 
aynmetrlc  wave  nmat  be  non-tilting.  When  the  ahear  wave  ia  stationary  in  a 
non-tilting  state,  the  upper  and  lower  fields  will  jointly  move  the  inter¬ 
face  deformations  intrlxisically  upstream  with  the  speed  U,  so  the  kinematic 
condition  at  the  upper  interface  of  the  stationairy  wave  is  (see  7>12) 


(«=  e"*) 


The  statically  stable  shear  layer  accordingly  has  two  stationary  neutral  v/aves, 
namely, 


(8.4) 


Stationary  a-wave:  n, (1 +  u)  =  1, 


n,  =  (1  -  er)/K, 


Stationary  b-wave:  nk(l+u)  =  1.  n>  =  (l+cr)/H. 
The  wave  numbers  of  these  stationary  waves  for  any  given  Richardson  number 
are  shown  by  the  two  lines  in  the  stability  diagram  in  fig.  8.  Both  lines 
have  the  same  asymptote,  u=k-1.  Ihe  a-line  starts  from  the  origin  with 
the  tangent  slope  H .  Near  the  origin  the  stationary  a-wave  is  very  nearly 
given  by 


u  =  )4h,  or  Bg/kVP  =  1.  (H=kd«l) 

If  the  depth  of  the  shear  layer  shrinks  to  zero  this  formula  is  exact  and 
gives  the  stationary  neutral  Helmholtz  wave  (2.2)  which  meirks  the  transition 
between  the  long  stable  and  the  short  unstable  waves  in  the  statically  stable 
vortex  sheet. 


The  b-line  starts  with  unit  slope  from  the  stationary  neutral  Rayleigh 
wave  which  marks  the  transition  between  the  long  unstable  and  the  short 
stable  waves  in  the  statically  neutral  shear  layer.  This  makes  it  plausible 
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that  the  shear  layer  has  unstable  waves  in  the  spectral  Interval  between  the 
two  lines,  as  indicated  in  the  figure. 


9.  Ihe  normal  modes  of  the  shear  layer.  -  In  an  arbitrary  state  the 
symmetric  wave  has  the  wave  elements 

Zi  =  A  cos  (kx-a)  =  Real  part 

=  (U/h)  a  cos  (kx-a)  =  Real  part  (Ze^^). 

The  complex  wave  parameters  Q  and  Z  have  the  values 

C  = 


(9.2) 


=  ?  -  ill. 

-ie 


Z  =  (U/k)  a  e  =  X-  iY. 


The  amplitude  A  of  the  gravity  wave  is  weighted  by  the  factor  (U/h).  This 

makes  A/A  the  amplitude  ratio  of  the  sliding  vorticity  to  the  shear  vorti- 
s 


city,  thus 


/Q  _  _A _ sliding  vorticity  amplitude 

~  sKear  vorticity  amplitu(ie  ' 

The  change  of  the  gravity  wave  is  governed  locally  by  the  overweight  of  the 
sloping  interface  filaments,  so  we  may  apply  this  condition  in  a  local  frame 
which  moves  with  the  filament.  The  individual  change  is  the  local  change  in 
this  frame.  In  the  units  U=k=l  the  dynamic  condition  (8.3)  for  the 
symmetric  wave  (9.1)  is  therefore  Z=  (u/H)iC,  which  has  the  real  and  imaginary 


parts 

X  =  (m/h)71. 

or 

• 

A 

• 

COB  0  -  0  A  sin  9  =  u  A  sin  o, 

(9.4) 

• 

A 

s 

-Y  =  (ii/h)5. 

or 

sin  9  9  A  cos  9  =  -  u  A  cos  a. 

s 

Elimination  of  0  gives  the  growth  rate  of  the  gravity  wave,  namely, 
(9*5)  A/A  =  (u/r)  sin  (o- 9). 

The  gravity  wave  grows  when  it  is  upstream  from  the  shear  wave  with  the 
sliding  vorticity  maximum  in  the  hill  upwind  from  the  interface  crest. 
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Ellfflination  of  A  from  (9*4)  givos  the  Intrlnalc  downwind  phase  veloeltj 
of  the  gravity  wave  relative  to  the  fluid  at  the  interface,  namely, 


9  a  -  (p/r)  cos  (o-  6). 

Since  the  fluid  moves  downstream  with  unit  speed  in  the  symmetric  frame,  the 
phase  velocity  In  this  frame  Is 


(9.6)  9  =  1  -  (p/r) cos  (a- 0). 

If  the  wave  has  a  state  of  stationary  phase,  9  =  9=0,  the  growth  rate  n 

and  the  amplitude  ratio  r  in  this  state,  obtained  from  the  ratio  and  the  sum 

8 

of  the  squares  of  (9*5)  and  (9*6),  are 

n  =  tan  (o  -  6  ), 

8  S 

(9.7) 


r|  =  uV(l  +  n*). 

The  behavior  of  the  shear  wave  in  an  arbitrary  state  is  not  quite  so 
simple  because  the  boundaries  of  the  shear  layer  are  moved  by  all  four 
Laplacean  fields.  However,  if  the  shear  wave  has  the  same  phase  as  the  gravity 
wave  (0=0),  It  behaves  at  that  moment  as  a  Rayleigh  wave  in  a  statically 
neutral  shear  layer,  only  with  the  field  intensities  augmented  by  the  factor 
1  +  r.  In  particular,  its  phase  velocity  in  a  state  of  o=  0  (see  7.14)  is 
(9.8)  o  =  1- (l  +  r)(n,coo’o  +  nfcein®a).  (o=0) 

In  the  non-tilting  states  both  the  gravity  wave  and  the  shear  wave  are  instan¬ 
taneously  neutral,  and  they  move  through  these  states  with  the  downwind  phase 
velocities 


0«  =  1  -  p/r,  =  n  tan  0.,  i  =  1  -  iv  (1  ♦  r,  )  »  n  tan  o_, 

0  6 

(9.9)  . 

0*  =  1  -  p/r*  =  -  n  cot  0  ,  o*  =  1  -  n* (1  +  r* )  =  -  n  cot  o  . 

8  6 

The  first  parts  of  these  formulas  apply  to  non-tilting  states  with  arbitrary 
amplitude  ratios.  If  the  amplitude  ratios  have  the  values  characteristic  of 
a  normal  mode,  the  non-tilting  phase  speeds  of  the  mode  are  related  to  the 
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growth  rate  and  stationary  phase  of  the  mode  by  the  formulas  In  (5.1}t  vdileh 
have  been  repeated  in  (9«9)«  Both  the  gravity  wave  component  and  the  shear 
wave  component  are  stationary  and  have  the  same  growth  rate  In  the  state  of 
stationary  phase  of  the  mode  (see  9*7)t  so  the  results  derived  In  section  5 
apply  to  both.  In  particular,  the  principle  of  constant  energy  transport  of 
the  component  waves  In  (3*5 )« 


Sa*  =  S  A»j 

8t  A 


iA‘  , 

a  sa* 


give  the  relation  between  amplitude  ratio  and  phase  velocities  In  an  arbitrary 
state  of  the  mode,  namely, 

{9*10)  (0/o)r*  =  (0t/a,  )r*  =  (Ok/o»)rf  =±r,rfc. 


or 

(9.11)  0,/Oa  *  Ofc/Sk  air^/r, . 

These  formulas  apply  to  both  unstable  modes  (n^>0)  and  stable  modes  (n^<0). 

In  both, the  phase  velocity  ratio  is  Inverted  from  one  non-tilting  state  to  the 
other.  This  means  that  the  stable  mode  has  a  tilting  state  wAiere  both  compo¬ 
nent  waves  move  with  equal  speeds  (d=±o).  This  state  of  the  stable  mode  is 
the  nearest  equivalent  to  the  state  of  stationary  phase  of  the  unstable  mode 
(zero  speed  of  both  waves).  We  shall  denote  it  by  the  subscript  m.  The 
amplitude  ratio  in  the  state  of  equal  wave  speeds  (from  9.10)  is  the  geometric 

mean  of  the  ratios  in  the  non-tilting  states  r^ar.rk.  The  phase  of  this  state 

n 

is  obtained  from  the  phase  equations  (3.^)  for  the  component  waves,  namely, 
tan  0  B  tan  0  thnt,  tan  a  =  tan  o  thnt. 


Their  ratio  gives 

(9.12)  tan  0/tan  o  *  9^/^  * 

The  tangents  of  the  phase  angles  of  the  component  waves  have  the  same  ratio  in 
all  states  of  the  normal  mode.  The  time  derivatives  of  the  phase  equations  in 
a  state  of  equal  speeds  have  the  ratio 
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(9.13)  coe*0|jj/co8*a^  =  ±  Oii/9b • 

Tht  product  of  those  equstioas  gives 

sin  29  T  sin  2o  =  2  sin  (9  To  )  cos  (•«*«)  *  0» 

n  n  Bunn  mm 

I I  ^  I  I  *  means  that 

l®ml  +  1‘^al  “ 

Thust  adiether  the  component  waves  of  a  stable  mode  move  in  the  same  or  in 
opposite  directions,  the  sum  of  the  distances  they  have  moved  from  a  non- 
tiltixig  state  to  the  state  of  equal  speeds  is  equal  to  the  distance  between 
the  non-tilting  states  (a  quarter  of  a  wave  length).  Thus,  in  the  state  of 
equal  speeds  the  ratio  in  (9*12)  becomes 

(9.1^)  =  ®o**®n*  ^^m"**m^ 

Finally,  the  value  of  the  equal  phase  speeds  "is  obtained  from  the  phase  equa¬ 
tion  (3*6),  namely, 

0  =  9,  cos* 9  ■••^sin*9  =  (4, -0k)cos*9  +0*, 
m  m  m  m 

with  the  value  of  cos  9^  substituted  from  (9*13) t  which  gives 

(9.15)  4  . 

These  formulas  define  the  state  of  equal  speeds  of  the  stable  mode  in  terms 
of  its  phase  velocities  in  the  non-tilting  states. 

Ihe  symmetric  waves  in  the  statically  stable  shear  layer  have  an  infinite 
number  of  non-tilting  states  depending  upon  the  value  of  the  amplitude  ratio 
of  the  component  waves,  which  may  be  assigned  arbitrarily.  To  construct  an 
arbitrary  non-tilting  state  as  a  resultant  of  normal  modes,  we  need  two  pairs 
of  modes.  The  frequency  equation  for  this  system  must  therefore  be  a  quadratic 
equation  in  n*.  We  know  already  that  the  equation  has  two  distinct  roots, 
n*  >  0,  for  the  wave  lengths  of  the  stationary  a-wave  and  the  stationary  b-wave 
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in  (8.4),  so  it  must  have  the  form 

[n*  -  1  + (1  +  li)]Ca*  “  1  +  n*  (1  +  ti)]  =  Bn®. 

The  value  of  B  Is  found  from  the  two  limiting  systems 

U-»0:  (n®/u+ n,  )(n*/pi+ n» )  *  (n*/u)(B/»i)  =  0, 

s  -•0;  n*  -  n*(B-  n,  -  n^  -t-  2)  (n.  -  l)(nk  -  1}  =  0, 

The  first  of  these  gives  the  frequencies  -  n®  s  m®  s  of  the  two  stable 
modes  in  the  resting  three-layer  system,  and  shows  that  B  is  Independent  of 
u.  One  mode  in  the  statically  neutral  shear  layer  (s=0)  is  the  Hayleigh 
mode  (see  7.14)  with  n®  s  -  (n,  -l)(nh  -1).  The  second  mode  has  therefore 
n®  =  -l.  It  is  stable  v/ith  the  frequence  m=l  and  the  amplitude  ratio  rs-1. 
The  shear  vortlclty  at  the  deformed  boundaries  is  everyidiere  neutralized  by 
an  equal  amoimt  of  sliding  vortlclty  along  the  boundaries.  Ihe  mode  is  a 
neutral  wave  idiieh  moves  downstream  with  the  fluid.  Die  coefficient  of  -n® 
in  the  frequency  equation  is  the  sum  of  the  roots,  so  B  has  the  value 
B  =  n,+n^-3  -  (n,  -  l)(nb-l)  =  -(n,  -2)(nk-2). 

With  this  value  substituted,  the  frequency  equation  takes  the  final  form 
(9.16)  Cn®  -  1  +  n,  (l  +  u)]Cn®  -  1  +  n>  (1  +  u)]  =  -  n®(n,  -  2)(n»  -  2). 
n® =  0  on  the  a-line  and  the  b-line  in  the  K,u-8tablllty  diagram  (fig.  8)  and 
for  no  other  wave  lengths,  so  only  on  these  lines  can  n*  change  sign.  Uie 
sign  is  Indeed  changed  for  one  of  the  n®-roots  when  these  lines  are  crossed, 
for  we  have  seen  that  the  a-line  near  the  origin  marks  the  transition  from 
long  stable  to  short  vmstable  Helmholtz  waves,  and  the  b-llne  marks  the 
transition  (for  ^=0)  from  the  long  unstable  to  the  short  stable  Hayleigh 
waves.  It  follows  that  one  family  of  modes  is  unstable  in  the  spectral 
interval  between  these  lines  and  stable  in  the  rest  of  the  spectrum.  The 
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Growth  roto  of  tho  onotoblo  modo 
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Other  family  la  stable.  The  growth  rate  of  the  imatable  modes  In  units  of 
U/d  is  shown  for  selected  values  of  the  Richardson  number  in  fig.  9. 

With  the  growth  rate  (or  frequency)  of  the  mode  known,  its  non-tilting 
amplitude  ratios  are  obtained  from  the  products  of  the  a-  and  b-equations  in 
(9.9),  namely, 

^  Bk  *  (1  -  u/r,  ) (1  -  a/r* )  =  -  n? 

(9.17) 

*  Cl-n,(l  +  r,)][l-nfc(l  +  rt)]  =  -n? 

The  non-tilting  amplitude  ratios  in  turn  determine  the  non-tilting  phase 
velocities  of  the  mode.  These  finally  determine  the  other  characteristic 
features  of  the  mode.  All  the  properties  of  the  normal  modes  are  therefore 
known  in  principle  from  their  growth  rate  (or  frequency).  An  arbitrary 
disturbance  may  be  represented  as  the  resultant  of  normed  modes,  so.  its  evolu¬ 
tion  is  known  in  principle  as  the  resultant  evolution  of  the  modes. 

To  become  familiar  with  the  physical  mechanism  of  this  evolution  we 
shall  examine  it  in  a  wave  where  the  mathematical  theory  is  very  simple, 
namely,  the  (nka2)-wave. 

10.  The  (nb  =  2)-wave.  -  This  wave  is  a  little  longer  than  the  wave  of 
maximum  growth  rate  in  the  statically  neutral  shear  layer  (7.16).  It  has  the 
wave  number 

(10.1)  n*  =  2:  H  =  Jid  +  e"*^)  =  0.74  (n,  =  2/h -2  =  0.7  ). 

The  value  of  n,  is  closer  to  O.7I,  but  the  slightly  smeller  value  in  (10.1)  is 
accurate  enough  for  a  rough  estimate. 

The  frequency  equation  (9.16)  gives  the  growth  rates  for  the  two 
(n^ s  2)-modes,  namely. 


(10.2) 


(i): 

(ii); 


n*  =  l-n,(l  +  u)  s  0.3 -pun,, 
n*  =  l-n^d  +  u)  =-(l  +  2w). 
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3%e  first  mode  is  unstable  below  the  a-llne,  and  stable  for  greater  static 

t 

stability.  The  second  Is  stable.  Let  \is  examine  the  evolution  of  these 
modes  from  their  non-tUtlng  states: 

(l)  w£tab^e_m£d^:  !Ihe  non-til  ting  amplitude  ratios  of  this  mode 
are  obtained  from  (9.17)t  namely 

1  -  n,  (1  +  u)  =  -  Cl  -n,  (1  +  r,  )][l  -  2(1  +  r»)3  »  -  (1  -  g/r,  )(1  -  g/r^). 

These  are  satisfied  by  s  u,  rb  =  0.  The  phase  velocities  In  the  non-tilting 
states  are  accordingly 

•  e 

r,  =  u:  0,  =  0,  o,  =  0.3 -un,, 

(10.3) 

r*  =  0:  9b  =  •,  Ob  =-l. 

The  gravity  wave  remains  permanently  stationary  in  the  non-tilting  a-state 
and  Is  absent  from  the  b-state.  The  shear  wave  behaves  formally  like  a 
Helmholtz  wave  (see  2.8),  but  the  working  mechanism  Is  different.  The 
Helmholtz  wave  moves  upstream  through  the  b-state  with  unit  speed  due  to  the 
kinematic  constraint  at  the  interface  nodes.  The  shear  wave  in  (10.3)  moves 
upstream  through  the  b-state  with  unit  speed  because  it  is  here  a  simple 
Hayleigh  wave  with  the  sliding  vorticities  momentarily  absent.  In  the  a-state 
the  mechanism  is  more  similar:  In  the  absence  of  static  stability  the  Helmholtz 
wave  moves  downstream  through  the  a-state  with  unit  speed  due  to  the  sliding 
vorticity  convergence  at  the  downstream  nodes.  The  statically  neutrad  shear 
wave  moves  downstream  with  the  speed  0.3,  resulting  in  inflation  of  the  shear 
layer  at  the  downstream  nodes.  The  static  stability  (and  the  sliding  vorti¬ 
city  it  excites)  counteracts  the  downstream  propagation  of  both  waves,  and 
becomes  the  dominating  factor  when  sufficiently  strong.  The  entire  evolution 
of  the  shear  wave  in  (10.3)  may  be  represented  graphicedly  by  a  diagram  very 
like  the  Helmholtz  diagram  in  fig.  2,  with  straight  a,-  and  Ob -lines. 

The  shear  wave  has  a  stationary  a-state  for  j.=  n^^  - 1  =  0.4,  where  the 
(nb  =  2)-llne  meets  the  a-line  (see  fig.  10).  For  weaker  static  stability 


Fig. to  The  unetoble  (n^>2)-mode 


-61- 


the  shear  vortlclty  dominates,  and  the  shear  wave  moves  downstream  throuj^  the 
a-state.  As  it  leaves  the  a-state  it  begins  to  grow  but  slower  than  the 
statically  neutral  Rayleigh  wave,  because  the  sliding  vortlclty  field  (which 
remains  behind  In  a  growing  a-state)  retards  the  growth.  As  the  static  sta¬ 
bility  Is  Increased  from  zero,  it  acts  In  two  ways  to  reduce  the  growth  of  the 
shear-wave:  (l)  It  reduces  the  downstream  speed  through  the  a-state  and  there¬ 
fore  gives  the  shear  wave  a  stationary  phase  closer  to  the  a-state  which  is 
less  favorable  for  efficient  growth.  (11)  It  retards  the  grovirth  directly 
because  the  crest  of  the  shear  wave  Is  downstream  from  the  sliding  vortlclty 
maximum  (see  fig.  10). 

For  u  >0.4  the  sliding  vorticities  are  strong  enough  to  move  the  shear 
wave  upstream  through  the  a-state.  It  is  a  stable  wave  which  moves  progres¬ 
sively  upstream,  while  the  gravity  wave  performs  standing  oscillations  in  the 
a-state.  The  shear  vortlclty  and  the  sliding  vortlclty  act  on  the  changes  of 
the  interface  amplitude  in  opposite  sense  in  these  waves.  If  a =1.85  the  two 
actions  are  equal,  and  the  shear  v/ave  Is  a  neutral  wave  which  moves  upstream 
with  constant  unit  speed  and  constant  amplitude.  For  weaker  static  stabilities 
(in  the  interval  0.4 <n <1.85)  the  shear  vortlclty  dominates,  so  the  shear  wave 
decays  and  speeds  up  from  the  a-state  to  the  b-state,  and  it  slows  dovm  again 
and  grows  as  the  wave  moves  on  to  the  next  a-state .  For  greater  static  stability 
the  sliding  vorticities  from  the  overweight  dominates,  and  the  shear  wave 
changes  speed  and  amplitude  in  opposite  rhythm. 

(11 )  8.t£bLe_modB :  The  amplitude  ratios  in  the  non-tilting  states 
of  this  mode  are  given  by  the  equations 

l  +  2u  =  [l  -  n,  (1 +r,  )][l  -  2(1  +  r^)]  =  (1  -  u/r,  )(l  -  n/r* ), 
which  give  r,  =  - 1,  r*  =  -  (1  +u).  The  mode  is  stable  and  moves  progressively 
downstream,  passing  the  nou-tilting  states  vdth  the  extreme  phase  velocities 


(10.4) 
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•  • 

-r,  ®  1:  •*  “  1  +  Ui  =  li 

-r*  =  1  +  u:  0*  ■  (1  +  2(i)/(l  +  u),  c»  =  l  +  2u. 

In  the  statically  neutral  shear  layer  this  wave  has  the  constant  amplitude 
ratio  rs  -1.  Ihe  sliding  vorticities  cancel  the  shear  vorticitles  so  both 
fields  drift  downstream  with  the  fluid.  In  statically  stable  shear  layers 
the  mode  has  the  same  general  behavior  for  all  values  of  u:  resultant 

field  is  zero  in  the  a-state,  so  the  shear  wave  drifts  downstream  with  the 
fluid  throu^  the  a-state,  while  the  gravity  wave  moves  downstream  faster 
than  the  fluid,  driven  by  the  overweight  at  the  downstream  nodes.  After  having 
left  the  a-state  the  sliding  vorticity  maximum  is  downstream  from  the  Interface 
crest  where  it  is  augmented  by  the  overweight,  so  the  gravity  wave  grows  emd 
slows  down.  The  shear-vorticlty  tends  to  augment  the  interface  amplitude,  but 
this  tendency  is  more  than  compensated  by  the  stronger  opposite  action  of  the 
sliding  vorticity.  The  shear  wave  decays  and  speeds  up  after  leaving  the 
a-state. 

As  the  gravity  wave  slows  down  and  the  shear  wave  speeds  up  they 
approach  the  state  of  equal  speed  with  the  gravity  wave  at  the  maximum  dis¬ 
tance  downstream  from  the  shear  wave.  Ihe  amplitude  ratio  heis  here  the  value 
r^  =  r, r>  =  /I  +  u.  The  phase  of  the  component  waves  is  given  by 

tan’Gn,  =  cot’a®  =  9,/o;  =  l  +  Ui  or  cos  2oin  =  -  cos 20in  =  w/(2  +  ii), 

and  the  equal  speed  of  the  waves  is  (2  +  3a)/(2  +  u).  The  (u  =  2)-mode  has  a 
very  simple  state  of  equal  speeds  (9^  =  60°,  651  =  2,  r5,  =  /5). 

After  the  state  of  equal  speeds  is  passed  the  sheaur  wave  has  the  faster 
speed.  Both  waves  arrive  together  in  the  b-state,  the  shear  wave  with  the 
maxlmtim  speed  1-)-2m.,  the  gravity  wave  with  the  minimum  speed  (1  2u)/(l u). 

The  evolution  from  an  arbitrary  state  of  the  (nks2)-wave  may  now  be 
found  as  the  resultant  evolution  of  the  normal  modes.  For  example  the  evolu- 
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tion  from  an  Initial  a-state  with  the  interface  amplitude  AgsAg^  and  no 
eliding  vortieity  along  the  interface  (AgsO)  may  be  represented  as  the 
resultant  of  the  a-state  of  the  unstable  mode  with  the  amplitude  Agi  »  Agg/(1  +  u) 
and  the  a-state  of  the  stable  mode  with  the  amplitude  Agg  =  Agaii/(1  +  ii).  The 
shear  nave  component  of  this  resultant  wave  moves  downstream  through  the 
a-state  as  a  simple  Saylel^  wave  with  the  speed  1  -ngs  0.3>  The  momentarily 
absent  gravity  wave  is  immediately  excited  with  the  negative  vortieity  maximum 
at  the  downstream  nodes.  In  terms  of  the  evolution  of  the  modes  the  appear¬ 
ance  of  this  sliding  vortieity  may  be  interpreted  as  the  dorastream  displace¬ 
ment  of  the  negative  vortieity  of  the  stable  mode,  while  the  positive  vortieity 
of  the  unstable  mode  remains  behind  in  the  a-state. 

11.  The  modes  on  the  a-  and  b-lines.  -  All  the  mathematical  formulas 
for  the  symmetric  waves  in  the  shear  layer  are  symmetric  in  the  subscripts  a 
and  b,  so  the  formulas  for  the  modes  on  the  a-  and  b-lines  are  the  same  with 
the  subscripts  interchanged. 

On  the  a-llne  the  roots  of  the  frequency  equation  (9.16)  are 
(i)  s  0,  r,  =  u, 

(11.1)  n,(l  +  u)=l 

(ii)  n®  =  1  -  n^d  +  u)  -  (o«  -  2)(nk  -  2). 

The  roots  on  the  b-llne  are  obtained  by  interchanging  a  and  b.  Consider  first 

the  (n®  =  0)-mode. 

(i)  The  (n^  s  0)-mode;  On  the  a-llne  this  mode  has  a  stationary  neutral 
a-state  with  the  amplitude  ratio  r^  z  u.  If  started  from  the  b-state,  the  mode 
moves  toward  the  stationary  a-state  and  approaches  this  state  with  an  asymptotic 
linear  growth  rate.  The  amplitude  ratio  for  the  b-state  of  this  mode  cannot  be 


obtained  from  (9>17}t  which  here  gives  n^  s  9,  so.  =0.  However,  it  will  be 
shown  in  the  next  section  (see  13.4)  that  the  non-tilting  phase  velocities  of 
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the  gravity  wave  in  a  normal  mode  are 

(11.2) 


1  _  n*  -  1  +  n.  (1  + 1*) .  A  _  n*  -  1  +  n*  (1  ♦  n) 

- sn - ’  **  - - 


These  values  satisfy  the  upper  equation  in  (9.17)f  and  they  reduce  to  the 
values  for  the  (nks2)-ittode8  in  section  10.  With  these  veilues  substituted, 
the  equations  (9.9)  give  the  non-tilting  amplitude  ratios  of  the  mode,  namely. 


(11.3) 


(11.4) 


2  -n, 
- 3-; 


£l 

u 


2-n, 


■L. 


li  1  +  +  n  u  1  +  un^  +  n* 

The  b-state  amplitude  ratio  of  the  (n^  s  0)-mode  on  the  a-line  is 


rt  _  2  -Hb  _  ^  nt  (1  +  u)  - 1  _  _  tan  fl 

u  1  +  im*  1  +  luib  ® 


(a-line) 


Since  r,  s  u.  this  formula  gives  the  b-state  phase  velocity  ratio  as  well 
(see  9.11|12).  The  b-state  amplitude  ratio  is  zero  in  the  (nb  =  2)-mode.  It 
is  positive  in  the  shorter  modes,  increaising  to  the  asymptotic  limit  r^  s  u  in 
the  very  short  modes.  It  is  negative  in  the  longer  modes,  decreasing  to  the 
asymptotic  limit  r^  -  *■>&  in  the  very  long  modes  (ant  ->1). 

From  (11.2)  and  (11.4)  the  b-state  phase  velocities  in  this  mode  are 


=  nb(l^a)-l  ^ 

(11.5)  (a-line) 

-ab=  n>(l.rb)-l. 

UUlb  -t- 1  Itflb  + 1 

On  the  a-line  n* (!  +  (*)>!,  or  tinb+l>2-nb,  so  the  shear  wave  moves  upstream 
in  all  the  (n*  =  0)-modes  on  this  line.  It  moves  with  unit  speed  in  the 
(nbs2)-fflode  where  the  sliding  vorticity  is  absent  from  the  b-state.  In  the 
longer  modes  its  upstream  speed  is  greater  than  unity,  but  less  than  the 
speed  of  the  Bayleigh  wave  because  it  is  retarded  by  the  negative  sliding 
vorticity.  The  gravity  wave  on  the  other  hand  moves  downstream  in  the  long 
modes,  faster  than  the  fluid,  driven  by  the  overwei^t  at  the  downstream  nodes. 
After  leaving  the  a-state  the  gravity  wave  grows  from  the  overweight  at  its 
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crests  and  slows  down.  Also  the  shear  wave  grows,  but  slower  than  the  Rayleigh 
wave  because  the  growth  Is  retarded  by  the  sliding  vorticity.  Both  waves  move 
on  in  opposite  directions  toward  asymptotic  approach  to  the  stationary  a-state 
which  they  never  quite  reach.  Ihe  gravity  wave  retains  a  residual  phase 
upstream  and  the  shear  wave  a  residual  i^ase  downstream  from  the  a-state. 

Since  the  energy  transport  is  constant,  their  growth  continues  indefinitely 

•  •  ♦  •  . 

at  a  linear  rate  (AsA^dk  Ags-Ag^Ok). 

In  the  shorter  (nk<2)-modes  both  component  waves  move  upstream  through 
the  b-state.  Ihe  speed  of  the  sheeir  wave  is  less  than  unity,  but  greater  than 
the  speed  of  the  Rayleigh  wave  due  to  the  positive  sliding  vorticity.  The 
gravity  wave  is  driven  upstream  by  the  overwei^t  at  the  upstream  nodes.  In 
modes  near  the  (ub  =  2)-mode  this  component  is  very  weak  and  moves  very  fast. 

In  the  shorter  modes  with  greater  static  stability  it  is  stronger  and  has  a 
slower  speed.  In  the  (nb-l)-mode  for  example  (see  7*9)  the  gravity  wave  has 
the  b-state  velocity  -9b  =&•  -1  =  0.77>  Since  rb<|i  in  the  short  modes 
(see  11.4),  the  gravity  wave  is  leading  and  the  shear  wave  trailing  all  the 
way  from  the  b-state  to  the  asymptotic  a-state.  The  gravity  wave  gi*ows  from 
the  overweight  at  its  crest.  'Die  shear  wave  grows  at  a  slower  rate  than  the 
Rayleigh  wave  because  its  growth  is  retarded  by  the  sliding  vorticity.  In 
the  asymptotic  a-state  the  mode  has  the  residual  downstream  phase  ratio 
a/9=  u/rb  >1.  The  waves  have  here  the  same  relative  position  as  in  the  longer 

s  e  •  • 

modes,  and  the  asymptotic  linear  growth  As  -AbOb,  As=  -Ag^Ob* 

On  the  b-llne  the  (n°  s  0)-mode  has  a  stationary  neutral  b-state  with 
the  amplitude  ratio  rb  =  p..  From  (11.4)  its  a-state  amplitude  ratio  is 

(11.6)  iJLslJliii-s  (b-line) 

U  1  pn,  1  pn*  6„  tan  9 

Here  Ob  (1 |i)  <1,  or  2  -  n,  >1  pn, ,  so  r,  >|i  in  this  mode.  From  (11.5)  its 


a-atate  phaae  velocities  are 


§  ,  i-iUiSL*  1-JL>0, 

2  -  a.  2  - 


(11.7) 


•  l-n,(l  +  tt)  2-1 


1  1  -«■  lU], 


(b-line) 


-1  =  1  - (1  +  r,  )  >0. 


Both  components  move  downstream  through  the  a-state  with  the  shear  wave 
leading  and  the  gravity  wave  trailing.  Ihe  shear  wave  moves  slower  than  the 
fluid,  but  faster  than  the  Kaylei^  wave  due  to  the  positive  sliding  vorticity. 
The  gravity  wave  moves  intrinsically  upwind  with  a  speed  less  than  unity 
(r,  '>u)  driven  by  the  overweight  at  the  upstream  nodes.  As  it  trails  behind 
the  shear  wave,  it  grows  from  the  overweight  at  its  crests,  while  the  growth 
of  the  shear  wave  is  retarded  by  the  sliding  vorticity.  In  the  asymptotic 
b-state  the  mode  has  the  residual  upstream  phase  ratio  Ok  -  9)/(>^  -  a)  =  rn/u.  >  1, 

e  e  e  •  . 

and  the  growth  continues  at  a  linear  rate  (A  =  A,9^,  AgsAg^^a,}. 

If  the  wave  length  is  much  shorter  than  the  depth  of  the  shear  layer, 
the  fields  from  the  lower  boundary  have  no  appreciable  effect  at  the  upper 
boundary  and  vice  versa.  In  the  short  modes  on  the  asymptotic  (a,b}-line  the 
shear  wave  and  the  gravity  wave  are  both  stationary  and  neutrad  in  an  arbitrary 
tilting  state  and  have  the  amplitude  ratio  u. 

The  (n^  =  0}-mode  on  the  a-  and  b-lines  hais  a  state  of  equal  speed 
(9  +  Os H a)  downstream  from  the  a-state.  The  phase  and  speed  of  this  state 
are  given  by  (9.1^)  and  (9. 15). 

£tab^e_mod^:  On  the  a-llne  this  mode  is  defined  by  the  second 
n®-root  in  (11.1).  With  this  value  substituted  (11.2)  gives  the  b-state 
phase  velocity  of  the  gravity  wave  in  this  mode,  namely, 

9b  =  2-  n,  =1- u/rb,  ao  u/rb  =  n,  - 1  =  -  u/(l  +  u). 

Its  b-state  amplitude  ratio  and  phase  velocities  are  therefore 
(11.8)  -Tb  a  1  +  U,  9b  =  (1  +  2u)/(l  +  u),  Ob  =  l  +  MJlb. 


(a-line) 
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The  amplitude  ratio  Is  the  same  as  In  the  stable  (nbs2)-fflode  (see  10.4). 

The  two  modes  have  a  rather  slmllair  behavior.  In  both  the  component  waves 
move  downstream  through  the  b-state  faster  than  the  fluid.  In  the  present 
mode  the  phase  velocity  difference  Is 

(11.9)  »*  -  =  uCnfcd  +  u)  -  ll/(l  +  n)  >0*  (a-llne) 

Also  here  the  shear  wave  Is  leading  the  gravity  wave  downstream  from  the 
b-state.  It  slows  down  and  grows  from  the  b-state  to  the  a-state,  speeds  up 
and  decays  from  the  a-state  to  the  b-state  and  so  on,  while  the  gravity  wave 
changes  speed  and  amplitude  In  opposite  rhythm. 

The  a-state  amplitude  ratio  and  phase  velocities  are  obtained  from 
the  formula  (9.11),  namely. 


(11.10) 


£L  =  2L  =  ilJi2L(i  +  ^)  =^=:iS£i>0, 


ft*  1  +  2H 


Oa 


tano 


(a-line) 


which  gives 

(11.11)  -r,  =  .LtitL,  e,  =l-n/r,,  o,  (a-line) 

l-t-n^u  »i  +  l 

The  a-state  amplitude  ratio  has  the  asymptotic  lower  limit  -r.  =)^ln  the  very 
long  modes  (n^u-*!).  Increasing  to  one  in  the  (nks2)-fflode  and  on  to  the 
asymptotic  upper  limit  u  in  the  very  short  modes. 

The  behavior  of  the  stable  mode  on  the  b-line  is  obtained  from  equa¬ 
tions  (11.8-11)  by  Interchanging  the  subscripts.  Also  here  the  amplitude 
ratio  Is  negative  with  -  r,  =  1  |i  and  -  r^  increasing  from  unity  in  the  stati¬ 
cally  neutral  shear  layer  to  the  asymptotic  upper  limit  u  in  the  very  short 
modes.  The  gravity  wave  therefore  moves  progressively  downstream  faster 
than  the  fluid,  leading  the  shear  wave  from  the  a-state  to  the  b-state: 

Oa-Oa  =  n[l  -  Ha  (1  +  ^)]/(l  +  >  0.  (b-llne) 

It  grows  and  slows  down  from  the  a-state  to  the  b-state,  speeds  up  and  decays 
during  the  next  quarter  period  and  so  on,  while  the  shear  wave  changes  speed 
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and  anplltude  In  opposite  rhythm,  precisely  as  in  the  stable  mode  on  the 
a-line.  Hie  behavior  of  the  gravity  wave  is  clearly  understood  in  terms 
of  the  overweight  at  its  crests  and  nodes.  The  shear  wave  responds  to  the 
joint  action  of  shear-  and  sliding  vortlclties,  with  the  sliding  vorticity 
dominating  in  all  states  of  this  mode. 

The  stable  mode  has  the  same  general  characteristics  on  the  (n^sS)- 
line  and  on  the  a-  and  b-llnes:  It  has  a  negative  amplitude  ratio  wi.th 
minimum  absolute  value  in  the  a-state  and  maximum  value  (-r^  ^  l-*-|i)  in  the 
b-state.  Its  component  waves  move  progressively  downstream.  They  move 
faster  than  the  fluid  in  all  modes  except  the  shear  wave  near  the  a-state  in 
the  long  (nb  >2)-modes.  The  gravity  wave  is  the  leading  component  from  the 
a-state  to  the  b-state.  It  is  the  dominating  component  in  all  states  of  the 
stable  mode.  Its  eliding  vorticity  governs  both  the  propagation  and  the 
growth  and  decay  of  the  shear  wave,  so  it  moves  intrinsically  downstream 
through  the  fluid  with  pulsations  opposite  to  the  action  from  its  own  field. 
It  is  plausible  that  the  stable  mode  has  the  same  behavior  in  all  parts  of 
the  spectrum. 

12.  The  unstable  mode  in  the  shear  layer.  -  We  have  seen  that  the 
statically  stable  shear  layer  has  an  unstable  mode  in  the  spectral  interval 
between  the  a-  and  b-lines  in  the  H,u-diagram.  Near  the  a-line  this  mode 
has  an  asymptotic  linearly  growing  state  with  the  shear  wave  a  residual 
phase  downstream  from  the  a-state.  Near  the  b-line  the  unstable  mode  has 
an  asymptotic  linearly  growing  state  with  the  shear  wave  a  residual  phase 
upstream  from  the  b-state.  At  intermediate  points  the  unstable  mode  heis 
an  exponentially  growing  state  of  statloneury  phase  with  the  phase  of  the 
shear  wave  downstream  from  the  a-state  Increasing  continuously  from  zero  on 
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th«  a-llne  to  90°  on  the  b-llne.  Ihe  gravity  wave  haa  the  proper  phase 
upatreaai  from  the  shear  wave  to  make  it  grow  at  the  rate  of  the  mode  (see  9<7)« 
The  unstable  mode  has  a  rather  simple  kinematic  structure  on  the 
straight  line  which  is  the  common  asymptote  of  the  a-  and  b-llnes,  namely, 

(12.1)  1  +  u  s  h:  n,  (l  +  ii)  =  1 -or;  n*  (1  +  u)  =  1 +  ar. 

On  this  line  the  non-tilting  phase  velocities  of  the  gravity  wave  in  the 
normal  modes  (11.2)  are 


0,  =  (af-n®)/(2-n,  )  *ntanea, 

(12.2) 

-0b  =  (ar+ n*)/(2  -  nb  )  =  n  cot  0a. 

Its  stationary  phase  in  the  unstable  mode  is  given  by 

(12.3)  tan»08  =  ~ ^  ^  2h -  1  ] -  [or»  4:n» {2k-  Ij}  ^ 

(or  +  n® )  (2  -  a,  )  er[n*  +  2k  - 1  ]  +  [or*  +  n*  (2h  - 1 )  ] 

where  n*  is  the  positive  root  of  the  frequency  equation 

(12.4)  n* -Of*  +  a*(2  - n,  )(2  - nb )  =0. 

In  shear  layers  with  very  weak  static  stability  (u<l)  this  wave  has  very 
nearly  the  wave  number  k  =  1,  and  its  frequency  equation  is 

I*<1;  n*-Qr*  +  n*(l-o(*)  =  (n*  +  l)(n* -or*)  =  0.  (k  =  1) 

The  vmstable  mode  has  here  the  growth  rate  n  =  or=e"^,  and  from  (9.7)  the  ampli¬ 
tude  ratio  ra  =  lid +  0(^)“’^<1.  The  very  weak  gravity  wave  has,  from  (12.2),  the 
stationary  phase  tan  ©g  =  (1 -Qr)/(l  +  or),  or  tan(45°  -  0b)  =  or=  n.  From  (9*7)  the 
stationary  phase  of  the  shear  wave  is  therefore  Og-  45°.  The  very  weak  gravity 
wave  is  15>2*’  upstream  from  the  shear  wave,  but  gives  no  significant  contribu¬ 
tion  to  its  growth.  The  growth  of  the  shear  wave  is  caused  by  the  shear  field, 
centered  at  the  opposite  interface,  fdilch  acts  with  its  full  strength  a  in  this 
mode. 


In  sheaur  layers  with  very  strong  static  stability  (|i>l)  the  mode  on 
the  a,b-asymptote  is  very  much  shorter  than  the  depth  of  the  shear  layer 
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(x»l),  80  the  component  fields  are  reduced  to  Insignificant  strength  at  the 
level  of  the  opposite  interface.  The  wave  at  each  interface  behaves  as  if  the 
wave  at  the  other  interface  were  absent.  The  frequency  equation  (12.4)  for 
these  very  short  modes  is 

U^^l:  n*  +  4n*  =  n*(n*  +  4)  =  0.  (k»1) 

One  mode  is  a  stable  neutral  wave  with  the  amplitude  ratio  r  s  which  moves 
intrinsically  upstream  through  the  fluid  with  unit  speed,  so  it  is  stationary 
in  the  symmetric  frame.  The  second  mode  is  a  stable  neutral  wave  with  the 
amplitude  ratio  r  =  - u  which  moves  intrinsically  downstream  through  the  fluid 
with  unit  speed,  so  it  moves  downstream  with  the  speed  and  frequency  m  -  2  in 
the  symmetric  frame. 

For  shear  layers  with  dominating  static  stability  (u>l)  a  good  approxi¬ 
mation  for  the  positive  root  is  obtained  by  writing  the  frequency  equation 

(12.4)  in  the  form 

= ^ _ 

(2  -  n,  )(2  -  n*  )  +  n®  [e''(2H  -  l)]®  +  (nne’*)®  -  1 

For  Hi 2,  the  term  e*‘(2H -  1)  >21,  so  practically  no  error  is  made  if  the  last 

term  in  the  denominator  (-1)  is  ignored.  The  approximation  is  even  better 

when  also  the  middle  term  is  ignored,  giving  the  growth  rate 

(12.5)  U^l;  n  w  Ke-V(2H-1).  (Ki2) 

Even  in  the  shear  layer  with  static  stability  and  sheeu:  of  equal  strength 
(|isl)  this  formula  gives  the  growth  rate  of  the  unstable  mode  on  the  a,b- 
asymptote  with  an  accuracy  better  than  one  per  cent.  From  (12.3)  the  station¬ 
ary  phase  of  the  gravity  wave  in  this  mode,  to  the  same  accuracy,  is  given  by 

(12.6)  uil:  cos  20,  =  Of  .  in^2) 

For  strong  static  stabilities  (u>^l)  the  growth  rate  is  n  =  )4ar,  the  gravity 
wave  has  the  stationary  phase  0,  s)((i(->6Qr),  and  from  (9.7)  the  shear  wave  has 


Fig,l2  The  unstoble  (/i=  k  <•  I )  -  mode 
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the  stationary  phase  OgS  amplitude  ratio  of  the  mode  Is  rg  =  |i« 

With  this  stationary  phase  the  shear  wave  would  grow  from  Its  own  field,  If 
acting  alone,  at  the  rate  or.  The  sliding  vortlcity  of  the  very  strong  gravity 
wave  with  its  maximum  or  upstream  from  the  interface  crest  reduces  this  growth 
rate  by  fifty  per  cent. 

In  conclusion  we  shall  as  a  typlcad  example  exhibit  the  unstable  (k=2) 
mode  in  the  shear  layer  with  static  stability  and  shear  of  equal  strength 
(ii.  =  l).  It  is  shown  in  fig.  12.  From  (12.5)  its  growth  rate  is 
U=l,  H=2:  n  =  2/(3e®)  =  0.09  mc=  0.l8U/d. 

From  (12.6)  and  (9.7)  its  component  waves  have  the  stationary  phase 
cos  20s  =7/(9e^)  =  0.104,  Og  =  42°,  Og  =  47.2®, 
and  the  amplitude  ratio  rg= 0.996*  If  the  temperature  difference  of  the 
outer  layers  is  5°  Kelvin  and  their  relative  translation  2U=10ms”^,  the 
shear  layer  between  them  must  have  the  depth  ds600mto  make  it  a  (u^^l)- 
layer  (see  8.2).  The  (H=2)-niode  in  this  layer  has  the  wave  length  l880m. 

In  the  asymptotic  state  of  stationary  phase  its  amplitude  doubles  every 
7.7  minutes  for  as  long  as  the  linear  theory  is  a  valid  approximation. 

13.  Shear  layer  between  bounded  outer  layers.  -  To  satisfy  the  kine¬ 
matic  conditions  of  zero  vertical  motion  at  the  rigid  boundaries  of  the  outer 
layers  (as  in  4.1),  the  streamf unction  of  the  field  associated  with  the 
periodic  vorticity  field  (u*  -u")  at  the  upper  boiuidary  of  the  shear  layer 
heis  above  and  below  that  interface  the  values 

t*  =  84k(h  -  z)/sh  kh,  (z>0) 

t’=  8hk(h  +  d  +  z)/ah  k(h  +  d).  (z<0) 

Here  h  denotes  the  depth  of  the  outer  layers  and  d,  as  before,  the  depth  of 
the  shear  layer.  The  tangential  components  of  the  field  at  the  interface  are 
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s  i^ikcthkh, 

-u*a  tikcthkCh  +  d). 

Th*  corresponding  vortlclty  field  at  the  interface  is  therefore 

u*-u*=  (2K/d)ti. 

Here  we  have  introduced  the  non-dimensional  parameter 
(13.1)  H  =  )4kd[cthkh  + cthk(h +  d)]. 

At  the  lower  boundary  of  the  shear  layer  the  field  is  reduced  to  the  value 
ail  t  where 

(13«2)  Of  =  shkh/sh  k(h  +  d). 

For  waves  which  are  much  shorter  than  the  depth  of  the  outer  layers  these 
parameters  reduce  to  the  earlier  imbounded  layer  values,  namely, 

kh>l:  H  =  kd,  or  =  e”**. 

From  the  dynamic  condition  in  (7*2)  the  shear  vorticity  field  associated 
with  the  deformation  Zj  of  the  upper  interface  is 

u*-u"=  (2U/d)zx  =  (2K/d)ti, 

so  the  streamf unction  of  this  field  has  the  Interface  level  value  (U/h)zx, 
as  in  (7«3)  for  the  unbounded  system.  The  sliding  vorticity  field  along  the 
interface  changes  in  accordance  with  the  dynamic  condition  in  (8.3),  namely, 

SO  the  streamf unction  of  this  field  has  the  change 

Let  us  apply  these  conditions  to  an  arbitrary  state  (9.1)  of  a  symme¬ 
tric  wave.  If  the  boundary  conditions  are  satisfied  at  the  upper  interface 
for  the  symmetric  wave,  they  are  automatically  satisfied  at  the  lower  inter¬ 
face,  so  it  is  sufficient  to  consider  the  upper  conditions.  The  non-dimen¬ 
sional  dynamic  condition  at  the  upper  interface,  namely. 
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\  dt  3x/ 


♦i  = 


with  the  wave  elements  substituted  from  (9>1)«  becomes 

Z-iZ=  (m/h)C*  (Cynamlc  condition) 

The  kinematic  condition  at  the  upper  interface,  sue  in  (7«12),  is 

^  +Zi/H  +  or(t|  -z,/h)]. 

For  the  symmetric  wave  in  (9*1)  this  condition  becomes 

C-iC  =  Z+  -  a(Z*  +  CVh).  (Kinematic  condition) 

Substitute  here  the  value  of  the  shear  wave  parameter  C  from  the  dynamic  con¬ 
dition.  The  resulting  equation 

Z-Z+ (l  +  u)(Z-(yZ*)/H  =  -2iZ  +  i(Z+erZ*)/H 
has  the  real  and  imaginary  parts 

X- [l-n,(l  +  u)3x  «  (n,-2)Y,  n,  =  (1  -  Qr)/H 

if-  Cl-n*(l  +  u)]Y  =  -  (ttk  -2)X.  nk  =  (l  +  £r)/H 

They  are  the  evolution  equations  for  the  gravity  wave  component  of  the 
symmetric  wave. 

In  a  state  of  stationary  phase  (9s  9^= const)  only  the  amplitude  of 
the  wave  changes  with  time,  so  the  evolution  equations  (15>5)  for  such  a 
state  are 


(13.5) 


A  -  [l -n,  (1  +  u)]A  =  (n,  -  2)A  tan  Sat 
A  -  [l  -  n*  (1  +  u)  ]A  =  -  (n*  -  2)A  cot  9g. 

These  are  different  forms  of  the  same  equation.  It  has  solutions  of  the 
fozmi  A~e^^  where  the  growth  rate  n  is  a  root  of  the  simultaneous  quadratic 
equations 

n®  - 1  +  n,  (1  +  u)  =  (n, -2)ntan9a  =  94(n, -2), 

(13.4) 

n®  - 1  +  n*  (1  +  |i)  s  -  (n*  -  2  )n  cot  9a  =  9*  (n*  -  2 ) . 

Elimination  of  9g  gives  the  bl-quadratlc  growth-rate  equation  in  (9.16), 
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whloh  accordlztgly  has  the  aame  form  for  bounded  and  for  unbounded  outer 
layers.  The  wave  has  a  state  of  stationary  phase  63  with  a  real  growth 
rate  n  In  the  spectral  Interval  where  the  growth  rate  equation  has  a  posi¬ 
tive  root  n^  >0.  The  sum  of  the  roots  is  never  positive,  so  the  growth 
rate  equation  has  a  positive  root  if  the  constant  term  is  negative,  that  is 
in  the  spectral  interval  where 

(13.5)  Cn,(l  +  u)-l][nfc(l  +  u)-l]<0.  (Instability) 

This  condition  is  satisfied  in  the  spectral  interval  between  the  stationary 
neutral  a-wave  and  the  stationary  neutral  b-wave  (8.4).  In  this  Interval  the 
shear  layer  has  one  pair  of  unstable  modes,  one  growing  and  one  decaying  at 
the  rate  n,  idiere  n  is  the  real  positive  root  of  the  growth  rate  equation. 

As  In  section  9  shall  regard  this  pair  as  the  two  states  of  stationary 
phase  of  a  single  mode. 


14.  The  instability  criterion  for  the  bounded  shear  layer.  -  The 
condition  for  instability  in  (13»5)  may  be  written 

(l4.l)  h/(1 -Of)  >  1  +  |i>H/(l  +  a),  (Instability) 

where  h  and  or  are  the  non-dimensional  depth  parameters  in  (13.1,2),  namely, 

H  =  )4kd[cth  kh  +  cthk(h  +  d)3. 


(14.2) 


a  a  shkh/sh  k(h  d). 


For  unbounded  outer  layers  K  =  kd,  ors  e~’*,  and  the  condition  gives  instability 
in  the  spectral  interval  between  the  a-  and  b-lines  in  Goldstein's  stability 
diagram  in  fig.  8.  For  bounded  outer  layers,  choosing  for  the  moment  the 
length  unit  }^dsl,  we  have 

H  _  ehk(h  ■♦•2)chkh-i-  shkh  chk(h-f2) 
k(l±o)  shkhC8hk(h  ■t-2)^shkh]  ’ 


so  the  lower  limit  in  (l4.l)  is 
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w  ^  sh2k(h-»l)  ohk(h  ■»!) 

kCl-fer)  ”  2  8hkhshk(h-«-l)ehk  *  shkhchk  * 

and  the  upper  limit  Is 

H _ ah2k(h-H)  ^  ahk(h-t-l) 

k(l-ar}  2  sh  kh  ch  k(h  l)shk  shkhshk  * 

Substitution  of  these  In  (l4.l)  with  arbitrary  length  imlts  restored  gives 

Instability  of  the  bounded  shear  layer  in  the  Interval 

(14.3)  cthkh +cthHkd>2(l +u)Ad>cthkh+ thHkd.  (Instability) 
The  upper  limit  is  a  stationary  neutral  a-wave  and  the  lower  limit  is  a 
stationiary  neutral  b-wave. 

If  the  depth  of  the  shear  layer  shrinks  to  zero,  no  wave  satisfies 
the  lower  b-wave  limit.  ^Qie  vortex  sheet  has  no  stationary  b-wave  whether 
the  layers  it  separates  are  bounded  or  not.  However,  as  kd-*0,  the  a-wave 
limit,  using  (8.2),  may  be  written 

cthkh-sgAtf*  =  (Hkd)”^ -cth>6kd- -kel/6-0, 
so,  as  shown  earlier  in  section  4,  the  vortex  sheet  between  layers  of  finite 
depth  is  unstable  for  waves  shorter  than 

(14.4)  jijj  =  (egAO®)thkh<l. 

Here  is  the  Richardson  number  of  the  Helmholtz  wave  in  (4.3). 

In  the  statically  neutral  (u=0)-Bhear  layer  no  wave  satisfies  the 
upper  a-wave  limit  in  (14.3).  l^e  condition  for  a  stationary  b-wave  is  here 
Stationary  b-wave;  cthkh  =  (34kd)"^  -  thyikd.  (u=0) 

The  longer  modes  are  unstable  and  the  shorter  stable.  This  condition  is 
represented  graphically  in  a  (kd,  kh) -diagram  by  a  line  which  leaves  the 
origin  with  the  slope  h  s  yi  d  and  has  ^he  vertical  aisymptote  kd  s  of  the 
stationary  b-wave  in  the  shear  layer  between  unbounded  outer  layers  (see  7.9). 
If  the  statically  neutral  shear  layer  fills  more  than  one  half  the  space 
between  the  outer  rigid  boundaries  (d>2h)  it  has  no  stationary  b-wave  and 


no  unstable  modes 
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Waves  which  are  very  much  longer  than  the  total  depth  between  the 
rigid  boundaries,  k(d-)-2h)^l,  obey  the  quasi-static  approximation  with 
good  accuracy  (see  section  5)*  ^e  instability  condition  (14.3)  for  the 
long  quasi-static  waves  in  the  boimded  shear  layer  may  be  written  in  the 
form 

(14.5)  1  >  sgh/U*=|ih  >  l-2h/d.  k(d  +  2h)<ci 

It  is  determined  by  the  Richardson  number  for  the  quasi-static  Helmholtz 
waves  (5 >6)  and  the  ratio  between  the  depths  of  the  shear  layer  and  the 
outer  layers.  If  shear  layer  has  no  unstable  quasi-static  waves. 

The  (ujj= l)-shear  layer  has  one  quasi-static  mode  with  a  stationary  neutral 
a-state  and  an  asymptotic  linearly  growing  a-state.  The  other  mode  is  stable. 
If  the  shear  layer  occupies  more  than  one  half  of  the  space  between  the  outer 
boundaries  (d>2h)  and  its  static  stability  has  the  value  ujj®l-2h/d,  it 
has  one  qiiasl-static  mode  with  a  stationary  neutral  and  a  linearly  growing 
b-state.  Finally  shear  layers  which  occupy  less  than  one  half  of  the  space 
between  the  outer  boundaries  have  quasi-static  instability  if 

15.  The  statically  stable  Th-shear  layer.  -  We  shall  now  examine  a 
quasi-symmetric  shear  layer  which  has  continuous  variations  of  velocity  and 
density  toward  asymptotic  values  far  above  and  below  the  central  level, 
namely,  . 

(15.1)  U  =  Uoth(2z/d),  p  =  Poexp[-s  th®(2z/d)]. 

We  shall  call  this  a  Th-shear  layer.  It  resembles  the  discontinuous  Taylor- 
Qoldstein  constant -shear  layer  in  section  8,  and  we  shall  show  that  it 
exhibits  similar  dynamic  characteristics.  Far  from  the  central  level  the 
riuid  moves  in  opposite  directions  with  the  same  constant  speed.  The  fluid 
ut  t)ic  central  level  has  the  maximum  vorticity  2Uo/d,  so  d  is  the  depth  of 
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the  "equivalent"  conatant-ahear  layer.  The  Ih-veloolty  profile  In  (1^.1) 
la  very  cloae  to  the  error-profile  which  developa  by  vortlclty  dlffiialon 
from  a  statically  neutral  vortex  aheet  (aee  fig.  6). 

The  ni-shear  layer  in  (15.1)  has  no  static  stability  at  the  central 
level.  The  variation  of  the  density  from  far  below  to  far  above  the  central 
level  is  a  small  fraction  of  the  density  at  the  eentrad  level,  so  the  laiyer 
la  quasl-aymmetrlc.  Let  the  aubscrlpts  1  and  2  denote  the  asymptotic  limiting 
densities  above  and  below  the  central  level.  Hiey  have  the  values 

Pi  =  Coe"®  a  Pod  -  a), 

(5«1) 

Pa  B  Poe®  B  Pod  +  s), 

80  the  stability  parameter  s  has  the  same  meaning  as  in  the  constant  shear 
layer  (see  8.1),  namely. 


(15.2)  6  =  (pa  -  Pi  )/(Dt  +  Pi )  «1. 

In  the  following  we  shall  use  units  of  length  and  time  such  that 
HdaUo  al,  and  Introduce  the  abbreviated  notations 

T  =  th(2z/d)  =  thz. 


(d=2) 


S  =  sech(2z/d)  =  sechz. 

These  are  connected  by  the  relations 

T'  a  S*  =  1-T“,  T"  =  2SS'  =-2TT'  b-2TS*. 

Here  and  In  the  following  the  accent  denotes  differentiation  with  reference 


to  height. 

The  static  stability  of  the  density  stratification  in  (15*1)  is 
represented  by  the  logarithmic  density  gradient,  namely, 

(15.5)  o  =  -dno)'  =  (Iner)'B  58T*S®. 

We  note  that  the  stratified  Th-shear  layer  in  (15.1)  has  maximum 
static  stability  and  maxlnrum  vortlclty  gradient  at  the  levels 

o'  =  0;  T®  ■  >4, 


U’"  B  0: 


T®  B  >4 


(z  B  O.hhd) 
(z  B  o.35d) 
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Both  11*  inside  the  boundaries  of  the  equivalent  constant-shear  layer  (z8>id). 
We  must  expect  some  difference  between  the  dynamic  behavior  of  the  7h-shear 
layer  and  the  constant-shear  layer  on  this  account. 

l6.  Ihe  wave  equations  for  the  Ih-shear  layer.  -  We  now  introduce  into 
the  Ih-shear  layer  an  arbitrary  two-dimensional  small  amplitude  wave  disturb¬ 
ance  in  the  planes  of  the  basic  flow  whose  velocity  field  ir(x,z)  is  repre¬ 
sented  by  a  streamfunction  as  in  (1.2).  The  resultant  motion 
satisfies  the  dynamic  equation 

DV 

~  =  (af  +  5)7(p  +  p), 

where  or  and  p  denote  the  undisturbed  static  fields  when  the  wave  is  absent, 
and  the  barred  symbols  denote  the  local  variations  of  these  fields  due  to  the 
wave.  We  linearize  this  equation,  substitute  for  the  static  pressure  gradient 
from  the  static  equation,  eliminate  the  dynamic  pressure  p  by 

vectoriad  differentiation.  The  result  is 

0  =  7 X  o(D3^t  +  j^/or)  a»7x  +  jp/or). 

The  approximation  indicated  in  the  last  member  on  the  rl^t  is  the  quasi- 
symmetry  approximation  vdilch  ignores  the  kinematic  effect  of  the  density 
gradient.  (See  section  4  and  also  the  theory  of  the  isothermal  atmosphere 
in  Chapter  II,  section  7*) 

For  the  two-dimensional  incompressible  flow  the  inertial  term  in 
the  vortlcity  equation  may  be  transformed  as  follows 
DV  DV  DQ 

7X^  =  7X  C.2:+QX  V)  =  _*  .  (Qa7XVaUM  +  7Xv) 

~Dt~Dt'*~  Dt  - -  - 

The  fluid  is  incompressible  so  the  local  change  of  the  specific  volvune 

ccmes  from  the  vertical  displacement  Zg  of  the  fluid  particles  from  their 
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undisturbed  levels.  Its  value  to  the  linear  approximation,  using  the  nota¬ 
tion  in  (13*3) I  Is  therefore 

a/a  B  -  (In  or) '  z,  a  -  azg» 

With  these  substituted,  the  quasi-symmetrle  vortieity  equation  above  takes 
the  form 


(16.1) 


D  -“a 


dz. 


where 

(16.2)  ‘I  =  98  + 9g 

is  the  vortlcity  field  of  the  wave.  The  vortlcity  equation  (l6.l)  is  the 
dynamic  equivalent  of  the  corresponding  condition  at  a  statically  stable 
interface  between  fluid  layers  of  different  densities  (see  2.4),  and  it 
represents  the  same  physical  principle:  The  vortlcity  of  the  fluid  particle 
is  changed  in  the  sense  of  its  overweight  idiich  is  proportional  to  the  static 
stability  and  the  slope  of  the  isosteres.  The  vortlcity  is  the  sum  of  the 
shear  'vorticlty  of  the  basic  flow  U*  and  the  vortlcity  which  is  added  by 
the  wave  motion.  As  In  the  constant  shear  layer  the  added  periodic  vorti¬ 
clty  field  q  in  the  wave  may  be  regairded  as  the  sum  of  two  component  fields 
as  indicated  in  (l6.2).  One  component  q^  comes  from  the  vertical  displace¬ 
ment  of  the  shear  vortlcity  U*  of  the  basic  flow.  We  call  this  component 
the  shear  wave.  The  second  component  q^  is  the  part  of  the  vortlcity  which 
is  changed  by  the  overweight  and  will  be  called  the  gravity  wave.  In  a 
statically  neutral  shear  layer  (osO)  the  gravity  wave  is  absent,  and  the 
dynamic  eqiiation  (l6.l)  gives  the  shear  vortlcity  change 

DU' 


(16.3) 


Dt 


»  -  U"w  a  -  U"  ■ 


Dt  - - W  ’ 

or,  when  Integrated,  the  shear  vorticlty  itself 
(16.4)  q^  a  -U"Zg. 


-Bo¬ 


lt  is  positive  (in  the  sense  of  the  basic  shear)  for  displacements  away 
from  the  central  level,  and  negative  for  displacements  toward  that  level. 
Evidently  the  same  relation  (16*4)  holds  for  the  shear  vortlclty  in  the 
statically  stable  shear  layer  so,  with  (16.3)  substituted,  the  dynamic 
vortlclty  equation  (16.1)  gives  the  change  of  the  vortlclty  field  In  the 
gravity  wave,  namely. 


(16.5) 


Dt  ^  ^8''  Dt 


Dq 

Dt 


og 


®®s 
8x  * 


Let  VIS  now  assume  that  the  shear  wave  and  the  gravity  wave  are  in  phase 
(or  In  opposite  phase)  at  all  levels,  with  the  amplitude  ratio 
(16.6)  r(z)  =  Qg/q*.  (Equal  phase) 

However,  the  wave  may  have  an  arbitrary  tilt.  In  such  a  state  of  the  wave 
the  extreme  values  of  the  vortlclty  coincide  at  all  levels  with  the  crests 
and  troughs  of  the  shear  wave,  so  the  gravity  wave  Is  momentarily  neutral 
in  a  state  of  equal  phase.  LetQ(z)  denote  the  intrinsic  upstream  speed  of 
the  gravity  wave  through  the  fluid  in  a  state  of  equal  phase.  In  a  frame 
which  moves  locally  with  the  gravity  wave  the  fluid  moves  downstream  with 
the  speed  Q  and  the  gravity  wave  has  at  the  moment  no  local  change,  so 
(16.5)  for  a  state  of  equal  phase  becomes 


dz. 


ag-£  =  -2££2£  , 

®  dx  U"  dx  ’ 


dx 


which  gives  the  intrinsic  upstream  phase  velocity 


(16.7) 


In  the  last  term  on  the  right  we  have  substituted  from  (15*1'’3) 


UU11  2  T*T' 


1.58g 


l*5»»i 


so  the  Richardson  nvunber  u  has  the  same  meaning  as  in  the  constant  shear 
layer  (see  8.2).  The  phase  velocity  formula  for  the  gravity  wave  in  (16.7) 
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resembltts  the  correspondla^  formula  (9*6 )  for  the  constant  shear  layer 
whiohi  for  a  state  of  equal  phase,  gives  the  intrinsic  upstream  speed  ^/r. 
However,  the  gravity  wave  in  the  Th-shear  layer  has  50  per  cent  more  vorti- 
city  than  a  wave  with  the  same  relative  speed  in  the  constant -shear  layer. 
The  reason  for  this  is  at  least  partly  that  n<I]^  at  all  levels  in  the 
Th-shear  layer,  so  the  gravity  wave  needs  less  overweight  to  have  the  same 
relative  speed. 

To  derive  a  corresponding  phase  velocity  formula  for  the  shear  wave 

we  must  specialize  further  to  a  non-tilting  state  of  equal  phase  where  also 

the  shear  wave  is  instantaneously  neutral.  Let  C(z)  denote  the  upwind 

intrinsic  speed  of  the  shear  wave  in  a  non-tilting  state.  Tbe  kinematic 

condition  for  such  a  state  may  he  written 

Pzs  c  8* 

Dt  dx  dx  * 

which  gives  the  intrizisic  upwind  phaise  velocity 


(16.8) 


C 

U 


Uz. 


0  q.  >Jq  q 


This  formula  resembles  the  corresponding  formula  for  the  constant -sheeur 
layer  in  (9«8),  which  gives  the  intrinsic  upstream  speed  (l-)-r}n,^t  for 
the  shear  wave,  where  the  factors  n^^  constants.  However  the  corres¬ 
ponding  factor  2S^j/q  in  (l6.8)  is  in  general  a  function  of  height.  The 
Ih-shear  layer  is  therefore  a  more  complex  dynamic  system. 

The  conditions  for  a  stationary  neutral  non-tilting  state  of  a  wave 
in  the  Th-shear  layer  are  from  (16.7-8) 


(16.9)  r  =  1.5u,  q  =  -7»t  =  (2  +  3»*)S®t.  (G=C=U) 


The  first  condition  makes  the  gravity  wave  stationary.  The  second  condi¬ 
tion  makes  the  shear  wave  stationary.  When  the  streamfunction  hsis  a  sinu¬ 
soidal  variation  in  the  x-directlon  this  condition  is  a  differential  equation 
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of  the  hypergeometrie  type«  It  has  been  studied  by  S.  V.  Garcia*  who  was 
able  to  find  all  the  stationary  neutral  waves  in  this  !ni-shear  layer.  We 
shall  derive  Garcia's  (still  unpublished)  solutions  of  (16.9)  in  section 
20  below.  But  first  we  shall  consider  a  few  simple  waves  which  may  be 
derived  with  less  mathematical  labor. 


17.  Some_£tationar^_i»ve8_j£jbhe^_^;^ear^_la22£»  ”  form  of  the 
stationary  wave  equation  in  (I6.9)  makes  it  plausible  that  it  may  have  solu¬ 
tions  whose  height  variations  are  simple  functions  of  T  and  S,  say 

(17.1)  ♦  =  AnT“s“cos  kx. 

For  the  moment  m  and  n  are  arbitrary  non-negative  constants.  IVom  (I6.9) 
the  corresponding  vortlclty  q  of  this  wave  is  given  by 

(17.2)  q/*  =  k* -♦"/♦. 

Logarithmic  height  differentiation  of  (17.1)  gives 

t'/t  «  nT*A  +  nS'/S  a  n/T-  (m  +  n}T, 


and  a  second  differentiation  gives 

♦"/♦  =  Cn/T- (m  +  n)T]' +  [n/T- (m  +  n)T]* 

=  -  C(m  +  n)(m  +  n+l)  -  n(n-  l)/T*  ]S®  +  m®. 

With  this  value  substituted.  (17.2)  becomes 

q/(S*t)  =  (m  + n)(m  + n  + 1)  -  n(n- 1)/T®  +  (k*  -  m*)/S®. 

We  shall  only  consider  waves  of  this  type  for  \idilch  m  a  k.  so  the  last  term 
on  the  right  is  zero.  Introducing  the  abbreviation 
(17.3)  Kn  a  (k  +  n)(k  +  n+l), 

these  waves  have  the  res\U.tant  fields 
♦n  =  AnT^sl'  coskx, 

9n  ®  AjiCKjjT“-n(n-l)T“"*]S^'*'*  coskx. 


(17.4) 


*  Garcia,  H.  V.,  1936:  Barotroplc  Waves  in  Straight  Parallel  Slow  with 
Curved  Velocity  Profile.  Tellus,  8,  pp.  82-93* 
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The  factor  In  (l6.8)  which  deteminea  the  phase  velocity  of  the  shear 
wave  component  of  these  waves  has  the  value 

(17.5)  q/(S»*)  *  Kn-n(n-l)/T®. 

This  factor  Is  Independent  of  height  In  the  (nsO)-wave  and  the  (n  =  l)-wave, 
so  the  Instantaneous  relative  phase  speeds  of  these  non-tlltlng  waves  are 
Independent  of  height  If  the  amplitude  ratio  is  the  same  at  all  levels.  They 
have  from  (16.7-8)  the  values 

(17.6)  Q/a  =  l.Sn/r,  C/U  -=  2(l  +  r)/(KoOrKi  ). 

From  (17.^)  these  families  of  waves  have  the  resultant  streamfunctions 

(17.7)  to  =  S^coskx,  ♦,  =  TS^tcoekx, 

with  an  arbitrary  constant  amplitude  factor  omitted. 

!Rie  (n=0)-family  is  symmetric  with  the  same  phase  above  and  below 
the  central  level.  It  is  a  b-state  of  a  symmetric  wave.  The  (n=l)-family 
is  symmetric  with  opposite  phase  above  and  below  the  central  level.  It  is 
an  a-state.  If  the  amplitude  ratio  has  the  value  1.5ui  the  gravity  wave  Is 
stationary  In  these  waves.  If,  furthermore,  Kq  has  the  value  2  3|A  also 
the  shear  wave  Is  stationary  In  the  (n =  0)-famlly.  If  has  this  value 
the  shear  wave  Is  stationary  In  the  (n= l}-famlly.  The  conditions  to  be 
satisfied  by  these  two  families  of  statloneu^  waves  are  accordingly 

(17.8)  stationary  a-wave:  r,  =  1.5»it  Kj  =  p,  (w  =  2  +  3>i). 

(17.9)  stationary  b-wave:  r^  =  1.5tii  K©  =  «•  Kq  =  (k  +  n)(k  +  n  +  1) 

These  families  represent  two  of  the  solutions  of  (16.9)  which  were  found 

by  Qarcla.  He  was  able  to  prove  that  these  are  the  only  stationary  waves 
In  the  Th-shear  layer  (15.I)  If  Its  static  stability  is  weak  (u<^/3). 

For  stronger  static  stability  there  are  other  stationary  waves  (see  sec¬ 
tion  20). 

Let  us  now  compare  the  a-  and  b-waves  in  (17*7)  with  the  correspond¬ 
ing  waves  in  the  constant  shear  layer. 
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l8.  Ih*  (n»l)»  -  This  ««▼«  has  from  (17.?  and  7)  an  lnstan> 

tanaous  a«state  with  ths  fisld 

(18.1)  ♦  ■  Ts'^coskx,  q/(S*t)  «  Ki . 

Its  Intrinsic  upstrsaa  phase  Tsloeitias  for  an  arbitrary  amplitude  ratio  are 

(18.2)  G^,/0  «  1.5»i/r.  ;  C,/n  »  2(1  +  r,  )Ai . 

It  is  a  stationary  neutral  a>wave  if  the  amplitude  ratio  has  the  value  l.?ti( 
and  the  parameter  (see  17.?)  has  the  value 

(18.3)  Kj  =  (k  +  l)(k  +  2)  =  2  +  3u. 

This  condition  is  represented  in  a  k.ii-dlagram  by  a  parabola  throi>gh  the 
origin. 

To  compare  this* parabola  with  Ooldsteln's  a>llne  for  the  stationary 
a>«aves  in  the  constant-shetur  layer,  «e  recall  that  k<B)&kdB)&H,  so  the 
stationary  a-wave  condition  in  (18.3)  may  be  written 

(18.4)  U  =  )4h  +  h*A2.  (H*kd) 

Near  the  origin  this  parabola  coincides  with  Goldstein's  a>line, 

us  k/(1  -  e'*’^)  >  1  to  the  order  of  k’.  For  shorter  waves  the  parabola 
lies  above  the  Goldstein  line,  but  the  two  curves  are  almost  coincident 
for  u<l  (see  fig.  19).  For  example  for  hsI  the  stationary  a-wave  in 
the  constant-shear  layer  has  the  static  stability  u=  (e  - 1)~^ s  0.382, 
while  the  same  wave  in  the  Th-shear  layer  has  the  static  stability 
U=7A2  =  0.583. 

For  other  waves  with  the  amplitude  ratio  1.5u  the  gravity  wave  is 
stationary  while  the  shear  wave  moves  intrinsically  upstream  with  the 
speed  C,  in  (l8.2).  Its  non-dimensional  downstream  speed  with  reference 
to  the  fluid  at  the  central  level  is 

a,  =  (U-C,)/Uo  =  T[1  -  (2  +  3i*)/Ki  ].  (Th-shear) 

Let  us  compare  this  speed  at  the  distance  )&dfrom  the  central  level 
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(Tb0.76~5A)  with  the  a-state  speed  in  the  constant  shear  layer,  namely, 

S',  »  1  -  (ii  + 1)(1  -  e"**)/*.  (Constant -shear) 

The  (Hail)-wave,  for  example,  has  here  the  a-state  speed 

y  s  HHi _ _ .]  B  (Th-shear) 

•  ^  (H+DOi+S)-*  20  •  ^  ^ 

8,  <=  1-  (|4+l)(l-e“^)  =  (Constant-shear) 

For  longer  waves  the  agreement  is  even  closer. 

Ihis  good  agreement  in  th.e  dynamic  a-state  behavior  of  the  symmetric 
wave  (l8.l)  in  the  Th-shear  layer  with  the  symmetric  wave  in  the  constant- 
shear  layer  suggests  that  the  two  layers  have  similar  dynamic  charMteristlcs, 
at  least  for  waves  longer  than  about  six  shear  layer  depths. 

19.  Ihe  b-wave  (naO).  -  This  wave  has  from  (17*5  and  7)  an  instan¬ 
taneous  b-state  with  the  field 

(19.1)  t  =  S^coskx,  q/(S»t)  »  Ko  =  k(k  +  l). 

Its  intrinsic  upstream  phase  velocities  for  an  arbitrary  amplitude  ratio  are 
Q*/U  =  1.5»»/rk ;  C*A  =  2(1  +  r*  )/Ke . 

This  wave  is  stationary  if  its  amplitude  ratio  is  1.5u  and  the  parameter 
Kq  has  the  value 

(19.2)  Ko  =  k(k  +  l)  =  2  +  5»i, 

which  is  represented  in  the  k,|t-diagram  by  a  parabola  identical  to  the 
a-parabola  in  (l8.3)  only  displaced  unit  distance  toward  increasing  k 
(see  fig.  19). 

This  parabola  departs  very  much  from  the  b-line  for  the  constant 
shear  layer.  The  resison  is  that  the  wave  in  (19.1)  has  the  maximum  vorti- 
city  near  the  central  level  where  the  wave  in  the  constant  shear  layer  has 
no  added  vortlcity.  The  parti«Q.  field  from  the  vorticlty  In  the  central 
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reglon  aiigments  the  upstream  propagation  of  the  diear  wave  at  all  levels, 
so  the  statioiwry  b-wave  In  (19*1)  is  shorter  than  the  stationary  b-«ave 
In  the  constant-shear  layer.  If  the  excess  vortlclty  in  the  central 
region  were  removed  from  the  wave,  the  shear  wave  would  move  downstream 
at  all  levels,  similar  to  the  same  wave  in  the  constant  shear  layer. 

It  is  imposaible  to  set  up  a  stationary  b-wave  in  the  Th-shear 
layer  unless  it  is  given  excess  vorticity  near  the  central  level,  where  no 
such  vorticities  can  be  generated  by  kinematic  advection  or  by  dynamic 
action  of  the  overweight.  However,  with  the  aid  of  the  field  of  the  a-wave 
in  (l8.1)  we  may  construct  a  b-wave  which  has  no  excess  vorticity  at  the 
central  level  simply  by  changing  the  sign  of  the  a-wave  vorticity  in  the 
lower  layer.  Accordingly  this  b-wave  has  the  vorticity  field 
(19.3)  q  =  Ki  iTlS^^^coskx. 

Above  the  central  level  this  wave  has  the  same  vorticity  as  the  a-wave  in 
(l8.l).  The  partial  field  from  the  vorticities  below  the  central  level  is 
Laplacean  in  the  upper  region  for  both  waves,  so  the  two  waves  differ  only 
by  a  Laplacean  field.  Above  the  central  level  the  streamfunction  of  the 
vorticity  wave  in  (19.3)  accordingly  has  the  form 

t  =  (TS^  +  Ae"*“)co8kx.  (z>0) 

Since  the  field  is  symmetric  with  respect  to  the  central  level  we  have 
t'(z*0)  =  0,  so  A  =  k”^.  The  streamfunction  of  the  vorticity  waive  (19.3) 
is  therefore 


(19.4)  t  =  (|TlS*'  +  k‘^e"*^l®bco8kx. 

If  this  wave  has  the  amplitude  ratio  1.3|i,  the  gravity  wave  is  stationary, 
and  from  (l6.8)  the  shear  wave  in  the  upper  region  has  the  relative  speed 


_ iLli _  Ci» _ 1 _ 

U  (k  +  l)(k  +  2)  kT(l+T)^ 


■). 


(19.5) 


(r  =  1.5>4) 


Fig.  19  Slotlonorir  womi  in  o  Th-ih«ar  lay*r 
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Its  upstream  speed  is  greater  than  In  the  corresponding  a-wave  (18.2).  For 
very  short  wave  lengths  it  has  the  ssme  asymptotic  behavior  as  the  a-wave. 
In  particular,  the  condition  for  a  stationary  short  wave  is  the  same  (l8.?) 
for  both  waves.  For  the  longer  b-waves  in  (19>^)f  on  the  other  hand,  the 
relative  speed  of  the  shear  wave  is  a  function  of  height,  so  for  no  value 
of  the  static  stability  are  these  waves  stationary  at  all  levels  simultan¬ 
eously.  The  wave  is  stationary  at  a  given  level  zsz(Ig}  if  the  static 
stability  has  the  value  of 

(19.6)  3u  +  2  =  (k  +  l)(k  +  2)/l  + - i - -V,  Cfc(T8)  =  U. 

(  kTgd  +  Ta)*'^ 

At  higher  levels  this  wave  moves  slower  than  the  air,  at  lower  levels 
faster.  The  wave  is  stationary  far  from  the  central  level  if 

3u  +  2=  (k+l)(k  +  2)Cl+ (k2*‘)’^r^.  C»  (Tg*  1)  =  Dq. 

This  line  is  shown  in  the  upper  u,K-diagram  in  fig.  19*  It  meets  the 
Qoldsteln  line  for  the  stationary  b-waves  in  the  constant  shear  layer  at 
H=2kBl.3  and  at  k=4.13.  In  the  spectral  interved  between  these  waves 
the  stationary  b-wave  in  the  constemt  sheaur  layer  has  a  little  more  static 
stability.  For  example  the  (H=2)-wave  in  the  constant  shear  layer  is 
stationary  if  u*  (e®  -  l)/(e*  + 1)  =  0.76,  while  the  same  wave  in  (19.4)  is 
stationary  far  from  the  central  level  if  |is2/3. 

We  recall  from  (13.4)  that  the  Hh-shear  layer  has  maximum  static 
stability  at  the  level  of  T’  c  .  With  this  value  of  Tg  substituted, 
equation  (19.6)  marks  the  waves  in  (19.4)  which  are  stationary  at  the 
level  of  maximum  stability.  The  corresponding  line  is  shown  in  the  upper 
diagram  in  fig.  19*  These  mves  have  considerably  less  static  stability 
than  the  stationary  b-waves  in  the  constant  shear  layer.  The  main  reason 
is  probably  that  the  vorticity  field  in  (19.3)  has  unrealistically  large 
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values  near  the  center  level  where  there  la  no  physical  mechanism  for 
the  generation  of  vorticity. 

Let  us  therefore  consider  a  third  class  of  b-waves  idvlch  has  very 

little  vorticity  near  the  central  level,  namely  waves  with  the  vorticity 
-  k+2 

distribution  q~!IrS  .  The  streamf unction  of  this  wave  is  obtained  by 
taking  the  resultant  of  two  of  the  waves  in  (17«^)*  Omitting  the  trigono¬ 
metric  factor,  the  waves  of  n=0,  and  n=2  are 


qo  -  KqS  , 

♦0  =  s*' 

(2) 

qi  =  (K,T®  -2)5^'^^; 

ts  = 

(Ko) 

The  resultant  of  these,  augmented  by  the  factors  to  the  right,  is 
(19.7)  q  =  ♦  =  (2  +  KqT*  )S^. 


If  the  amplitude  ratio  of  this  wave  is  l.^Ui  Its  gravity  wave  is  station¬ 
ary  and  its  shear  wave  moves  with  the  intrinsic  upstream  speed 


(19.8) 


Si.  -  3u-t-2 

U  (k  +  2)(k  +  5) 


(1  + 


k(k  +  1)1* 


). 


(r=  1.5(ji) 


This  wave  has  maximum  vorticity  at  the  levels  of  T®  =)4(l+)4k)~\  For  the 
very  long  waves  these  coincide  with  the  levels  of  maximum  stability.  For 
k=2  they  coincide  with  the  levels  of  maximum  vorticity  gradient.  For 
k=0.8  the  levels  of  maximum  vorticity  are  in  the  middle  between  the 
levels  of  maximum  stability  and  maximum  vorticity  gradient.  The  wave 
in  (19.7)  is  stationary  at  the  level  of  maximum  stability  (T®=:J4)  if 
the  stability  has  the  value  of 


(19.9)  3m- 2  =  Ck(T«=)4)  =  U 

k"  k  +  4 

In  a  Ufh-diagram  (see  lower  diagram  in  fig.  I9)  this  line  meets  the 
b-line  for  the  constant  shear  layer  at  the  points  (u=0,  HmKqs  1.278...) 
and  (ti=1.75i  h=2,9).  Between  these  points  the  wave  in  (19.9)  has  a 
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llttle  1«SB  statio  stability,  but  otherwiat  ita  propartiaa  are  quite 
aifflilar  to  thoae  of  the  atatioaary  b-wave  in  the  oonatant  ahear  layer. 
Like  the  atationary  a-wavea  in  aeotion  18  thia  ainilarity  probably 
refleeta  a  fundamental  dynamio  ainilarity  between  the  Th-ahear  layer 
and  the  oonatant  ahear  layer.  However,  whereM  the  oonatant  ahear  layer 
haa  only  theae  two  famillea  of  atationary  wavea,  the  !Bi-ahear  layer  has 
an  infinite  number  of  suoh  families. 

20.  The  stationary  wavea  in  the  Pi-ahear  layer.  (Qarcla's  solu¬ 
tions.  )  -  From  (16.9}  the  general  condition  for  a  stationary  wave  in  the 
Th-shear  layer  with  reference  to  the  fluid  at  the  central  level  (Q  =  UsC) 
la 

(20.1)  q  =  -7»t=aS»t.  (i>  =  3|i  +  2) 

When  the  streamfunction  haa  a  sinusoidal  variation  in  the  x-direction, 
Geircla  noticed  that  thia  differential  equation,  by  suitable  transforma¬ 
tion,  is  the  hyper-geometric  differential  equation.  As  such  ita  solu¬ 
tions  may  be  represented  by  the  series 

ra  =  n  . 

(20.2)  t  ^  Aa!I®S*'coskx, 

OB  0 

where  the  coefficients  Ag,  are  the  coefficients  of  the  hyper-geometric 
series.  These  coefficients  are  in  the  present  case  easily  determined 
by  the  formulas  in  (17.4):  The  field  in  (20.2)  has  the  vorticity 

(20.3)  q  =  ^AniCK|BT”-m(m-l)T'°''^3S^'*'^coskXs  ^»A-T“S^'*’ ^  cos  kx. 

mBO  mnO 

In  a  stationary  wave  this  vorticity  must  have  the  value  in  (20.1),  as 
indicated,  The  coefficient  for  T*”  gives 


“  (m  +  l)(n  +  2)A„,^., 
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Rep«at«d  use  of  this  recursion  formula  gives 

(20.4)  A„  =  AqCKq-wXK, (m  even) 

(20.5)  A||j  =  At(Ki -u)(Ka -v)....(K^_,-v)(ml)T^.  (n  odd) 

Thus  i  in  (20.2)  is  a  solution  of  the  differential  equation  if 

(20.6)  u  =  Kjj,  or  3|i  +  2  *  (k  + n)(k  + n  + 1), 

where  n  is  any  non-negative  integer.  If  n  is  an  even  integer,  Aj  =  0  and 
the  coefficients  A^  have  the  values  in  (20.4).  If  n  is  an  odd  integer, 

Ag  =  0  and  the  coefficients  A^  have  the  values  in  (20.5). 

The  stationary  wave  conditions  in  (20.6)  are  represented  in  the 
k,U“diagram  by  a  family  of  identiCBO.  paLrabolas,  each  displaced  unit 
distance  toward  decreasing  k  from  the  next  lower  member.  The  (n=  0)- 
parabola  marks  the  stationary  b -waves  (see  19.2)  which  have  no  nodal 
plane.  The  (n=  l)-parabola  marks  the  stationary  a-waves  (see  l8.5)  which 
have  one  nodal  plane  at  the  central  level.  The  (n= 2)-parabola  marks  a 
family  of  stationary  waves  with  two  nodal  planes  and  so  on.  If  n  is  odd, 
the  central  level  is  a  nodBui  plane  and  the  field  has  opposite  phase  at 
equal  heights  above  and  below.  The  wave  is  of  the  a-type.  If  n  is  even, 
the  wave  has  the  same  phase  and  amplitude  at  equal  heights  above  and  below 
the  interface.  The  wave  is  of  the  b-type. 

By  examining  all  the  solutions  of  the  hyper-geometric  equation 
Qarcia  foimd  that  the  solutions  in  (20.6)  represent  all  the  stationary 
waves  in  the  I3i-shear  layer  in  (15.1)  which  have  bounded  and  continuous 
velocity  fields. 
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